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PREFACE 


The  work  reported  herein  was  conducted  as  part  of  the  Water  Quality 
Research  Program  (WQRP) ,  Work  Unit  32368.  The  WQRP  is  sponsored  by  the  Head¬ 
quarters,  US  Army  Corps  of  Engineers  (HQUSACE) ,  and  is  assigned  to  the  US  Army 
Engineer  Waterways  Experiment  Station  (WES)  under  the  purview  of  the  Environ¬ 
mental  Laboratory  (EL) .  Funding  was  provided  under  Department  of  the  Army 
Appropriation  No.  96X3121,  General  Investigation,  The  WQRP  is  managed  under 
the  Environmental  Resources  Research  and  Assistance  Programs  (ERRAP) , 

Mr.  J.  L.  Decell,  Manager.  Mr.  Robert  C.  Gunkel  was  Assistant  Manager,  ERRAP, 
for  the  WQRP.  Technical  Monitors  during  this  study  were  Messrs.  Pete  Juhle 
and  James  Gottesman  and  Dr.  John  Bushman,  HQUSACE. 

Principal  Investigators  for  the  study  were  Dr.  Mark  S.  Dortch  and 
Ms.  Dorothy  H.  Tillman  of  the  Water  Quality  Modeling  Group  (WQMG) ,  Ecosystem 
Research  and  Simulation  Division  (ERSD),  EL.  This  report  was  prepared  by 
Dr.  Dortch,  Ms.  Tillmcin,  and  Dr.  Barry  Bunch  of  the  WQMG  under  the  direct 
supervision  of  Dr.  Dortch,  Chief,  WQMG,  and  under  the  general  supervision  of 
Mr.  Donald  L,  Robey,  Chief,  ERSD,  and  Dr.  John  Harrison,  Director,  EL.  The 
report  was  reviewed  by  Mr.  Tom  Cole  of  WQMG;  Drs.  Robert  Kennedy  and 
James  Brannon  of  the  Aquatic  Processes  and  Effects  Group,  ERSD;  and 
Dr.  Joe  Nix,  Department  of  Chemistry,  Ouachita  Baptist  University,  Arka- 
delphia,  AR, 

At  the  time  of  publication  of  this  report.  Director  of  WES  was 
Dr.  Robert  W.  Whalin.  Commander  av.d  Deputy  Director  was  COL  Leonard  G. 
Hassell,  EN. 

This  report  should  be  cited  as  follows: 

Dortch,  Mark  S.,  Tillman,  Dorothy  H. ,  and  Bunch,  Barry  W.  1992. 

"Modeling  Water  Quality  of  Reservoir  Tailwaters,"  Technical  Report 

W-92-1,  US  Army  Engineer  Waterways  Experiment  Station,  Vicksburg,  MS. 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non- SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 
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PART  I :  INTRODUCTION 

Background 

1.  Thermal  stratification  in  reservoirs  is  often  accompanied  by  deple¬ 
tion  of  dissolved  oxygen  (DO)  in  the  bottom  waters  (hypolimnion) .  Oxygen 
depletion  and  the  establishment  of  reducing  conditions  in  the  hypolimnion 
increase  mobilization  from  the  sediments  of  dissolved  nutrients  (i.e.,  ammo¬ 
nium  and  inorganic  phosphorus),  sulfide,  reduced  metals  (e.g.,  iron  and  manga¬ 
nese),  and  organic  substances  (i.e.,  simple  organic  acids  and  methane).  These 
substances  can  accumulate  in  the  hypolimnion,  thus  impacting  in-pool  and 
release  water  quality.  Reservoir  releases  that  are  low  in  DO  and  high  in 
reduced  substances  can  threaten  aquatic  life,  cause  water  treatment  problems 
for  downstream  water  supply,  and  can  be  obnoxious  to  downstream  recreational 
users . 

2.  When  water  is  released  to  the  downstream  environment,  stream  reaera¬ 
tion  occurs  and  reduced  substances  begin  to  oxidize  (Stumm  and  Morgan  1981) . 
Water  quality  improves  as  the  water  moves  downstream,  eventually  recovering  to 
a  more  natural  stream  condition.  The  recovery  distance,  which  depends  on 
site-specific  physical,  chemical,  and  biological  processes,  is  often  on  the 
order  of  tens  of  miles  (Nix  et  al.  1991). 

3 .  A  better  understanding  of  the  recovery  mechanisms  and  chemical 
transformations  in  tailwaters  is  needed  to  be  able  to  address  issues  concern¬ 
ing  reseiT'^oir  tailwater  quality  and  to  better  manage  tailwater  quality  prob¬ 
lems.  As  a  result  of  this  need,  a  research  work  unit,  entitled  "Techniques 
for  Evaluating  Water  Quality  of  Reservoir  Tailwaters,"  was  initiated  in  the 
Water  Quality  Research  Program.  Since  surface  releases  generally  contain  DO 
and  are  void  of  troublesome  reduced  substances ,  this  research  focused  on 
degraded  water  quality  downstream  of  dams  with  deep,  anoxic  releases. 

Oblectives 


4.  The  objectives  of  this  research  were  to; 


a.  Develop  an  improved  understanding  of  chemical  transformations  in 
reservoir  tailwaters. 

Provide  guidance  on  sampling  and  analysis  of  tailwater  quality. 

c.  Develop  an  easy-to-use  personal  computer  (PC)  model  to  predict 
water  quality  of  reservoir  tailwaters  given  in-pool  or  release 
concentrations . 

This  report  is  the  second  of  two  reports  produced  from  this  research  work  unit 
and  addresses  Objectives  a  and  c  above.  Objective  b  was  addressed  in  a  report 
by  Nix  et  al.  (1991).  Nix  et  al.  (1991)  present  field  study  results  for  four 
study  sites.  One  of  the  study  sites,  Nimrod  tailwater,  was  revisited  in  1989 
to  gain  additional  information.  Those  data  are  presented  in  Appendix  A. 
Results  from  the  four  field  study  sites  are  analyzed  herein  to  obtain  an 
improved  capability  for  predicting  the  time -dependent  changes  in  reduced  sub¬ 
stances  (e.g.,  reduced  iron  and  manganese)  in  reservoir  tailwaters.  A  numeri¬ 
cal  water  quality  model  for  reservoir  tailwaters  was  developed  and  is  present¬ 
ed  in  Part  IV  and  Appendix  B.  Appendix  C  presents  the  kinetic  rate 
regressions . 


Approach 

5.  Effective  water  quality  management  of  reservoir  tailwaters  requires 
assessment  of  existing  conditions  and  prediction  of  future  conditions  result¬ 
ing  from  structural  and/or  operational  modifications  of  the  dam  or  tailwater 
system.  For  example,  the  conversion  of  a  nonhydropower,  deep  release  to  a  hy¬ 
dropower  release  can  lead  to  water  quality  problems.  Releases  through  hydro- 
power  turbines  are  not  subjected  to  the  usually  high  degree  of  reaeration  to 
which  nonhydropower  releases  are  subjected  (Bohac,  Harshbarger,  and  Lewis 
1983);  thus,  low  DO  may  persist  for  a  number  of  miles  downstream,  which  is 
detrimental  to  downstream  aquatic  habitat. 

6.  Mathematical  water  quality  modeling  is  a  cost-effective  tool  for 
predicting  future  conditions  resulting  from  human  actions  (McCutcheon  1989) . 

A  model  can  be  used  to  estimate  downstream  water  quality  for  a  proposed 
release  condition  at  a  dam.  However,  water  quality  models  are  limited  in  the 
context  of  process  descriptions.  For  example,  consider  a  first-order  loss  or 
decay  of  a  substance  (e.g.,  biochemical  oxygen  demand).  First-order  loss  or 
decay  occurs  when  the  time  rate  of  change  of  a  constituent  concentration  C 
is  directly  proportional  to  the  concentration,  which  is  mathematically  stated 
as 
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dC 

■Jt 


-KC 


(1) 


where  t  is  time.*  Equation  1  has  been  successfully  used  to  describe  the 
kinetics  of  numerous  water  quality  constituents.  The  difficulty  in  applying 
Equation  1  in  practice  is  estimating  or  calibrating  the  reaction  coefficient, 

K  .  The  modeler  usually  tries  to  fit  the  model  to  observed  concentrations  by 
varying  K  during  calibration.  In  reality,  K  can  vary  with  environmental 
conditions,  such  as  temperature,  pH,  DO,  etc.  One  of  the  greatest  difficul¬ 
ties  of  water  quality  modeling  is  that  information  on  the  processes  that 
affect  K  is  limited.  Therefore,  empirical  observations  are  required  to 
develop  K  values  for  site-specific  conditions.  Improved  understanding  of 
the  physicochemical  processes  is  required  for  more  general  application  of 
water  quality  models. 

7.  An  improved  understanding  of  water  quality  processes  can  be  obtained 
through  carefully  designed  field  study  because  laboratory  studies  may  fail  to 
include  all  environmental  factors  that  affect  water  quality  kinetics.  Short¬ 
term  intensive  field  studies  were  conducted  at  four  reservoir  tailwater  sites 
where  physical  and  water  quality  conditions  were  measured  to  quantify  chemical 
transformations.  The  study  sites  were  Little  Missouri  River,  Arkansas,  below 
Lake  Greeson;  Fourche  La  Fave  River,  Arkansas,  below  Nimrod  Reservoir;  Rough 
River,  Kentucky,  below  Rough  River  Reservoir;  and  Guadalupe  River,  Texas, 
below  Canyon  Reservoir.  Narrows  Dam  (i.e..  Lake  Greeson)  has  hydropower; 
Canyon  Dam  was  retrofitted  for  non- Federal  hydropower  shortly  after  this 
study.  The  other  sites  are  nonhydropower.  All  sites  had  deep  releases  char¬ 
acterized  by  the  presence  of  reduced  substances  with  low  DO.  The  study  sites 
provided  a  range  of  physical  and  environmental  conditions,  such  as  varying 
stream  slope,  substrate,  and  pH. 

8.  Samples  were  collected  for  water  quality  analyses  at  a  number  of 
stations  extending  over  the  tailwater  reaches.  The  study  reaches  extended 
from  the  dam  to  about  10  to  36  river  kilometers  downstream,  depending  on  flow 
conditions  and  stream  characteristics,  both  of  which  impact  particle  travel 
time.  Particle  travel  time  is  the  time  required  for  a  parcel  of  water  to  move 


*  For  convenience,  symbols  and  abbreviations  are  listed  in  the  Notation 
(Appendix  D) . 
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from  the  release  point  to  a  given  clstance  downstream.  The  study  reaches  were 
generally  limited  to  a  travel  time  of  about  2  days,  because  most  of  the  chemi¬ 
cal  reactions  occurred  within  this  time  period. 

9.  The  release  flow  rate  from  the  dam  was  held  constant  during  each 
field  study  to  provide  steady  flow  and  nearly  steady- state  conditions  (with 
the  exception  of  dlel  effects) .  This  approach  greatly  simplifies  interpreta¬ 
tion  of  results.  Lagranglan  sampling,  which  tracks  a  water  parcel  as  it  moves 
downstream,  and  snapshot  sampling,  in  which  samples  are  collected  at  all 
stations  at  or  near  the  same  point  in  time,  yielded  similar  results  for  steady 
flows  (Nix  et  al.  1991),  which  means  the  systems  were  approximately  at  steady 
state.  Dye  studies  were  conducted  to  determine  the  travel  time  to  each  sta¬ 
tion.  With  concentrations  measured  at  the  stations  and  with  known  travel  time 
to  the  stations,  it  is  possible  to  evaluate  the  reaction  kinetics  of  Equation 
1  for  various  processes,  such  as  the  loss  of  reduced  manganese.  Part  II  of 
this  report  describes  tailwater  quality  processes  based  upon  literature 
review.  Part  III  presents  the  analyses  of  the  field  studies  for  use  in  model 
development,  and  Part  IV  presents  the  numerical  model  and  its  application  to 
the  field  study  sites. 
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PART  II;  DESCRIPTIONS  OF  TAILWATER  QUALITY  PROCESSES 


Model  Basis 

10.  It  is  useful  to  first  present  the  fundamental  model  that  will  be 
used  to  model  water  quality  of  reservoir  tailwaters.  The  mathematical  model 
is  based  on  the  one - dimens ional  (1-D)  mass  conservation  equation  for  streams 


ac 

It 


(2) 


where 

U  -  stream  mean  velocity 
X  -  distance  downstream 
D  -  longitudinal  dispersion 

S  -  rate  of  change  in  concentration  resulting  from  transformation  or 
chemical  reactions 

Equation  2  assumes  completely  mixed  conditions  over  the  depth  and  width  of  the 
stream.  The  dispersion  term,  which  is  the  first  term  on  the  right  side  of 
Equation  2,  is  usually  much  less  than  advection  (second  term  on  the  left  side 
of  Equation  2)  for  streamflow  (Fischer  et  al.  1979).  Neglecting  dispersion, 
assuming  steady* state  conditions,  and  changing  from  an  Eulerian  to  a 
Lagrangian  reference  frame,  where  U  -  dx/dt  is  the  velocity  of  a  parcel  of 
water.  Equation  2  can  be  simplified  to 


dC 

•at 


±  S 


(3) 


If  S  in  Equation  3  is  a  first-order  loss  rate  (i.e.,  -KC) ,  Equation  3  is 
identical  to  Equation  1. 

11.  The  Eulerian  viewer  observes  concentrations  at  multiple  stream 
stations  as  the  water  flows  by;  the  Lagrangian  viewer  follows  a  parcel  of 
water  as  it  flows  dowtistream  and  observes  its  concentration.  The  analyses  in 
the  following  sections  use  the  Lagrangian  view  because  it  conveniently  relates 
results  from  all  field  study  sites  to  a  common  variable,  travel  time  of  a 
parcel  of  water.  However,  the  numerical  model  discussed  in  Part  IV  uses  the 
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Eulerlan  view,  but  assumes  steady  flow  and  steady-state  conditions.  There¬ 
fore,  the  two  views  are  still  equivalent. 

Manganese 

12.  Anoxic  conditions  in  the  hypolimnion  of  deep  reservoirs  induce  the 

formation  of  dissolved,  reduced  manganese,  Mn*^  or  Mn(II).  Reduced  manganese 
can  accumulate  in  the  reservoir,  affecting  in-pool  as  well  as  release  water 
quality.  Mn*^  accumulates  in  the  water  column  from  sediment  release  and 
reduction  of  Mn(IV)  in  the  water  column  under  anoxic  conditions,  and  is 

mostly  in  the  dissolved  state  with  particulate  Mn'*^^  being  negligible  (Wilson 
1980) . 

13.  Little  work  has  been  conducted  to  evaluate  the  fate  of  reduced 
substances  after  they  are  released  to  the  tailwater  system  (Gordon  1983).  The 
oxidation  of  reduced  manganese  in  the  tailwater  can  be  influenced  by  many 
environmental  factors  such  as  temperature,  pH,  presence  of  DO,  electron 
acceptors,  degree  of  mixing,  etc.  (Gordon  1989).  In  an  effort  to  observe  the 
fate  of  reduced  manganese  spatially  in  the  Duck  River,  Gordon  conducted  a 
field  study  below  Normandy  Dam,  near  Tullahoma,  TN  (Gordon  1989).  Gordon 
found  from  the  Duck  River  field  data  that  the  manganese  removal  rate  was 
first-order  with  time  of  travel,  and  the  reaction  rate  was  much  faster  than 
that  typically  found  in  laboratory  studies.  He  also  found  that  the  removal 
rate  was  a  function  of  river  discharge. 

14.  The  removal  of  reduced  from  reservoir  tailwaters  can  involve 
physical  processes  (i.e.,  adsorption)  as  well  as  oxidation.  Nix  (1986), 

Gordon  (1989),  and  Gordon  and  Burr  (1989)  have  found  that  a  primary  removal 
process  of  can  involve  adsorption  of  Mn*^  onto  an  oxide-coated  substrate. 

Stumm  and  Morgan  (1981)  suggest  that  once  the  Mn*^  is  adsorbed  onto  the  sub¬ 
strate,  oxidation  on  the  manganese  oxide  surface  is  a  slow  process.  Hsiung 
(1987)  reported  that  reduced  manganese  is  removed  primarily  by  oxidation  for 
pH  above  9.0  and  by  surface  catalysis  at  lower  pH.  Wilson  (1980)  also  sug¬ 
gests  that  the  removal  of  Mn^^  during  the  turnover  of  lakes  is  probably 
through  adsorption  onto  particulate  matter  and  co-precipitation  followed  by 
slow  oxidation. 

15.  Work  by  Nix  (1986)  showed  that  gravel  coated  with  hydrous  manganese 
oxides  from  the  tailwater  of  Narrows  Dam  (Lake  Greeson)  effectively  removed 
reduced  manganese  in  the  laboratory,  whereas  removal  without  the  gravel  was 


9 


slow.  Data  from  the  study  sites  report  by  Nix  et  al.  (1991)  indicate  that 
total  and  dissolved  manganese  both  consistently  decrease  in  the  downstream 
direction,  suggesting  an  adsorption  removal  mechanism.  If  only  water  column 
oxidation  was  occurring  without  adsorption  onto  the  bottom,  one  would  think 
that  the  stream  turbulence  would  be  sufficient  to  keep  particulate  (i.e., 
oxidized)  manganese  suspended,  thus  holding  the  total  manganese  concentrations 
relatively  constant  rather  than  decreasing  as  observed. 

16.  Gordon  (1989)  found  through  laboratory  experiments  that  reduced 
manganese  removal  was  associated  with  a  "slime”  found  on  rock  surfaces  in  the 
Duck  River  below  Normandy  Dam.  Gordon  and  Burr  (1989)  were  able  to  effec¬ 
tively  remove  reduced  manganese  using  columns  packed  with  glass  marbles  accli¬ 
mated  with  the  slimy  black  coating  characteristic  of  the  Duck  River  stones. 

The  results  of  these  studies  all  Indicate  the  need  to  consider  the  substrate 
as  a  primary  removal  mechanism  of  Mn*^  in  the  tailwater. 

17.  During  some  earlier  work,  Dortch  and  Hamlin  (1988)  hypothesized 
that  flow  rate  Impacted  the  manganese  removal  rate  in  the  Lake  Greeson  tail- 
water.  The  reasoning  behind  this  hypothesis  was  that  if  adsorption  onto  the 
substrate  is  a  primary  removal  mechanism,  then  altering  the  streamflow  rate 
changes  the  ratio  of  bottom  contact  area  to  cross-sectional  flow  area,  which 
is  equivalent  to  altering  the  hydraulic  depth.  Thus,  higher  flow  rates,  with 
greater  depths ,  should  be  less  effective  in  removing  dissolved  manganese ;  this 
was  found  to  be  true  for  the  Greeson  tailwater. 

18.  Hess,  Byung,  and  Roberts  (1989)  developed  a  manganese  model  using 
the  Duck  River  data  and  subsequently  applied  the  model  to  the  Chattahoochee 
River  below  Buford  Dam,  Georgia.  This  model  contained  two  removal  mechanisms, 
homogeneous  and  heterogeneous  oxidation.  Homogeneous  oxidation  is  direct 
oxidation  of  Mn*^  in  the  water,  and  heterogeneous  removal  involves  catalytic 
activity  (i.e.,  sorption  and  oxidation)  with  available  surfaces  (i.e.,  hydrous 
manganese  oxides  and  suspended  particulates) .  The  Hess  model  for  manganese 
reaction  is  based  on  the  work  of  Morgan  (1967),  which  states 

=  -  (  ko[OH-]2[D01  +  k[0H-]2[D0]  [Mn(IV)]  )  [Mn(II)] 
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where 


[Mn(ll)]  -  dissolved  (reduced)  manganese  concentration,  moles/i 
ko  -  homogeneous  reaction  rate,  day"^  (moles/i)'^ 

[OH"]  -  hydroxide  ion  concentration,  moles/i 
[DO]  -  dissolved  oxygen  concentration,  moles/i 

k  -  heterogeneous  reaction  rate,  day"^  (moles/i)"* 

[Mn(IV)]  -  particulate  (oxidized)  manganese  concentration,  moles/i 

Morgan  (1967)  found  that  the  rate  constants  are  temperature  dependent.  Hess, 
Byung,  and  Roberts  (1989)  reported  that  pH  has  a  much  greater  effect  on  the 
manganese  oxidation  rate  than  dissolved  oxygen  concentration. 

19.  The  major  limitation  of  the  Hess  model  is  that  it  does  not  include 
sorption  of  manganese  onto  the  channel  bottom  (substrate).  The  ideas  devel¬ 
oped  in  Equation  4  provide  a  good  starting  point  for  analyzing  the  data 
obtained  within  the  present  study,  provided  that  the  effects  of  substrate 
adsorption  are  included. 


Iron 

20.  The  rate  of  oxidation  of  reduced  iron  (i.e.,  ferrous  iron,  Fe"'"'^, 
Fe(II))  is  first-order  with  respect  to  concentrations  of  Fe(II)  and  O2  and 
inverse  second- order  with  respect  to  the  iT  concentration  (i.e.  ,  [H'*']  in 
moles/i  where  [H"*"]  -  10"^®)  in  solutions  of  pH  >  5  (Stumm  and  Morgan  1981). 
This  rate  law  is  expressed  as 


<l[Fe(II)]  ,  ,  [Fe(Il)][D0] 

[ipp 


(5) 


where  [Fe(II)]  is  the  concentration  in  moles/i  of  ferrous  iron  and 
kp,  -  3.0  X  10"^^  moles/i  min"^  at  20  "C.  Equation  5  is  in  the  form  of 
Equation  1  with  the  overall  reaction  rate  for  ferrous  iron  described  as 


Kp,  =  kp. 


[DO] 

Tirp 


(6) 
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The  reaction  rate  Increases  with  temperature;  using  an  activation  energy  of 
23.0  kcal/mole  in  the  Arrhenius  equation  (Stumm  and  Morgan  1981),  the  tempera 
ture  dependence  can  be  described  with 


Kj  -  Kao  (1.14)"-“> 


(7) 


where 

Ki  -  reaction  rate  at  temperature  T,  *0 

K20  -  reaction  rate  at  20  "C 

Sxmg  and  Morgan  (1980)  concluded  that  temperature  dependence  is  primarily  a 
result  of  the  change  in  the  tT*'  concentration  due  to  the  temperature  dependence 
of  the  Ionization  constant  of  water. 

21.  Sung  and  Morgan  (1980)  also  reported  that  the  rate  constant  (l.e., 
kp,  of  Eqtiatlon  5)  Is  Inversely  proportional  to  Ionic  strength  and  Individual 
anions  (e.g.,  Cl"  and  SO*"^) .  For  example,  in  the  pH  range  of  6.5  to  7.2,  sul¬ 
fate  Ions  have  a  significant  retarding  influence  on  the  oxygenation  rate.  For 
pH  of  about  7  or  greater,  data  Indicate  an  autocatalytic  rate  expression, 
l.e.,  homogeneous  and  heterogeneous  oxidations  with  the  latter  Involving 
adsorption  of  ferrous  iron  by  ferric  hydroxide  (Sung  and  Morgan  1980) . 

Sarlkaya  (1990)  proposed  recycling  of  ferric  sludge  at  treatment  plants  to 
provide  effective  catalytic  removal  of  ferrous  Iron,  thus  reducing  volumes 
required  for  aeration  tanks. 

22.  As  reducing  conditions  progress  In  anoxic  reservoir  bottom  sedi¬ 
ments,  sulfate  reduction  follows  iron  reduction,  and  Iron  sulfide  forms  (Stiiimn 
and  Morgan  1981,  Wetzel  1975).  Although  iron  sulfide  Is  Insoluble  at  the  pH 
conditions  typically  found  In  reservoir  waters,  iron  sulfide  may  be  temporari¬ 
ly  suspended  In  the  water  column  and  released  to  the  tallwater. 

23.  The  presence  of  Iron  sulfide  In  the  tallwater  can  complicate  the 
understanding  of  ferrous  Iron  and  sulfide  reactions.  Generally,  not  enough 
H2S  Is  available  for  pyrlte  (FeS2)  formation;  thus.  Iron  sulfide  should  exist 
as  FeS  in  fresh  water  (Berner  1980) .  The  oxidation  of  FeS  is  rapid.  Connell 
(1966)  showed  that  120  ppm  total  sulfides  In  a  soil  was  reduced  to  5  ppm  after 
2  hr  of  exposure  to  air;  this  translates  Into  an  oxidation  rate  of  about  38.0 
day"^ . 

24.  The  reactions  for  FeS2  should  be  somewhat  similar  to  those  that 
occur  In  streams  subjected  to  acid  mine  drainage  as  explained  by  Stumm  and 
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Morgan  (1981).  The  sulfide  of  FeS2  Is  oxidized  to  sulfate-releasing  dissolved 
ferrous  iron  and  acidity  (i.e.,  H''^) .  Subsequently,  ferrous  iron  is  oxidized 
to  ferric  iron,  which  hydrolyzes  to  form  insoluble  ferric  hydroxide,  releasing 
more  acidity.  Ferric  iron  can  be  rapidly  reduced  by  FeS2  to  oxidize  sulfide, 
releasing  more  dissolved  ferrous  iron  and  acidity. 

25.  Because  iron  sulfide  reactions  in  acid  mine  drainage  usually  occur 
with  very  low  pH,  one  would  expect  the  oxidation  of  ferrous  iron  to  ferric 
iron  to  be  slow.  In  fact,  the  rate  of  oxidation  of  ferrous  iron  for  low  pH 
(i.e.,  pH  below  4  to  5)  is  controlled  primarily  by  the  action  of  autotrophic 
iron- oxidizing  bacteria  rather  than  chemical  oxidation  (Stumm  and  Morgan  1981; 
Noike,  Kanji,  and  Jun' Ichiro  1983).  However,  Nolke,  Kanji,  and  Jun' Ichiro 
(1983)  found  ferrous  iron  oxidation  rates  between  25.9  and  49.0  day*^  in  a 
low-pH  stream  (pH  1.5  to  2.9)  receiving  acid  mine  drainage.  These  rates  are 
comparable  to  those  that  occur  at  the  higher  pH  values  found  in  natural 
waters.  The  oxidation  of  sulfide  by  oxygen  or  hydroxides  also  occurs  rapidly. 
Stumm  and  Morgan  (1981)  report  half-times  on  the  order  of  20  to  1,000  min  for 
oxidation  of  iron  sulfide  by  ferric  hydroxide;  this  corresponds  to  reaction 
rates  of  50  day'^  to  1  day*^. 


Svlf.idg 

26.  Dissolved  sulfide,  S[-1I],  can  exist  as  hydrogen  sulfide  (H2S) , 
bisulfide  ion  (HS*^,  and  S**.  Negligible  amounts  (less  than  0.5  percent)  of 
dissolved  sulfide  exist  as  for  the  pH  range  of  most  streams  and  lakes;  the 
distribution  between  HS-  and  H2S  varies  with  pH,  with  high  pH  favoring  HS‘  and 
low  pH  favoring  H2S  (American  Public  Health  Association  et  al.  1981).  At  a  pH 
of  7.0,  sulfide  exists  as  about  half  hydrogen  sulfide  and  half  bisulfide  ion. 
Hydrogen  sulfide  can  be  removed  from  the  water  through  volatilization  and 
oxidation,  although  oxidation  of  H2S  occurs  very  slowly  compared  with  oxida¬ 
tion  of  HS*  (Millero  1986).  HS"  is  removed  primarily  by  oxidation. 

27 .  The  general  consensus  from  the  literature  is  that  oxidation 
kinetics  of  sulfide  in  natural  waters  is  very  complex  and  poorly  understood 
(Chen  and  Morris  1972a,  Stumm  and  Morgan  1981,  Millero  1986).  Based  upon  the 
review  by  Millero  (1986),  the  oxidation  of  sulfide  can  be  described  as  first- 
order  with  respect  to  HS*  with  a  rate  that  increases  with  increasing  tempera¬ 
ture,  pH,  and  concentrations  of  DO,  metal  ions,  and  initial  sulfide.  There  is 
also  general  agreement  that  sulfur  bacteria  play  an  important  role  in  sulfide 
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oxidation  in  natural  and  waste  waters.  These  Influencing  factors  help  to 
explain  why  the  first-order  sulfide  oxidation  rates  reported  vary  between 
0.26  and  about  55.0  day*^,  which  correspond  to  half-lives  between  65.0  and 
0.3  hr  (Millero  1986). 

28.  The  oxidation  rate  of  sulfide  is  strongly  dependent  on  pH.  For  an 
increase  in  pH  from  about  6.0  to  8.0,  there  is  an  eightfold  increase  in  the 
rate  (Chen  and  Morris  1972b) .  Millero  (1986)  attributes  this  to  the  shift  in 
sulfide  concentration  from  predominantly  nonreactive  H2S  to  predominantly 
reactive  HS'  with  increasing  pH.  This  reasoning  suggests  that  HS*  concentra¬ 
tion,  based  on  pH,  should  be  used  in  the  kinetic  equation. 

29.  Almgren  and  Hagstrom  (1974)  found  a  second-order  rate  dependence 
with  the  oxidation  dependent  on  [S(-II)]  and  [O2] .  However,  Wilmot  et  al. 
(1988)  reported  that  the  order  of  reaction  with  respect  to  [O2]  was  only  0.2 
rather  than  1.0;  they  preferred  to  use  the  Michaelis -Menton  law  to  describe 
the  influence  of  DO,  which  is  stated  as 


dS  ^  DO 

•at  ^  TXii,;  +  T50 


(8) 


where 

S  ••  sulfide  concentration  in  the  form  of  HS*,  mg/i 
k,  -  sulfide  oxidation  rate,  day*^ 

DO  -  dissolved  oxygen  concentration,  mg/i 
DOi/2  ~  dissolved  oxygen  half- saturation  constant,  mg/i 
A  value  of  3.0  was  recommended  by  Wilmot  et  al.  (1988)  for  DO1/2  •  O'Brien 
and  Birkner  (1977)  reported  that  the  order  of  reaction  with  respect  to  [O2] 
was  about  0.8.  Chen  and  Morris  (1972a)  also  found  that  the  variation  of  rate 
with  DO  was  less  than  first-order  (e.g.,  order  of  about  0.5).  These  discrep¬ 
ancies  are  probably  related  to  the  different  experimental  conditions  and 
methods  used  (Millero  1986).  In  the  absence  of  better  information.  Equation  8 
is  appealing  from  a  modeling  standpoint. 

30.  Ions  present  in  natural  waters  can  enhance  sulfide  oxidation  rates. 
Almgren  and  Hagstrom  (1974)  found  fairly  high  rates  of  sulfide  oxidation  for 
seawater  (k,  on  the  order  of  25.0  day'^) .  O'Brien  and  Birkner  (1977)  reported 
that  sulfide  oxidation  increases  with  increasing  ionic  strength.  Paschke, 
Hwang,  and  Johnson  (1977)  found  that  the  use  of  manganese  salts  as  catalysts 
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and  adjusting  the  pH  to  8.5  were  effective  In  removing  approximately 
90  percent  of  the  sulfide  within  1  to  2  hr  for  wastewater  from  a  cellulose 
sponge -making  operation.  This  result  translates  Into  sulfide  oxidation  rates 
between  about  27  and  55  day"^.  Chen  and  Morris  (1972a)  also  found  strong  cat¬ 
alytic  effects  of  metal  Ions,  e.g.,  orders  of  magnitude  reduction  In  oxidation 
rates . 

31.  Sulfide -oxidizing  bacteria  are  also  considered  to  Increase  the 
oxidation  rate  (Chen  and  Morris  1972a,  Ullmot  et  al.  1988).  Wllmot  et  al. 
(1988)  reported  oxidation  rates  between  about  14  and  130  day'^  for  wastewaters 
with  as  much  as  half  of  the  oxidation  attributed  to  bacteria.  Wllmot  et  al. 
(1988)  also  determined  that  larger  Initial  sulfide  concentrations  resulted  In 
greater  oxidation  rates. 

32.  The  sulfide  oxidation  rate  approximately  doubles  for  a  15.0  ”C  tem¬ 
perature  Increase  (Ullmot  et  al.  1988).  This  dependency  on  temperature  can  be 
represented  as 


(k,)^  =  (k,)^^  (1.05)<T-20)  (9) 

where  (^*^1)20  oxidation  rate  at  20  *C. 

33.  Hydrogen  sulfide  (H2S)  Is  a  highly  volatile  dissolved  gas  that  oxi¬ 
dizes  very  slowly  (Chen  and  Morris  1972a,  1972b).  Therefore,  the  major  loss 
mechanism  Is  through  volatilization  from  the  water.  Methods  described  by 
Lyman  et  al.  (1982)  were  used  for  determining  volatilization  rates  of  H2S. 
Henry's  constant  at  25  *C  for  H2S  Is  1.01  x  10’^  atm  mVniol.  Therefore,  the 
transfer  Is  liquid  phase  controlled,  and  volatilization  can  be  computed  with  a 
first-order  rate  constant  from  Information  on  the  dissolved  oxygen  reaeratlon 
rate  In  streams. 


K,  =  K, 


D" 

D" 


(10) 


where 

K,  -  volatilization  rate  of  H2S,  day'^ 

K2  -  DO  stream  reaeratlon  rate,  day'^ 

D"  -  molecular  dlffuslvlty  of  H2S,  craVs®^ 


15 


D*  -  molecular  diffusivity  of  DO,  cmVsec 
-  molecular  weight  of  DO,  g/mole 
M*  -  molecular  weight  of  H2S,  g/mole 

Since  the  molecular  weights  of  hydrogen  sulfide  and  DO  are  about  equal,  K, 
is  approximately  equal  to  K2  .  With  little  or  no  H2S  in  the  atmosphere,  the 
volatilization  (mg/i/day)  is  computed  as 


J  -  -KjC. 


(11) 


where  C,  is  the  concentration  of  H2S  in  the  stream.  The  stream  reaeratlon 
rate  can  be  estimated  from  a  variety  of  mathematical  formulations  as  noted  in 
Appendix  B. 


Other  Water  Quality  Constituents 

34.  Other  water  quality  constituents  that  are  of  primary  interest  in 
reservoir  tallwaters  are  temperature,  dissolved  oxygen,  and  nutrients.  Heat 
exchange  mechanisms  are  relatively  well  understood,  and  the  capability  to 
accurately  model  steam  temperature  is  well  established  (Dortch  and  Martin 
1989).  Therefore,  there  was  no  emphasis  in  this  research  on  predicting  steam 
temperature.  Dissolved  oxygen  is  required  for  aquatic  life;  thus,  it  is  of 
prime  importance  in  reservoir  tallwaters,  and  the  sources  and  sinks  of  oxygen 
must  be  considered.  The  nutrients  nitrogen  and  phosphorus  can  exist  in  vari¬ 
ous  forms:  organic,  inorganic,  particulate,  and  dissolved.  The  only  nutrient 
form  that  presents  an  Immediate  potential  problem  *  tallwaters  dominated  by 
anoxic  hypo limnetic  releases  is  ammonia.  Ammonia  ..s  mportant  because  oxygen 
is  used  during  nitrification  to  nitrate.  The  photosynthesis,  respiration,  and 
nutrient  uptake  of  phytoplankton  and  macrophytes  were  not  considered  in  this 
work  because  of  the  relatively  minor  influence  on  tailwater  quality  processes 
at  the  sites  studied  (Nix  et  al.  1991). 

35.  Oxidation  of  reduced  substances  in  the  tailwater  is  a  major  sink, 
or  loss  mechanism,  of  oxygen  (Nix  et  al.  1991).  The  primary  source  of  oxygen 
is  stream  reaeration.  The  reactions  and  their  products  must  be  known  in  order 
to  determine  how  much  oxygen  is  taken  up  during  oxidation.  Once  this  is 
known,  the  stoichiometry  for  the  amount  of  DO  taken  up  per  unit  of  substance 
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oxidized  can  be  specified  within  the  model.  Each  of  the  reduced  substances 
discussed  above,  plus  ammonia,  Is  considered  below. 

36.  The  oxidation  of  Hn(II)  may  be  represented  (Benefield,  Judkins,  and 
Weand  1982)  by 


2Mn**  Oj  ♦  2HaO  «>  2MnOa  ♦  4H*  (12) 

This  reaction  requires  0.29  mg/i  of  DO  to  oxidize  1.0  mg/i  of  Mn.  However, 
Stumm  and  Morgan  (1981)  point  out  that  experimental  findings  Indicate  the 
extent  of  Mn(II)  removal  is  not  accounted  for  by  the  stoichiometry  of  the 
oxidation  reaction  alone,  and  the  products  of  oxygenation  are  nonstolchlo- 
metric,  showing  various  degrees  of  oxidation  ranging  from  about  30  to  90  per¬ 
cent  oxidation  to  Mn02.  Thirty  percent  oxidation  requires  about  0.10  mg/i  of 
DO  per  1.0  mg/i  Mn.  A  DO  uptake  stoichiometry  of  0.29  mg/i  per  mg/i  of  Mn 
removed  can  be  used  as  a  worst  case  for  modeling. 

37.  The  oxidation  of  Fe(II)  may  be  represented  (Stumm  and  Morgan  1981; 
Benefield,  Judkins,  and  Weand  1982)  by 

4Fe*2  +  O,  +  lOHjO  *>  4Fe(OH)3  ♦  8H*  (13) 


This  reaction  requires  0.14  mg/i  of  DO  to  oxidize  1.0  mg/i  of  Fe(II).  Based 
on  stoichiometry,  the  total  oxidation  of  FeS  to  ferric  iron  and  sulfate 
requires  about  0.73  mg/i  of  DO  per  1.0  mg/i  of  FeS,  and  FeS2  requires  about 
1.0  mg/i  of  DO  per  1.0  mg/i  of  FeS2. 

38.  Oxidation  of  sulfide  can  produce  sulfur  (S) ,  thiosulfate  (8203'^), 
sulfite  (SOa"^),  or  sulfate  (SO^"^)  according  to  the  following  reactions  (Chen 
and  Morris  1972b) : 


HS'  +  1/20,  +  H*  =>  H,0  S  (l'^) 

2HS-  ♦  20,  +  *>  H,0  ♦  8,03"*  (15) 
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HS-  ♦  3/20a  «>  SOj*  +  H 


(16) 


HS-  20j  ->  so;®  +  H* 


(17) 


The  above  reactions  (l.e.,  Equations  14-17)  require  0.5,  1.0,  1.5,  and 
2.0  mg/i  DO  per  1.0  mg/i  S“^,  or  0.48,  0.97,  1.45,  and  1.94  mg/i  DO  per 
1.0  mg/i  HS",  respectively.  Equation  16  Is  an  Intermediate  step  for  sulfate 
production,  where  another  half  molecule  of  O2  Is  required  to  complete  the 
reaction.  However,  the  oxidation  of  sulfite  to  sulfate  has  been  observed  to 
occur  slowly  (O'Brien  and  Blrkner  1977,  Chen  and  Morris  1972b).  Thiosulfate 
Is  considered  a  relatively  stable  oxidation  product  (O'Brien  and  Blrkner  1977) 
that  Is  slowly  oxidized  to  sulfate  with  an  observed  half  time  of  8  days 
(Avrahaml  and  Golding  1968).  At  low  pH,  elemental  sulfur  has  been  observed  as 
a  precipitate  (l.e..  Equation  14),  whereas  for  high  pH,  thiosulfate  (Equa¬ 
tion  15)  Is  the  primary  product  (Chen  and  Morris  1972a).  At  neutral  pH,  a 
high  ratio  of  sulfide  to  oxygen  results  In  precipitation  of  sulfur,  while  a 
low  ratio  results  In  direct  oxidation  to  thiosulfate  (Chen  and  Morris  1972b) . 

39.  A  literature  review  by  O'Brien  and  Blrkner  (1977)  Indicated  that 
most  experiments  observed  sulfite,  thiosulfate,  and  sulfate  as  oxidation  pro¬ 
ducts  of  sulfide.  They  proposed  the  hypothetical  reaction  model 

4HS-  5  1/2  Oj  ->  SaOj"®  ♦  SOj'®  ♦  SO*"®  +  2H*  +  HjO  (18) 


This  reaction  requires  1.38  mg/i  O2  per  1.0  mg/i  S"®,  or  1.33  mg/i  O2  per 
1.0  mg/i  HS".  Experimental  results  of  O'Brien  and  Blrkner  (1977)  Indicated  an 
average  stoichiometric  value  of  1.36  mg/i  O2  per  1.0  mg/i  S"®  at  pH  -  7.55. 
Resplrometrlc  studies  by  Chen  and  Morris  (1972a)  showed  that  about  one  mole¬ 
cule  of  oxygen  was  required  per  atom  of  sulfur  oxidized.  Indicating  a  net 
oxidation  primarily  to  thiosulfate  for  a  pH  range  of  6.7  to  13.5.  In  view  of 
the  variation  In  results  and  the  complications  of  sulfide  reactions,  the  value 
of  1.38  mg/i  O2  per  1.0  mg/i  S"®  Is  recommended  to  provide  a  reasonable  yet 
conservative  estimate  of  DO  uptake  requirements. 

40.  Ammonia  can  exist  In  un- Ionized  form,  NHs,  or  In  Ionized  form,  NH^'*^ 
(ammonium) .  Models  usually  consider  total  amnonia  and  do  not  distinguish  the 
two  forms.  The  distribution  between  Ionized  and  un- Ionized  ammonia 


18 


concentration  varies  with  temperature  and  pH,  where  un- ionized  ammonia 
increases  with  temperature  and  pH.  However,  un- ionized  ammonia  usually  exists 
in  small  quantities  in  most  natural  waters.  For  example,  for  a  pH  of  8.0  and 
a  temperature  of  20.0  "C,  less  than  4.0  percent  of  the  total  ammonia  is  NH3 
(Bowie  et  al.  1985). 

41.  Nitrifying  bacteria  convert  ammonivun  (NH*'*’)  to  nitrate  (NOs")  where 
the  overall  reaction  is  described  as  follows  (Wetzel  1975) 

% 

NH4  ♦  20j  ->  NOj*  ♦  HjO  ♦  2H* 

thus  requiring  4.57  mg/i  O2  per  1.0  mg/i  anuaonium  nitrogen  (i.e.,  NH*"*^  as  N) . 
Nitrite  (N02')  is  an  intermediate  product  that  is  highly  labile  to  oxidation. 
Field  studies  of  reservoir  tallwaters  (Nix  et  al.  1991)  Indicate  that  almost 
all  of  the  loss  in  ammonia  can  be  accounted  for  by  nitrification  to  nitrate, 
i.e.,  loss  in  ammonia  nitrogen  approximately  equaled  gain  in  nitrate  nitrogen 
with  some  loss  of  nitrate  due  to  denitrification. 

42.  The  rate  of  nitrification,  or  the  oxidation  of  ammonium,  can  vary 
widely  and  is  considered  to  be  first-order  dependent  on  ammonium  concentration 
and  dependent  on  temperature,  DO,  pH,  and  stream  hydraulics,  such  as  wetted 
perimeter  and  depth  (Bowie  et  al.  1985).  Measured  nitrification  rates  between 
0.1  and  0.5  day'^  are  typical  for  the  deeper,  larger  water  bodies,  and  rates 
greater  than  1.0  day'^  are  not  uncommon  for  shallow  streams  (Thomann  and 
Mueller  1987).  A  temperature  correction  is  usually  used  in  the  form  of  Equa¬ 
tion  7  with  a  value  of  1.08  for  the  base  of  the  power  term.  Much  work  remains 
to  be  done  before  there  exists  a  capability  to  predict  the  nitrification  rate. 

43.  Nitrification  is  a  source  for  nitrate  nitrogen,  and  denitrification 
and  algal  uptake  are  sinks.  Most  stream  water  quality  models,  such  as  QUAL2E 
(Brown  and  Barnwell  1987),  do  not  provide  for  denitrification  when  the  stream 
is  aerated  (e.g.,  DO  greater  than  about  0.5  to  2.0  mg/i).  However,  it  is  well 
established  that  denitrification  can  still  occur  in  the  bottom  sediments  of 
well -oxygenated  streams,  thus  removing  nitrate  from  the  water  column  (Hill 
1979) .  If  proper  representation  of  the  nitrate  balance  is  desired,  then  the 
model  should  Include  nitrification  and  denitrification,  regardless  of  the  DO 
concentration  of  the  water  column. 

44.  Sediment  oxygen  demand  (SOD)  can  have  a  significant  Impact  on  DO  in 
surface  waters  (Hatcher  1986) .  Benthic  respirometers  were  used  at  several  of 
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the  tailwater  study  sites  in  this  research.  SOD  values  of  about  1.0  to  2.0 
g/mVday  were  measured  in  the  Nimrod  and  Rough  River  tailwaters  where  the  bed 
consisted  predominantly  of  fine  sediments  (Nix  et  al.  1991).  These  rates 
exert  an  oxygen  demand  that  is  comparable  to  the  demands  caused  by  the  oxida¬ 
tion  of  reduced  substances  and  on  the  same  order  of  magnitude  as  reaeration 
for  sluggish  streams.  Therefore,  SOD  should  be  considered  for  a  tailwater  DO 
model . 

45.  Most  surface  waters  experience  oxygen  uptake  due  to  the  decomposi¬ 
tion  of  organic  matter  (i.e.,  organic  carbon),  in  addition  to  the  other  oxygen 
sinks  discussed  above.  There  was  an  attempt  to  quantify  this  demand  during 
the  tailwater  field  studies  of  Nix  et  al.  (1991).  Standard,  seeded  5-day  bio¬ 
chemical  oxygen  demand  (BOD)  measurements  Indicated  that  BOD  in  the  tailwaters 
ranged  from  0.8  to  greater  than  7.0  mg/i.  However,  this  measurement  includes 
any  nitrogen  demand  and/or  chemical  oxidation,  such  as  oxidation  of  iron  and 
iron  sulfide.  Therefore,  these  measurements  were  of  little  use  for  determin¬ 
ing  organic  carbon  oxygen  demand. 

46.  Ten-day,  unseeded  BOD  studies  were  also  conducted  where  multiple 
bottles  were  analyzed  at  1-  to  2 -day  intervals  for  each  sample  location  (Nix 
et  al.  1991).  During  these  studies,  various  constituents  that  can  exert  an 
oxygen  demand  (ammonium,  iron,  etc.)  were  measured,  in  addition  to  DO  and 
total  organic  carbon  (TOC) .  The  hope  was  to  stoichiometrlcally  account  for 
the  various  components  of  known  oxygen  demand  and  calculate  the  organic  carbon 
demand  by  difference.  Although  there  was  a  definite  oxygen  demand  in  each  10- 
day  BOD  study,  the  TOC  did  not  show  any  trend,  and  the  studies  failed  to  re¬ 
veal  any  definitive  quantification  of  the  organic  carbon  demand. 

47.  The  QUAL2E  model  has  carbonaceous  BOD  (i.e.,  CBOD)  in  it  as  a 
modeled  state  variable.  This  variable  can  be  used  to  represent  the  organic 
carbon  demand  since  the  other  oxygen  demands  will  be  accounted  for  directly 
through  modeled  state  variables.  The  problem  remains  of  how  to  estimate  a 
value  for  CBOD  released  from  the  reservoir  that  properly  represents  the 
organic  carbon  oxygen  demand. 

48.  Reaeration  is  considered  to  be  the  primary  source  of  DO  in  reser¬ 
voir  tailwaters  fed  by  deep  releases  because  production  of  oxygen  through 
photosynthesis  decreases  with  temperature,  and  deep  releases  tend  to  be  cold. 
Additionally,  deep  releases  tend  to  be  void  of  phytoplankton.  Light  and  dark 
bottle  DO  studies  indicated  that  this  is  a  good  assumption  (Nix  et  al.  1991). 
Numerous  stream  reaeration  models  have  been  developed  (Bowie  et  al.  1985),  and 
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most  stream  water  quality  models  allow  the  user  to  select  from  a  number  of 
these  models.  Although  this  area  of  stream  water  quality  modeling  Is  well 
developed,  there  Is  still  not  a  unified  theory  or  equation  to  predict  reaera- 
tlon  for  a  variety  of  stream  conditions.  It  was  beyond  the  scope  of  this 
research  to  Investigate  Improved  methods  for  modeling  steam  reaeratlon. 
Existing  formulae  are  used  for  the  model. 
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PART  III:  DATA  ANALYSES 


49.  A  list  of  water  quality  constituents  collected  during  the  field 
studies  of  Nix  et  al.  (1991)  Is  provided  In  Table  1.  The  analyses  discussed 
herein  do  not  address  each  constituent  but  are  limited  to  those  of  Interest 
for  the  water  quality  model.  Presently,  the  model  focuses  on  DO,  the  vari¬ 
ables  that  Impact  DO,  and  the  most  problematic  substances  (l.e.,  reduced  Iron, 
reduced  manganese,  and  sulfide).  Eventually,  the  model  may  be  expanded  to 
Include  other  processes,  such  as  pH  dynamics,  and  the  analysis  of  other  data 
collected  during  the  field  studies  may  become  necessary.  The  sections  that 
follow  In  Part  III  describe  analyses  of  field  data  to  support  development  of 
the  kinetic  algorithms  of  the  numerical  model  discussed  In  Part  IV. 


Table  1 

Water  Quality  Constituents  Measured 
During  Field  Studies 


Temperature 
Dissolved  oxygen 
Specific  conductance 
pH 

Free  carbon  dioxide 

Alkalinity 

Total  organic  carbon 

Carbonaceous  biochemical  oxygen  demand 
Iron  (total  and  dissolved) 

Manganese  (total  and  dissolved) 

Sulfate 

Sulfide 

Nitrate -nitrogen 

Ammonia - nl trogen 

Total  KJeldahl  nitrogen 

Total  phosphorus 

Soluble  reactive  phosphorus 

Chloride 

Turbidity 


Reduced  Manganese  Removal 


SO.  Data  taken  from  the  four  tallwater  study  sites  (Nix  et  al.  1991; 
Appendix  A)  plus  data  from  the  earlier  Greeson  tallwater  study  (Nix  1986)  and 
Duck  River  (Gordon  1989)  were  analyzed  for  dissolved  manganese  remova"*  rate. 
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A  general  equation  to  calculate  the  first-order  removal  rate  (K^^)  for  Hn^^  was 
developed  based  on  stream  characteristics.  The  procedure  for  development  of 
an  equation  to  predict  was  accomplished  in  several  steps. 

51.  The  first  step  was  to  determine  the  observed  dissolved  manganese 
loss  rate  from  the  field  data  of  dissolved  manganese  and  travel  time. 

The  SAS  procedure  PROC  REG  (SAS  Institute,  Inc.  1988)  was  used  to  obtain  the 
least  squares  best  fit  of  the  natural  logarithm  transform  of  the  observed 
concentrations  versus  time.  This  regression  equation  form,  which  is  shown  in 
Equation  20,  is  the  natural  logarithm  of  Equation  1,  where  the  slope  of  the 
line  is  ,  or 


In  C  *  In  Cq  —  b 


(20) 


where 

C-  Mn*^  observed  at  a  downstream  station  with  travel  time,  t,  mg/t 
Co  -  initial  (release)  Mn'*’^  concentration,  mg/i 

-  Mn«  loss  rate  for  a  specific  field  condition,  1/day 
t-  travel  time  to  station,  days 

The  regression  models  developed  from  each  observed  data  set  to  predict  concen¬ 
tration  of  Mn*^  versus  time  (including  regression  coefficients,  equation  sig¬ 
nificance  values,  and  correlation  coefficients)  are  listed  in  Table  2.  The 
initial  concentration  Cg  was  a  free  parameter  to  be  estimated  (i.e.,  the 
intercept  of  the  regression).  If  the  significance  probability  of  the  equation 
was  greater  than  0.05,  an  estimated  was  not  accepted.  Thus,  unaccept¬ 
able  values  were  not  used  further  in  the  determination  of  a  general 

equation  for  .  All  the  data  sets  (i.e.,  observed  dissolved  manganese 

concentration  versus  time)  and  the  regression  lines  are  plotted  in  Appendix  C. 
It  should  be  noted  that  the  manganese  data  for  Rough  River  wer<;  separated  and 
analyzed  as  upper  (i.e.,  the  Falls  of  Rough  pool  caused  by  the  mill  dam)  and 
lower  (i.e.,  free-flowing  stream)  reaches  since  the  data  showed  different 
slopes  with  a  change  in  slope  at  the  Falls  of  Rough  Dam. 

52.  The  next  step  was  to  try  to  relate  the  observed  (through  regres¬ 
sion)  values  to  other  variables,  such  as  stream-average  ambient  water 

quality  conditions.  The  initial  approach  used  stepwise,  multiple  linear 
regression  analysis  (FROG  STEPWISE,  SAS  Institute,  Inc.  1988)  to  identify 
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Table  2 


Regression  Results  for  Concentration  Versus  Time  from  Observed  Data 


Data  Set 

Independent 

Regression 

Equation 

_£i- 

Duck  River 

1 

Intercept 

0.45 

0.0001 

0.98 

Time 

-1.19 

Duck  River 

2 

Intercept 

0.44 

0.0001 

0.96 

Time 

-1.65 

Duck  River 

3 

Intercept 

0.63 

0.0001 

0.96 

Time 

-1.72 

Little  Missouri 

1 

Intercept 

-1.38 

0.2810 

0.36 

River 

Time 

-0.87 

Little  Missouri 

2 

Intercept 

-1.36 

0.0177 

0.79 

River 

Time 

-2.09 

Little  Missouri 

3 

Intercept 

-1.35 

0.0002 

0.92 

River 

Time 

-2.34 

Little  Missouri 

4 

Intercept 

-0.88 

0.0002 

0.95 

River 

Time 

-4.37 

Little  Missouri 

5 

Intercept 

-1.13 

0.0746 

0.86 

River 

Time 

-2.27 

Little  Missouri 

6 

Intercept 

-2.52 

0.0007 

0.92 

River 

Time 

-1.56 

Little  Missouri 

7 

Intercept 

-2.43 

0.0001 

0.97 

River 

Time 

-2.89 

Little  Missouri 

8 

Intercept 

-1.67 

0.0051 

0.94 

River 

Time 

-3.94 

Little  Missouri 

9 

Intercept 

-1.76 

0.0146 

0.97 

River 

Time 

-4.18 

Fourche  La  Fave 

1 

Intercept 

1.05 

0.0014 

0.94 

River 

Time 

-1.52 

Fourche  La  Fave 

2 

Intercept 

0.82 

0.1136 

0.79 

River 

Time 

-0.70 

(Continued) 


*  The  data  set  number  Is  shown  on  the  plots  In  Figure  Cl. 
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Table  2  (Concluded) 


Data  Set 
Number* 

Independent 

Variables 

Regression 

Coefficients 

Equation 

Significance 

_E2_ 

Fourche  La  Fave 

3 

Intercept 

1.21 

0.0001 

0.99 

River 

Time 

-0.71 

Fourche  La  Fave 

4 

Intercept 

0.12 

0.0001 

0.99 

River 

Time 

-1.79 

Fourche  La  Fave 

5 

Intercept 

0.21 

0.0004 

0.97 

River 

Time 

-1.81 

Rough  River 

1 

Intercept 

0.75 

0.0855 

0.56 

(upper) 

Time 

-0.05 

Rough  River 

2 

Intercept 

0.60 

0.2178 

0.89 

(lower) 

Time 

-0.30 

Rough  River 

3 

Intercept 

0.72 

0.1801 

0.40 

(upper) 

Time 

-0.03 

Rough  River 

4 

Intercept 

0.67 

0.0142 

0.99 

(lower) 

Time 

-0.43 

Guadalupe  River 

Intercept 

-2.08 

0.0066 

0.87 

(Canyon) 

Time 

-4.45 

Independent  variables  that  were  the  best  predictors  of  K(.h  .  Many  of  the 
variables  tested  (i.e.,  pH,  [DO],  [OH*],  [€„],  temperature,  etc.)  have  been 
cited  by  other  researchers  as  affecting  manganese  oxidation  (Stiiimn  and  Morgan 
1981,  Hess  1984,  Hslung  1987).  Table  3  lists  the  variables  and  combinations 
of  variables  tested  using  PROC  STEPWISE.  Average  values  for  the  tallwater 
reach  were  used  for  each  variable  In  Table  3.  Hess'  equation  (Equation  4)  was 
the  first  equation  examined  since  It  had  been  used  to  describe  manganese  oxi¬ 
dation  In  tallwaters.  However,  the  Hess  equation  produced  unrealistic  rela¬ 
tionships  (i.e.,  negative  coefficients)  for  some  of  the  Independent  variables, 
thus  Incorrectly  describing  the  processes.  As  a  result,  variations  of  this 
equation  and  many  others  were  examined  to  try  to  explain  the  mechanisms  driv¬ 
ing  the  oxidation  and  loss  of  Mn'*^^.  As  with  the  Hess  equation,  unrealistic 
relationships  were  obtained  for  many  of  the  Independent  variables  examined  in 
the  other  equation  forms.  The  PROC  STEPWISE  analysis  did  not  produce  an  equa¬ 
tion  that  adequately  described  the  rate  coefficients  for  Mn'*'^  removal. 
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Table  3 


Variables  Tested  In  Km,  Rearesslon  Analysis 


Variable 

Variable  Descrlotlon 

[0H-]2 

Hydroxide  Ion  concentration  squared 

[H^]2 

Hydrogen  Ion  concentration  squared 

[Co] 

Initial  Mn**  concentration,  moles/i 

Ratio 

Top  width  (W)  over  hydraulic  depth  (H) ,  ft/ft 

[DO] 

Dissolved  oxygen  concentration,  moles/i 

Po 

Partial  pressure  of  oxygen,  atm 

1/H 

Inverse  of  hydraulic  depth  (H) ,  1/ft 

u«/H 

Ratio  of  shear  velocity  (u*)  to  H,  1/sec 

[Co]/H 

Ratio  of  {C(,]  to  H,  moles/i-ft 

COMBI 

Combination  of  [OH‘]^*Ratlo*Co 

C0HB2 

Combination  of  [OH"]^Ratlo*Co*[DO] 

C0MB3 

Combination  of  [OH‘]**[DO] 

C0MB4 

Combination  of  Ratlo*Co 

COMBS 

Combination  of  [OH*]**Po 

C0MB6 

Combination  of  [OH']^[DO]*Co 

COMB? 

Combination  of  [OH‘]^*Co 

COMBS 

Combination  of  [C^] [DO] 

C0MB7/H 

Ratio  of  tOH']^Co  to  H 

C0MB6/H 

Ratio  of  [0H*)2*[D01*Co  to  H 

COMBS/ [H‘^]  2 

Ratio  of  [Co] [DO]  to  [M"]* 

[Co]/[H^]^ 

Ratio  of  [Co]  to  [M*^]^ 

1/[H^]2 

Inverse  of  [H*]^ 

Note:  All  variables 

are  averages  for  the  tallwater  reach. 

53.  As  discussed  In  Part  II,  many  researchers  suggest  that  the  primary 
removal  mechanism  for  dissolved  manganese  Is  adsorption  which  occurs  rapidly, 
followed  by  slow  oxidation  on  particulates  or  the  substrate.  If  adsorption 
onto  the  substrate  Is  a  primary  removal  mechanism  for  dissolved  manganese.  It 
makes  sense  that  variables  describing  the  rate  of  oxidation  would  have  little 
meaning  In  explaining  the  observed  loss  rate.  Consider  the  mass  balance  of 
dissolved  manganese  In  a  control  volume  of  fixed  volume  V  and  bottom  planar 
area  A  .  The  rate  of  decrease  In  dissolved  manganese  through  adsorption  onto 
the  bottom  can  be  described  by 
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-V,  A  (Mn  -  Mnj,) 


(21) 


V 


dMn 

■HF 


where 

Mn  -  dissolved  (i.e.,  reduced)  manganese  concentration,  mg/i 

V,  -  mass  transfer  velocity  across  diffusive  sublayer  separating 
water  and  substrate,  L/T 

Mn),  -  reduced  manganese  concentration  In  sediment  pore  water,  mg/i 
Assuming  that  the  manganese  In  the  pore  water  Is  rapidly  adsorbed  onto  sedi¬ 
ment  particles,  Mn^,  can  be  assumed  to  be  Insignificant,  and  Equation  21 
becomes 


dMn 


(22) 


where  H  Is  the  hydraulic  depth,  I.e.,  V/A  . 

54.  Equation  22  looks  like  Equation  1  where  the  loss  rate  K^i  Is 
V,/H  .  This  reasoning  helps  to  explain  why  observed  data  Indicate  that  dis¬ 
solved  manganese  Is  rapidly  removed  In  streams  according  to  a  first-order  rate 
law.  From  the  PROC  STEPWISE  regressions,  1/H  was  found  to  be  the  most  sig¬ 
nificant  Independent  variable,  having  a  coefficient  of  determination  (I.e., 

R^)  of  0.62,  or  accounting  for  62  percent  of  the  variation  In  the  dependent 
variable . 

55.  The  next  step  was  to  find  a  way  to  estimate  V,  .  This  led  to  an 
approach  based  on  the  results  published  by  Boudreau  and  Gulnasso  (1982)  that 
relate  the  mass  transfer  of  dissolved  substances  across  the  diffusive  sublayer 
of  the  ocean  bed.  Several  relationships  from  Boudreau  and  Gulnasso  (1982) 
were  used  to  compute  V,  .  These  computed  V,  values  for  each  equation  were 
graphically  compared  with  measured  (I.e.,  based  on  observed  manganese  data) 

V,  values,  which  were  calculated  from  V,  -  H  ,  where  values  were 

taken  from  the  slope  of  the  manganese  versus  time  regressions  (Table  2). 

56.  All  of  the  relationships  from  Boudreau  and  Gulnasso  (1982)  require 
water  viscosity  i/  and  molecular  dlffuslvlty  of  the  solute  D  for  computing 
V,  .  Therefore,  V,  values  were  computed  using  values  of  v  and  D  corre¬ 
sponding  to  the  observed  average  stream  temperature,  T  .  The  value  for  Dj 
(I.e.,  at  temperature  T  )  can  be  computed  (Thibodeaux  1979)  from 


27 


(23) 


D,-D 


T»'25 


where 

Dj  -  molecular  diffusion  coefficient  of  Mn'*’*  at  stream  temperature  T  , 
cmVsec 

D23  -  molecular  diffusion  coefficient  of  Mn'*’*  at  25  “C,  6.88 
X  10'®  cmVsec 

T  -  ambient  stream  temperature,  “C 

1^25  “  kinematic  viscosity  of  water  at  25  "C,  0.93  x  10"®  cmVsec 
i/f  -  kinematic  viscosity  at  stream  temperature  T  ,  cmVsec 
A  regression  for  viscosity  (t/j)  versus  temperature  was  developed  from 
tabulated  values,  resulting  in 

i/j  «  0.000001674  -  3.061667  x  lO*®  T  (24) 

57.  The  following  equation  was  chosen  to  calculate  (i.e.,  predict)  V, 
since  its  absolute  mean  error  (i.e.,  mean  |predicted  •  measured! )  less 
than  the  absolute  mean  error  of  the  other  relationships  in  Boudreau  and 
Guinasso  (1982) : 


D  u,  Sc^'® 
- 7^ - 


(25) 


where 

u*  -  shear  velocity  (VgHs) ,  L/T 
g  -  gravitational  acceleration,  L/T® 
s  -  water  surface  or  streambed  slope,  nondimens ional 

Sc  •  Schmidt  Number  (v/D) ,  nondimens ional 
and  values  for  D  and  1/  are  taken  at  ambient  stream  temperature  T  . 

58.  Graphs  of  predicted  versus  measured  V,  indicated  that  the  data 
could  be  divided  into  two  groups,  one  representing  cobble  bed  streams  and  the 
other  representing  sediment  (i.e.,  fine  grain,  silt  and  clay)  bed  streams.  A 
two-tailed  mean  analysis  was  performed  on  both  groups  to  see  if  there  was 
enough  difference  between  the  mean  V,  for  each  group  to  Justify  having 
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separate  equations  for  V,  for  the  two  streambed  types.  To  perform  the  two- 
tailed  mean  analysis,  the  assumption  had  to  be  verified  that  the  variances  of 
V,  for  the  two  streambed  types  were  equal.  This  assiimption  was  verified 
using  the  F  test  analysis  at  the  significance  level  of  a  -  0.02  .  Upon 
verification  that  the  variances  were  equal,  the  two-tailed  mean  analysis  was 
tested  at  a  significance  level  of  a  -  O.OS  .  It  was  concluded  from  the  anal¬ 
ysis  that  the  mean  V,  values  for  the  two  streambed  types  were  different 
enou^  to  develop  separate  equations  for  V,  .  It  is  not  known  how  a  sandy 
substrate  would  Impact  the  manganese  removal  rate  since  none  of  the  tallwaters 
studied  here  had  significant  quantities  of  sandy  material. 

59.  Linear  regression  analysis  with  PROC  REG  (SAS  Institute,  Inc.  1988) 
was  used  to  relate  the  measured  and  computed  (i.e.,  computed  with  Equation  25) 
Vg  values  for  both  streambed  types,  or 


V.  =  a  V, 


(26) 


where 

V,  -  measured  mass  transfer  rate,  L/T 

O 

a  -  regression  coefficient,  dimensionless 

-  mass  transfer  coefficient  computed  with  Equation  25,  L/T 
Figure  1  shows  the  regression  lines  fitting  the  computed  versus  measured  V, 
for  the  cobble  and  sediment  streambed  types,  respectively.  Regression  coef¬ 
ficients,  significance  value,  and  R^  for  the  Equation  26  model  are  listed  in 
Table  4  for  both  streambed  types.  It  is  reasonable  that  the  cobble  bed  re¬ 
gression  results  in  larger  values  for  V,  since  larger  roughness  elements 
associated  with  cobbles  would  tend  to  increase  turbulence,  thus  increasing  the 
diffusive  sublayer  mass  transfer  velocity. 

60.  The  general  method  to  predict  K^h  can  be  summarized  as  follows. 
Equations  23  and  24  are  first  used  to  obtain  D  and  v  for  the  ambient  water 
temperature.  Next,  V,  is  estimated  with  Equation  25  using  the  previously 
corrected  values  for  D  and  u  and  the  stream  hydraulics  to  obtain  u*  . 
Equation  26  is  used  to  obtain  an  improved  estimate  for  V,  using  the  value 
computed  by  Equation  25.  Finally,  the  result  obtained  with  Equation  26  is 
divided  by  H  to  obtain  K^i  •  A  plot  of  predicted  versus  observed  K^, 
using  this  methodology  is  shown  in  Figure  2.  The  line  in  Figure  2  is  the 
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Figure  1.  Measured  versus  computed  V,  for  Mn  data  divided 
into  cobble  and  sediment  substrates 


Table  4 

Regression  Results  for  Equation  26 


Streambed 

Independent 

Regression 

Equation 

HHl 

Cobble 

V. 

1.654 

0.0001 

0.89 

Sediment 

V. 

0.579 

0.0057 

0.88 
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Figure  2.  Predicted  versus  observed  dissolved 
manganese  removal  rate, 

least  squares  best  fit.  The  slope  of  this  line  is  .68  and  has  an  R*  value  of 
0.74.  The  regression  slope  is  not  significantly  different  (o  -  0.05)  from  a 
perfect  correlation  slope  of  1.0. 

Reduced  Iron  Removal 

61.  It  was  difficult  to  interpret  the  reduced  iron  data  taken  at  the 
four  study  sites.  Reduced  iron,  l.e.,  ferrous  iron,  was  defined  as  dissolved 
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Iron,  or  the  iron  that  passed  through  a  0.45-/i  filter.  The  dissolved  Iron 
data  from  the  Greeson  and  Canyon  tallwaters  were  very  low  (less  than  0.1  mg/i) 
and  showed  little  change  downstream.  Higher  concentrations  of  total  and  dis¬ 
solved  Iron  were  measured  at  Nimrod  and  Rough  River  tallwaters.  The  dissolved 
Iron  showed  little  change  In  the  downstream  direction  of  the  Nimrod  tallwater, 
whereas  at  Rough  River  there  was  a  downstream  decline.  Nix  et  al.  (1991) 
speculated  that  fine  floe  of  ferric  Iron  could  have  passed  the  0.45-/i  filter, 
yielding  artificially  high  dissolved  Iron  concentrations  In  the  downstream 
reaches  of  the  tallwaters.  Gordon,  Bonner,  and  Milligan  (1984)  recommended 
that  a  0. 10- fi  filter  be  used 'to  filter  particulate  and  colloidal  Iron. 

Gordon,  Bonner,  and  Milligan  (1984)  did  not  detect  sulfides  below  Normandy 
Dam,  but  found  total  Iron  to  be  considerably  higher  than  dissolved.  Indicating 
the  possible  presence  of  Iron  sulfide.  At  Nimrod  and  Rough  River  tallwaters, 
the  total  Iron  was  about  two  to  three  times  greater  than  the  dissolved  Iron, 
and  there  was  an  Increase  In  sulfate  In  the  downstream  direction.  Therefore, 
It  was  concluded  that  Iron  sulfide  was  probably  released  Into  both  tallwaters, 
thus  complicating  the  Iron  dynamics  (Nix  et  al.  1991). 

62.  The  Information  yielded  from  the  Iron  data  was  far  less  than  desir¬ 
able  for  evaluation  of  Iron  oxidation  (l.e.,  removal)  rates.  However,  three 
of  the  data  sets  that  clearly  showed  a  downstream  decline  In  dissolved  Iron 
concentration  were  analyzed  for  order  of  magnitude  rates  of  oxidation  and  were 
compared  with  predicted  values.  These  three  data  sets  were  Greeson,  September 
1983  (Nix  1986);  Greeson,  October  1987  (Nix  et  al.  1991);  and  Rough  River, 
August  1988  (Nix  et  al.  1991).  The  observed  changes  In  dissolved  Iron  concen¬ 
tration  and  travel  time  were  used  to  estimate  an  "observed"  oxidation  rate  of 
dissolved  Iron  using  the  first-order  rate  law  (Equation  1) .  The  "predicted" 
oxidation  rate  at  20  "C  was  computed  with  Equation  6  using  the  average  DO  and 
pH  observed  throughout  the  tallwater  reach  from  which  the  "observed"  oxidation 
rate  was  estimated.  The  predicted  rates  were  then  corrected  for  the  tempera¬ 
ture  observed  In  the  reach  using  Equation  7. 

63.  The  following  results  were  obtained  for  the  three  data  sets: 


Observed  Kp,  Predicted  Kp. 

Date  dav~^  dav'^ 


Greeson 

Sep  1983 

.  — 

0.54 

.  ■■  J 

0.12 

Greeson 

Oct  1987 

6.65 

1.16 

Rough  River 

Aug  1988 

1.25 

150. 
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The  fact  that  the  predicted  rates  of  iron  oxidation  do  not  agree  with  those 
observed  is  not  surprising  considering  the  complex  Interactions  surrounding 
iron  and  iron  sulfide  dynamics,  as  discussed  in  Part  II.  These  results  at 
least  show  that  Equation  6  can  provide  the  proper  order  of  magnitude  estimate 
for  Kp,,  as  it  did  for  the  Greeson  site.  Gordon,  Bonner,  and  Milligan  (1984) 
measured  dissolved  iron  below  Normandy  Dam  that  exhibited  removal  rates 
between  about  1.8  and  18.0  day'^. 

64.  The  observed  rates  were  higher  than  the  predicted  rates  at  the 
Greeson  site,  which  might  be  explained  by  the  presence  of  other  trace  metals 
(Stumm  and  Morgan  1981,  Nix  1986).  Higher  rates  were  observed  and  predicted 
for  the  October  1987  Greeson  study  than  for  the  earlier  Greeson  study;  this 
may  be  partially  due  to  the  higher  pH  observed  in  the  later  study.  A  high 
oxidation  rate  was  predicted  for  the  Rough  River  study  because  of  much  higher 
pH  (i.e.,  pH  -  7.5  versus  about  5.6  and  6.2  for  the  Greeson  studies).  The 
relatively  low  oxidation  rate  observed  at  Rough  River  may  not  be  a  true 
indication  of  the  actual  oxidation  rate  due  to  the  suspected  presence  of  iron 
sulfide  and  the  associated  sampling/analysis  problems  discussed  above.  If 
colloidal  iron  had  passed  the  filter  during  the  dissolved  iron  measurements, 
the  dissolved  iron  measurements  would  be  artificially  high,  thus  lowering 
estimates  for  the  oxidation  rates.  Additionally,  the  sulfate  concentrations 
at  Rough  River  were  about  three  times  higher  than  those  at  Greeson.  As  men¬ 
tioned  in  Part  II,  high  sulfate  concentrations  can  have  a  retarding  influence 
on  the  iron  oxidation  rate.  There  was  simply  not  enough  information  to  allow 
definitive  analysis  of  iron  sulfide  dynamics. 

Sulfide  Removal 


65.  Very  little  information  could  be  extracted  from  the  field  studies 
for  sulfide  removal.  Sulfide  was  detected  in  small  quantities  at  Greeson 
tailwater  in  1983  (Nix  1986) ,  Nimrod  tailwater  in  1989  (Appendix  A) ,  and  Can¬ 
yon  tailwater  in  1988  (Nix  et  al.  1991).  Analyses  of  these  three  data  sets 
are  discussed  below. 

66.  The  1983  Greeson  data  indicate  that  sulfide  decreased  from  about 
0.3  mg/i  to  0.04  mg/i  over  the  first  10  km  of  the  tailwater.  The  travel  time 
for  this  reach  and  flow  rate  reported  by  Nix  (1986)  are  now  considered  to  be 
in  error;  the  correct  travel  time  is  approximately  3.8  hr.  This  decrease  in 
sulfide  over  3.8  hr  translates  into  a  loss  rate  of  about  12.7  day"^.  The  pH 
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during  this  study  was  approximately  5.5,  which  means  that  most  (i.e.,  about 
95  percent)  of  the  sulfide  existed  as  H2S  and  was  lost  by  volatilization.  A 
hydrogen  sulfide  volatilization  rate  of  12.7  day~^  is  of  the  same  order  of 
magnitude  as  the  reaeration  rate  of  about  7.0  day~^  that  was  measured  (Dortch 
and  Hamlin  1988) .  Recall  from  Part  II  that  the  reaeration  rate  can  be  used  to 
estimate  the  hydrogen  sulfide  volatilization  rate. 

67.  Sulfide  was  found  to  decrease  at  Nimrod  in  1989  from  about 

0.017  mg/i  just  below  the  dam  to  about  0.012  mg/i  at  the  next  downstream  sam¬ 
pling  station,  which  is  about  1.8  hr  travel  time  at  the  high-flow  discharge. 
About  75  percent  of  the  measured  sulfide  was  H2S  for  the  observed  pH  of  6.5. 
Therefore,  HS"  decreased  from  0.0043  to  about  0.0030  mg/i.  This  decrease 
results  in  an  overall  sulfide  oxidation  rate  of  about  4.6  day"^.  Correcting 
for  DO  and  temperature  (see  Equations  8  and  9,  respectively)  provides  an  esti¬ 
mated  k,  in  Equation  8  of  about  4.9  day"^  at  20  ”C.  This  rate  is  on  the  low 
side  of  the  range  of  reported  rates  discussed  in  Part  II. 

68.  Nix  et  al.  (1991)  reported  that  the  sulfide  at  Canyon  tailwater 
decreased  from  about  0.04  mg/i  just  below  the  dam  to  background  levels  at  the 
next  downstream  station.  The  travel  time  to  the  first  station  was  70  min. 
Based  on  the  data,  it  was  assumed  that  background  levels  were  about  0.01  mg/i. 
For  the  observed  pH  of  7.6,  80  percent  of  the  sulfide  is  HS".  This  translates 
into  an  overall  oxidation  rate  for  HS"  of  about  24.0  day"^.  Correcting  for 
temperature  and  DO  as  above,  k,  at  20  "C  is  36.7  day"^.  This  value  falls 
about  in  the  middle  of  the  reported  ranges. 

Nitrification 

69.  The  ammonia  data  reported  by  Nix  et  al.  (1991)  were  analyzed  in  a 
manner  similar  to  the  manganese  data.  Ten  data  sets  (Table  5)  that  clearly 
showed  a  downstream  decrease  in  ammonia  were  used  in  the  analysis.  Linear 
regressions  were  developed  for  the  natural  log  of  the  ammonia  concentrations 
versus  travel  time  so  that  the  overall  ammonium  loss  (i.e.,  nitrification  rate 
K,  ,  day"^)  could  be  determined  from  the  slope  of  the  regression.  Plots  of 
the  data  and  the  regression  lines  are  presented  in  Appendix  C.  The  K,  value 
obtained  from  the  regressions  (i.e.,  measured  K,) ,  the  correlation  coeffi¬ 
cient,  and  significance  of  the  regression  are  shown  in  Table  5.  Data  were 
excluded  from  further  analysis  if  the  significance  probability  of  the  equation 
was  greater  than  0.05.  Separation  of  the  ammonia  data  into  upper  and  lower 
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Table  5 

Regression  Results  for  HHT  Concentration  Versus  Time  from  Measured  Data 


Data  Set 
Number* 

Independent 

Variables 

Regression 

Coefficients 

Equation 

Significance 

_Ri_ 

Little  Missouri 

Intercept 

-2.69 

0.2194 

0.61 

River 

Time 

-4.33 

Fourche  La  Fave 

1 

Intercept 

-0.46 

0.0007 

0.96 

River 

Time 

-1.35 

Fourche  La  Fave 

2 

Intercept 

-0.86 

0.1587 

0.71 

River 

Time 

-0.77 

Fourche  La  Fave 

3 

Intercept 

-0.57 

0.0097 

0.92 

River 

Time 

-0.66 

Fourche  La  Fave 

4 

Intercept 

-1.32 

0.0001 

0.99 

River 

Time 

-1.79 

Fourche  La  Fave 

5 

Intercept 

-1.25 

0.0007 

0.96 

River 

Time 

-1.32 

Rough  River 

1 

Intercept 

0.13 

0.0029 

0.74 

Time 

-0.84 

Rough  River 

2 

Intercept 

0.08 

0.0116 

0.62 

Time 

-0.13 

Rough  River 

3 

Intercept 

0.09 

0.0012 

0.80 

Time 

-0.15 

Guadalupe  River 

Intercept 

-2.18 

0.0003 

0.90 

Time 

-2.39 

*  The  data  set  number  is  shown  on  the  plots  in  Figure  C2. 


reaches  for  Rough  River  resulted  in  insignificant  regression;  thus,  the  ammo¬ 
nia  data  were  not  separated. 

70.  As  discussed  in  Fart  II,  nitrification  is  considered  to  be  affected 
by  stream  hydraulics,  i.e.,  contact  with  the  substrate.  Therefore,  the  meth¬ 
ods  outlined  earlier  for  computing  the  diffusive  sublayer  mass  transfer  veloc¬ 
ity  V,  were  used  here,  with  the  exception  that  D20  for  ammonia  is  1.76 
X  10'’  cm^/sec,  and  a  reference  teiaperature  of  20  *C  rather  than  25  *C  is  used 
in  Equation  23.  The  quantity  V,/H  was  correlated  with  the  measured  K,  , 
yielding  a  regression  coefficient  of  0.136  with  an  R^  value  of  0.83  and  a 
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significance  of  0.0006.  Thus,  aononia  is  nodeled  with  a  first-order  removal, 
with  the  removal  (nitrification)  rate  prescribed  by 

K,  -  0.136  ^  (27) 

71.  A  plot  of  predicted  (with  Equations  23-25  and  27)  versus  observed 
is  shown  in  Figure  3.  The  slope  of  the  least  squares  best  fit  regression 
line  shown  in  Figure  3  is  0.65  and  has  an  value  of  0.77.  This  slope  is  not 
significantly  different  (a  -  0.05)  from  the  perfect  correlation  slope  of  1.0. 

A  temperature  correction  was  applied  to  K,  with  an  equation  of  the  same  form 
as  Equation  7  and  a  value  of  1.08  for  the  base  of  the  power  term. 
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Figure  3.  Predicted  versus  observed  ammonla-N 
removal  rate,  K, 


PART  IV:  NUMERICAL  MODEL 


Description 


General 

72.  The  focus  of  the  Tailwater  Quality  Numerical  Model  (TWQM)  is  to 
predict  the  downstream  transformation  of  problem  constituents,  such  as  DO  and 
several  of  the  reduced  substances  (i.e.,  iron,  manganese,  sulfide,  and 
ammonia) .  It  is  intended  that  the  model  require  as  little  input  data  as  pos¬ 
sible,  allowing  modest  resources  for  application.  The  obvious  way  to  reduce 
input  requirements  is  to  restrict  the  model  dimensionality;  the  model  is  1-D 
in  the  streamwise  direction.  Another  way  to  reduce  input  requirements  is  to 
restrict  the  model  to  steady- state  applications  so  that  time -varying  input 
files  are  not  required.  Therefore,  TWQM  does  not  allow  time -varying  flow  or 
water  quality  Inputs,  and  a  steady-state  solution  is  performed. 

73.  Steady-state  solutions  for  1-D  problems  can  often  be  performed  ana 
lytically,  rather  than  numerically  (Dortch  and  Martin  1989).  However,  a 
numerical  modeling  approach  was  selected  for  the  TWQM  to  allow  flexibility  in 
interactions  of  state  variables.  Additionally,  with  the  microcomputer 
resources  available  today,  numerical  solutions  are  easily  and  rapidly 
attainable . 

74.  TWQM  is  based  on  the  US  Environmental  Protection  Agency  (EPA)  1-D 
stream  water  quality  model,  QUAL2E  (Brown  and  Barnwell  1987).  The  major  modi 
fication  to  QUAL2E  was  to  include  subroutines  to  model  the  reduced  substances 
in  tallwaters  and  a  user  interface  to  provide  easier  model  application. 

QUAL2E  was  selected  for  several  reasons:  (a)  the  model  is  widely  used  and 
accepted,  (b)  the  model  provides  for  steady-state  solutions,  (c)  the  model 
contains  other  water  quality  processes  that  could  be  used  if  desired,  and 

(d)  the  most  recent  version  of  QUAL2E  had  been  developed  for  microcomputers 
with  simpler  input  requirements. 

75.  The  TWQM  combines  sfc\ 'ral  components  within  a  user-friendly  inter¬ 
face  that  utilizes  menus  for  selecting  alternatives.  In  addition  to  the 
QUAL2E-based  tailwater  component,  there  is  a  reservoir  release  component. 
Reservoir  release  water  quality  is  required  for  the  upstream  boundary  condi¬ 
tion  of  the  stream  component.  The  user  may  specify  this  information  from 
observations  or  can  predict  the  release  concentrations  based  upon  observed  in 
pool  concentrations  using  the  SELECT  model  (Davis  et  al.  1987).  Given  the 
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in-pool  temperature  stratification,  the  outlet  features,  release  flow  rate, 
and  vertical  distribution  of  water  quality  constituents,  SELECT  computes  the 
release  concentrations.  Also  Included  in  SELECT  is  a  structural  reaeration 
component  that  predicts  uptake  of  dissolved  oxygen  as  flow  passes  through  the 
release  structure.  Release  concentrations  are  provided  by  SELECT  to  the  tail- 
water  component  for  predicting  downstream  concentrations.  Except  for  DO, 
neither  SELECT  nor  the  tallwater  model  alters  concentrations  of  water  quality 
variables  as  flow  passes  through  the  outlet  structure;  in  other  words,  vari¬ 
ables  are  treated  as  conservative  as  they  pass  through  the  structure  because 
of  the  short  residence  time  in  the  structure.  SELECT  does  allow  the  option  of 
structural  reaeration. 

76.  The  combination  of  SELECT  with  the  tallwater  model  allows  the  user 
to  evaluate  the  impact  on  downstream  water  quality  of  various  release  schemes, 
such  as  hydropower  retrofit  or  moving  the  vertical  location  of  the  intakes  of 
the  outlet  structure. 

77.  TWQM  has  the  capability  of  simulating  up  to  17  water  quality  con¬ 
stituents.  Constituents,  which  can  be  modeled  in  any  combination  by  the  user, 
are: 

a.  DO. 

fe.  CBOD. 

£ .  Temperature . 

d.  Algae  as  chlorophyll  a. 

e.  Total  organic  nitrogen. 

£.  Ammonia  nitrogen. 

g.  Nitrite  and  nitrate  nitrogen. 

h.  Total  organic  phosphorus. 

1.  Dissolved  Inorganic  (orthophosphate)  phosphorus. 

i.  Dissolved  (reduced)  iron. 

ic.  Dissolved  (reduced)  manganese. 

i.  Total  dissolved  sulfide  (HS'  and  H2S) . 

S.  Iron  sulfide. 

E.  Arbitrary  nonconservative  constituent. 

£.  Two  conservative  constituents. 

Most  of  the  above  constituents  were  originally  in  QUAL2E;  the  reduced  sub¬ 
stances  iron,  manganese,  sulfide,  and  iron  sulfide  were  added.  The  model's 
compartmental  diagram  (Figure  4)  Illustrates  the  various  state  variables  and 
how  they  interact.  The  model  solves  the  steady- state  mass  balance  equation 
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SEDIMENTS 

S  - 

SETTLING 

UO  - 

OXYGEN  UPTN<E  FROM  DECAY 

R  - 

RESPNATION 

NO  - 

OXYGEN  UPTAKE  FROM  NITRIFICATION 

P  - 

PHOTOSYNTHESIS 

UG  - 

ALGAE  UPTAKE  FOR  GROWTH 

H  - 

HYDROLYSIS 

RO  - 

ALGAE  RELEASE  FROM  DEATH,  EXCRETION 

N  - 

NITRIFICATION 

BR  - 

BENTHIC  RELEASE 

OR  - 

OXYGEN  REAERATION 

OR  - 

OXYEN  DEAERATION 

AD  - 

ADSORPTION 

SOO  - 

SED.  OXYGEN  DEMAND 

CO  - 

CHEMICAL  OXIDATION 

V  - 

VOLITILIZATION 

ON  - 

SEDIMENT  DENITRIFICATION 

LO  - 

LOST  THROUGH  OXIDATION 

Figure  4.  Tailwater  quality  model  compartmental  diagram 
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(l.e.,  mass  transport  equation  or  advection-diffuslon  equation  with  mass 
sources/slnks)  for  each  state  variable.  The  mass  balance  equation  Is  referred 
to  as  the  energy  equation  when  temperature  Is  the  state  variable. 

78.  The  tallwater  Is  physically  described  by  subdividing  the  stream 
system  Into  reaches  (the  basic  division  of  the  model).  Reaches  represent 
portions  of  the  river  having  similar  channel  geometry,  hydraulic  characteris¬ 
tics,  and  chemlcal/blologlcal  coefficients.  Reaches  are  further  divided  Into 
equally  spaced  units  called  computational  elements,  or  nodes.  An  example 
schematic  of  a  modeled  tallwater  system  Is  Illustrated  In  Figure  5.  Each 
computational  element  has  Inputs,  outputs,  and  reaction  terms.  The  energy  and 
mass  balance  equations  are  solved  simultaneously  (Implicitly)  for  all  computa¬ 
tional  elements.  This  solution  Is  repeated  In  an  Iterative  fashion,  using  the 
previous  solution  results,  until  convergence  for  the  steady- state  solution  Is 
reached  (l.e.,  the  concentrations  quit  changing  within  a  small  tolerance). 

MOST  UPSTREAM 
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(RESERVOIR  RELEASES) 
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Figure  3.  Example  schematic  of  a  model  stream  system 
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79.  The  model  assumes  steady,  nonuniform  flow.  Hydraulic  conditions 
(velocity  and  depth)  used  within  the  energy  and  mass  balance  equations  are 
entered  for  each  reach  and  assumed  uniform  throughout  the  reach.  Steady  flow 
implies  that  the  flow  rate,  velocity,  width,  and  depth  at  any  point  in  the 
stream  are  constant  with  time.  Nonuniform  flow  is  allowed,  which  means  flow 
rate,  velocity,  width,  and  depth  can  change  in  the  longitudinal  direction  from 
reach  to  reach. 

80.  Model  processes  of  state  variables  that  were  added  or  kinetics  that 
were  modified  are  described  in  the  sections  that  follow.  As  with  any  numeri¬ 
cal  solution,  errors  can  be  associated  with  the  solution  scheme.  A  descrip¬ 
tion  of  the  numerical  error  inherent  in  this  model  and  the  remedy  are 
discussed  in  the  next  section  on  numerical  diffusion. 

Numerical  diffusion 

81.  The  version  of  QUAL2E  from  which  TWQM  was  developed  required  the 
user  to  specify  the  ntunber  of  computational  elements  that  a  reach  would  be 
divided  into.  QUAL2E  allows  a  maximum  of  20  elements  per  reach,  which  is  not 
always  adequate.  Computational  element  length  can  affect  the  results  of  a 
simulation.  Computational  element  lengths  that  are  too  long  can  result  in 
numerical  diffusion. 

82.  Numerical  diffusion  occurs  when  the  mathematical  solution  scheme 
used  in  a  model  incorrectly  contributes  to  the  diffusion  of  modeled  sub¬ 
stances.  For  a  model  to  be  representative,  predicted  movement  of  substances 
should  be  a  function  of  the  advective  velocity  of  the  water  and  the  dispersion 
characteristics  of  the  stream  (Fischer  et  al.  1979).  For  example,  if  the 
movement  of  a  tracer  in  a  river  is  to  be  modeled,  the  computational  elements 
should  have  a  length  that  causes  the  advective  velocity  (i.e.,  displacement 
velocity  or  displacement  distance  divided  by  elapsed  time)  of  the  tracer  mean 
concentration  to  be  the  same  as  the  mean  flow  velocity  of  the  stream.  Also, 
the  variance  of  this  tracer  cloud  should  increase  at  a  rate  of  about  two  times 
the  dispersion  coefficient  for  a  normally  distributed  cloud.  When  these  con¬ 
ditions  are  not  met,  the  distribution  of  the  tracer  is  not  only  a  function  of 
the  dispersion  coefficient  but  also  of  the  numerical  solution  used. 

83.  To  decrease  the  effects  of  numerical  diffusion  in  TWQM,  the 
computational  element  length  was  determined  by  equating  the  estimated  longitu¬ 
dinal  dispersion  of  the  stream  with  numerical  dispersion.  The  stream  disper¬ 
sion  can  be  estimated  from  Fischer  et  al.  (1979). 
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where 


u  -  average  stream  velocity,  m/sec 
W  -  top  width  of  stream,  m 

Numerical  dispersion  for  steady-state  conditions  can  be  computed  (Roache  1972) 
from 


^  - 


u  Ax 


(29) 


where  Ax  Is  the  computational  element  length  (meters).  Variables  in  Equa¬ 
tions  28  and  29  are  shown  with  the  appropriate  metric  units.  TWQM  has  an 
option  that  allows  expression  of  distances,  flows,  velocities,  and  temperature 
In  non-metrlc  units.  Solving  Equations  28  and  29  yields  a  value  for  Ax  of 


.  2(0.011)u2w2 
- IMT. - 


(30) 


Substituting  for  u*  and  expressing  W  in  terms  of  flow  rate  (Q) ,  depth,  and 
average  velocity  yields  Equation  31, 


Ax  = 


2(0.011)Q2 

uH3-5gO-5s<»-5 


(31) 


84.  The  total  reach  length  Is  divided  by  the  value  of  Ax  determined 
from  Equation  31  to  determine  the  total  number  of  elements  per  reach.  This 
number  Is  rounded  to  the  nearest  integer  so  as  not  to  have  a  fraction  of  an 
element  within  a  reach.  The  Integer  number  of  elements  Is  divided  Into  the 
reach  length  to  define  the  computational  element  length  for  that  reach.  The 
maximum  allowable  number  of  computational  elements  per  reach  was  Increased  In 
TWQM  to  100,  which  limits  the  minimum  length  of  any  computational  element  to 
1/100  of  the  reach  length.  This  limit  was  arbitrarily  set  to  restrict  storage 
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requirements  of  TWQM,  but  the  user  can  easily  increase  this  limit  through  a 
parameter  statement  in  the  source  code. 

85.  To  test  the  effectiveness  of  these  changes,  simulations  were  per¬ 
formed  in  which  a  nonconservative  tracer  was  continuously  injected  at  the  dam 
and  its  steady-state  concentration  was  predicted  at  distances  downstream.  In 
the  case  where  the  computational  element  length  was  determined  (Figure  6a) , 
the  predicted  concentrations  were  very  close  to  those  obtained  using  the 
analytical  solution  (Fischer  et  al.  1979). 

[|(^TT.i)jj  (32) 


where 

Co  -  release  concentration  at  x  -  0,  mg/i 

k  -  reaction  rate  constant,  day"^ 

4D  k 

u 

For  the  siiaulation  in  which  the  computational,  element  length  was  set  at  an 
arbitrary  length  of  2  miles  (Figure  6b) ,  the  predicted  concentration  is  not  as 
close  to  the  analytical  solution.  Both  simulations  were  performed  using  the 
same  Initial  concentrations,  velocities,  and  channel  geometry.  The  only  dif¬ 
ference  between  the  two  simulations  was  the  computational  element  length. 

86.  QUAL2E,  and  therefore  TWQH,  require  that  the  length  of  computa¬ 
tional  elements  be  the  same  in  all  reaches.  In  simulations  performed  using 
multiple  reaches,  the  minimum  computational  element  length  determined  in  any 
reach  is  used  for  all  reaches.  In  this  case  the  reach  length  of  the  other 
reaches  might  be  increased  or  decreased  slightly  to  prevent  a  reach  from  con¬ 
taining  a  fraction  of  an  element. 

87.  QUAL2E  requires  the  user  to  input  the  element  type  (headwater, 
standard,  point  load,  end)  for  each  computational  element.  In  TWQM  this  pro¬ 
cess  has  been  automated  due  to  the  potentially  large  number  of  elements  used 
per  reach.  The  only  user  input  requirements  of  this  type  are  the  locations  in 
river  miles  of  point  loads  (including  waste  discharges  and  tributaries)  or 
withdrawals  within  a  reach. 
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Figure  6.  Results  of  stream  transport  test 


Dissolved  manyaneae 

88.  Dissolved  (l.e.,  reduced)  nanganese  removal  Is  modeled  as  a  first- 
order  reaction  using  the  equation 

^  -  -K^  (33) 

where 

Mn  -  Mn"''^  concentration,  mg/i 
-  Mn«  loss  rate ,  day"^ 

The  manganese  loss  mechanism  Is  assumed  to  be  predominantly  through  adsorption 
onto  the  bed  as  discussed  In  Part  HI.  The  user  has  a  choice  of  having  the 
model  calculate  the  loss  rate  K^q  as  described  In  the  Reduced  Manganese 
Removal  section  of  Part  III  (l.e.,  Equations  23-26),  or  specifying  Ki^  .  A 
default  value  Is  used  If  the  user  does  not  elect  to  have  the  model  pre¬ 

dict  K(^  or  does  not  specify  Kmq  .  If  the  user  elects  to  model  the  removal 
rate,  the  streambed  type  (l.e.,  cobble  or  sediment)  of  the  tallwater  must  be 
specified.  If  a  K^,  value  Is  specified  or  the  default  value  Is  used,  the 
value  entered  Is  assumed  to  be  at  20  *C,  and  the  following  equation  Is  used  to 
correct  for  temperature  dependence: 

Kr  »  (34) 


where 

Kj  -  Mn'*’*  removal  rate  at  stream  temperature  T 
K20  -  Mn'*'*  removal  rate  at  20  "C 
8  -  temperature  correction  constant 

Hess,  Byung,  and  Roberts  (1989)  recommend  a  value  of  1.054  for  $  for  tem¬ 
peratures  below  25  "C.  Evaluation  of  data  reported  by  Hslung  (1987)  suggested 
a  value  of  1.105  for  8  for  temperatures  above  25  *C. 

Dissolved  Iron 

89.  Dissolved  (l.e.,  reduced)  Iron  removal/oxldatlon  Is  also  modeled  as 
a  first-order  reaction  with  respect  to  dissolved  Iron  concentration.  The  loss 
rate  Kp,  can  be  calculated  by  the  model  (l.e.,  using  Eqviatlons  6  and  7)  or 
specified  by  the  user.  Again,  If  the  user  does  not  specify  a  value  or  does 
not  elect  to  let  the  model  compute  the  loss  rate,  a  default  value  will  be 
used.  When  the  default  value  Is  used  or  a  value  Is  specified,  the  temperature 
dependence  Is  taken  Into  account  using  Equation  7. 
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Dissolved  sulfide 

90.  Total  dissolved  sulfide  is  a  modeled  state  variable  composed  of  HS~ 
and  HjS.  The  fractions  of  H2S  and  HS~  present  are  estimated  based  on  the  pH- 
dependent  relationship  in  Standard  Methods  (American  Public  Health  Association 
et  al,  1981).  For  the  pH  range  of  5  to  7,  the  fraction  of  total  dissolved 
sulfide  as  H2S  (FEIACT)  can  be  calculated  as 


FRACT 


-8.3813X3  ♦  135.898X2  -  742. 31X  ♦  1,461.46 

- run - 


(35) 


where  X  is  the  pH  of  the  stream.  For  the  pH  range  of  7  to  9,  the  fraction 
of  total  dissolved  sulfide  as  H2S  is  calculated  as 

TOArT  »  -8.25406X3  +  213.821X2  _  1.852. 79X  +  5,374.11  ^6) 

run 

The  fraction  of  dissolved  sulfide  as  HS~  is  then  calculated  as  1-FRACT.  If  pH 
is  less  than  5  or  greater  than  9,  FRACT  is  set  to  1  or  0,  respectively.  An 
insignificant  amount  of  dissolved  sulfide  is  made  up  of  S~2  (American  Public 
Health  Association  et  al.  1981). 

91.  Removal  of  HS'  through  oxidation  is  modeled  using  Equation  8.  The 
oxidation  rate  k,  can  be  specified  by  the  user;  otherwise,  a  default  value 
of  25.0  day'^  is  used.  The  temperature  correction  for  k,  is  made  with 
Equation  9. 

92.  Loss  of  H2S  gas  through  volatilization  is  modeled  using  Equation  11 
where  the  volatilization  rate  K,  is  set  equal  to  the  reaeration  coefficient 
K2  ,  as  discussed  in  Part  II.  The  reaeration  coefficient  is  modeled  according 
to  the  formula  selected  by  the  user,  as  discussed  in  Appendix  B. 

Iron  sulfide 

93.  Although  the  oxidation  and  loss  of  iron  sulfide  may  be  complex,  a 
simple  approach  for  modeling  is  proposed  here,  considering  the  lack  of  better 
Information  and  supporting  data  at  this  time.  The  loss  of  iron  sulfide  from 
the  tallwater  is  accounted  for  by  two  mechanisms,  first-order  oxidation  (i.e.. 
Equation  1)  in  the  water  column  and  settling  to  the  bottom.  Thus,  the  kinetic 
equation  for  iron  sulfide  is 
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d(FeS) 

at 


(37) 


[Xfs  ♦  5  j  (FeS) 


where 

(FeS)  -  iron  sulfide  concentration,  mg/i 
-  first-order  oxidation  rate,  day'^ 

U.  -  settling  rate  of  FeS,  m/day 

The  settling  and  oxidation  rates  are  either  specified  by  the  user  or  set  to 
the  default  values  of  1.0  m/day  and  38.0  day"^,  respectively. 

Hitrification  and  denitrification 

94.  QUAL.2E  contains  algorithms  for  the  nitrogen  cycle,  including  miner¬ 
alization  of  organic  nitrogen  to  ammonium  nitrogen,  nitrification  from 
ammonium  to  nitrite  and  from  nitrite  to  nitrate  nitrogen,  and  denitrification. 
These  algorithms  were  retained  for  use  in  the  TVQM  with  some  minor  modifica¬ 
tions.  The  method  for  computing  the  nitrification  rate  discussed  earlier  in 
Part  111  (i.e..  Equation  27)  was  included  as  a  user  option.  Additionally, 
nitrite  nitrogen  was  removed  as  a  state  variable  since  it  oxidizes  rapidly  and 
its  presence  only  further  complicates  the  model.  Therefore,  ammonium  nitrogen 
is  oxidized  directly  to  nitrate  nitrogen.  A  first-order  removal  mechanism  for 
nitrate  nitrogen  was  added  to  represent  denitrification  by  sediments.  This 
denitrification  mechanism  is  Independent  of  DO  concentration,  and  the  user  can 
specify  the  first-order  rate  or  use  the  default  value. 

Dissolved  oxygen 

95.  The  equation  for  dissolved  oxygen  in  QUAL2E  had  to  be  modified  to 
allow  proper  representation  of  the  various  DO  sinks  caused  by  the  oxidation  of 
reduced  substances  in  the  tailwater.  The  sinks  of  DO  are: 

fi.  Algal  respiration. 

Decomposition  of  organic  carbon  (i.e.,  CBOD  decay). 

£.  Nitrification  of  ammonium  nitrogen  to  nitrate  nitrogen. 
j|.  Oxidation  of  reduced  manganese  sorbed  onto  the  bed. 

£.  Oxidation  of  reduced  (i.e.,  dissolved)  iron. 

£.  Oxidation  of  iron  sulfide,  both  in  the  water  column  and  in  the 
bed. 

g.  Oxidation  of  HS'. 
b-  Sediment  oxygen  demand  (SOD). 
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96.  CBOD  Is  used  to  represent  the  oxygen  uptake  due  to  decomposition  of 
organic  carbon.  The  kinetics  and  oxidation  of  CBOD  decay  and  algal  respira¬ 
tion  In  QUAL2E  were  not  modified.  The  oxygen  uptake  associated  with  nitrifi¬ 
cation  was  altered  to  reflect  that  the  Intermediate  step  of  nitrification  to 
nitrite  was  eliminated.  Oxygen  uptake  due  to  oxidation  of  dissolved  Iron, 

Iron  sulfide  in  the  water  column,  and  HS"  were  Included  by  subtracting  the 
three  corresponding  stoichiometric  sink  terms  from  the  right  side  of  the  oxy¬ 
gen  mass  balance  equation.  QUAL2E  contained  a  loss  mechanism  of  DO  through 
SOD.  The  SOD  for  TWQM  Is  considered  to  be  caused  by  the  decomposition  of 
organic  material  In  the  bottom  sediments. 

97.  The  oxidation  of  Inorganic  matter  on  the  bottom  (l.e.,  sorbed, 
reduced  manganese  and  Iron  sulfide  that  settles  to  the  bottom)  must  also  be 
Included  In  the  DO  balance.  Matter  lost  from  the  water  column  accumulates  In 
the  bed  for  subsequent  oxidation.  The  concentration  of  manganese  and  Iron 
sulfide  In  the  bed  can  be  derived  from  the  steady- state  mass  balance  of  these 
two  substances  In  a  bed  control  volume 


wA,C 


(38) 


where 

C),  -  concentration  of  reduced  manganese  or  Iron  sulfide  In  bed,  mg/i 

w  -  loading  rate  (l.e.,  mass  transfer  rate  or  settling  rate)  to 
bed,  m/day 

A,  -  surface  area  of  bed  control  volume,  m^ 

C  -  concentration  of  reduced  manganese  or  Iron  sulfide  In  water 
column,  mg/f 

Kb  **  oxidation  rate  of  reduced  manganese  or  Iron  sulfide  In  bed, 

1/day 

Vb  -  volume  of  bed  control  volume,  m^ 

The  oxidation  of  these  substances  creates  an  additional  sediment  oxygen  demand 
that  can  be  calculated  by 


ASOD 


~K- 


awC 


(39) 
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where 


ASOD  -  additional  SOD  resulting  from  manganese  or  Iron  sulfide 
oxidation  in  bed,  g  02/(mVday) 

a  -  stoichiometric  conversion  for  mg/i  of  O2  uptake  per  mg/i 
of  substance  oxidized,  nondlmenslonal 

Therefore,  the  steady-state  additional  SOD  depends  on  the  loading  rate  from 
the  water  column  to  the  bed.  Water  column  oxygen  uptake  resulting  from  SOD  Is 
accomplished  In  the  DO  balance  by  multiplying  the  SOD  by  the  surface  area  of 
the  bed  and  dividing  by  the  volume  of  the  water  column,  which  is  the  same  as 
SOD/H  .  Therefore,  the  water  column  DO  sink  arising  from  Equation  39  is  com¬ 
puted  by  a(w/H)C  ,  where  w  for  reduced  manganese  Is  actually  V,  computed 
from  Equations  25  and  26 . 

98.  The  additional  sink  terms  arising  from  oxidation  reactions  of 
reduced  substances  that  were  Included  in  the  DO  balance  are  summarized  as 
follows : 


Dissolved  iron  = 

Dissolved  sulfide  =  a^K„(HS~) 
Reduced  manganese  =  (Mn**) 

Iron  sulfide  »  (FeS) 


(40) 


where  the  substance  quantities  In  parentheses  are  concentrations  In  milligrams 
per  liter.  The  default  values  for  the  parameters  in  the  above  relations  are 
shown  in  Table  6.  The  sources  of  DO,  reaeratlon  and  algal  photosynthesis, 
were  not  modified.  It  will  be  necessary  to  model  algae  In  many  tallwater 
applications  of  TWQH. 
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Table  6 

Default  Values  for  Stoichiometry  and  Kinetics 


Parameter 


Default  Value 


Op,,  DO  uptake  conversion  for  Fe"*^  oxidation 
Ohs,  do  uptake  conversion  for  HS~  oxidation 
Otto,  DO  uptake  conversion  for  Mn"*^  oxidation 
opaSi  uptake  conversion  for  FeS  oxidation 
Kp,,  oxidation  rate  of  Fe"*^ 

Kfis,  oxidation  rate  of  HS~ 

loss  rate  of  Mn'*'''  from  water  column 
Kp,s,  oxidation  rate  of  FeS  in  water  column 
V, ,  settling  rate  of  FeS  from  water  column 


0.14  mg/i  O2  per  mg/i  Fe'*'^ 
1.38  mg/i  O2  per  mg/i  S'^ 
0.29  mg/i  O2  per  mg/i  Mn'*'* 
0.73  mg/i  O2  per  mg/i  FeS 
(day"^)* 

(day*^)* 

(day*^)* 

38.0  day"^ 

1.0  m/day 


*  Parameter  predicted  by  the  model  for  default. 

Application 


General 

99.  TW(^  was  applied  to  one  data  set  for  each  of  the  four  field  sites 
(Lake  Greeson,  September  1987;  Nimrod  Reservoir,  August  1989;  Rough  River 
Reservoir,  August  1988;  and  Canyon  Reservoir,  September  1988)  discussed  in  Nix 
et  al.  (1991).  All  four  data  sets  were  obtained  with  the  steady-state  method 
of  sampling  (Nix  et  al.  1991).  Water  quality  constituents  modeled  during  each 
application  were  DO,  CBOD,  ammonium  nitrogen,  organic  nitrogen,  nitrate  nitro¬ 
gen,  dissolved  iron  (i.e.,  Fe'*’*) ,  dissolved  manganese  (i.e.,  Mn'*'*) ,  and  total 
dissolved  sulfide  (if  measured).  Iron  sulfide  was  also  modeled  for  Rough 
River.  Temperature  and  pH  were  not  modeled  but  were  input  for  each  reach. 

One  data  set  from  each  site  was  selected  for  model  verification. 

100.  Reaction  rates  for  Mn*^,  Fe*^,  and  NH*  removal  were  initially 
modeled  using  equations  discussed  in  previous  sections  of  this  report.  Con¬ 
centrations  of  Mn^^,  Fe"*^*,  and  NH*  were  compared  with  observed  results  to 
determine  how  well  removal  of  the  reduced  species  was  predicted  throughout  the 
tallwater.  If  predictions  did  not  compare  well  with  observations,  the  reac¬ 
tion  rate  was  changed  to  yield  concentrations  similar  to  those  observed.  All 
other  reaction  rates  (e.g.,  organic  nitrogen  hydrolysis)  were  set  to  values 
recommended  by  Brown  and  Barnwell  (1987) .  Table  7  contains  final  values  used 
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for  the  Input  parameters,  such  as  hydraulic  data,  headwater  boundary  con¬ 
ditions,  ambient  constants,  and  reaction  rates,  for  each  field  site. 

101.  The  following  section  presents  application  results.  The  order  of 
presentation  Is  to  discuss  the  results  for  each  water  quality  constituent  as 
modeled  for  all  sites  (e.g.,  results  of  modeling  Mn'*'^  for  all  sites),  rather 
than  to  discuss  the  results  of  all  water  quality  constituents  for  each  site. 
Results 

102.  Final  calibration  results  for  all  four  sites  are  shown  In  Fig¬ 

ures  7  - 10 .  Final  results  for  are  sho%m  in  Figures  7b  (Greeson) ,  8b 
(Nimrod) ,  9b  (Rough  River) ,  and  10b  (Canyon) .  Mn'*’*  removal  was  successfully 

modeled  at  three  of  the  sites  using  the  predictive  equations  for  K(ti  (l.e.. 
Equations  22-26),  resulting  In  acceptable  estimations  of  Mn'*'^  concentrations 
compared  with  observed  values.  However,  for  the  Rough  River  application,  the 
predictive  equations  for  K^n  caused  Mn'''^  to  be  removed  too  quickly,  espe¬ 
cially  in  the  upper  reaches.  The  first  four  reaches  (approximately  the  upper 
10.0  km)  of  the  modeled  Rough  River  tailwater  are  in  the  pool  created  by  the 
mill  dam,  referred  to  as  Falls  of  Rough.  Conditions  within  this  pool  are  much 
different  from  the  downstream  reach  below  Falls  of  Rough.  By  specifying  low 
Kmt,  values  that  are  consistent  with  those  estimated  within  the  mill  dam  pool, 
much  better  results  were  obtained  (Figure  9b).  During  data  analysis 

(Part  III) ,  it  was  observed  that  the  Kf*,  values  measured  in  the  pooled 
reaches  of  Rough  River  were  much  less  than  those  measured  at  the  other  sites. 

A  higher  value  of  Kj*,  ,  similar  to  values  observed  at  the  other  sites,  was 
used  in  the  free -flowing  reach  below  the  mill  dam. 

103.  The  predictive  equations  for  the  manganese  removal  rate  1^ 
probably  did  not  perform  well  for  the  pooled  reaches  of  the  Rough  River  tail- 
water  because  of  the  conditions  that  normally  exist  during  the  summer  months 
at  this  site.  The  releases  from  Rough  River  Dam  are  usually  quite  low  (e.g., 
20  cfs*) .  A  low-flow  and  a  high-flow  condition  were  sampled  in  this  study. 

The  results  presented  here  are  for  the  high  flow  of  200  cfs.  For  this  flow, 
the  water  in  the  pool  was  well  mixed  from  surface  to  bottom.  However,  at  the 
normal  low  flows,  the  Fall  of  Rough  pool  slightly  stratifies,  and  the  release 
water  plunges  under  a  warmwater  wedge.  This  plunging  phenomenon  prevents  the 
release  water,  which  is  rich  in  reduced  substances  (e.g.,  dissolved 


*  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  3. 
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Table  7 


Greeson  Hydraulic  Depth  (ft) 

Velocity  (ft/sec) 
Start  of  reach  (km) 
Slope 


3.11  2.50 

1.65  1.50 

170.27  166.41 

0.0016 


2.80  2.63 

0.88  0.75 

161.26  159.65 


Headwater  Flow  (cfs)  500.00 

DO  (mg/1)  4.40 

CBOD  (mg/1)**  2.00 

Fe**  (mg/1)  0.04 

Mn+*  (mg/1)  0.09 

pH  6.15 

Org-N  (mg/1)  0.18 

NH4-N  (mg/1)  0.02 

NO3-N  (mg/1)  0.17 


Ambient 

Temperature  (*C) 

14.04 

15.40 

16.00 

15.50 

constants 

pH 

6.15 

6.40 

6.30 

6.80 

Reaction 

CBOD  (1/day) 

0.10 

0.10 

0.10 

0.25 

rates 

SOD  (^m*  day‘‘) 

0.00 

0.00 

0.00 

2.0 

K2  value 

6.24 

7.03 

5.62 

0.62 

Org-N  hydrolysis  (1/day) 

0.20 

Org-N  settling  (1/day) 

0.10 

NH4-N  (1/day) 

Modeled 

Fe**  oxidation  (1/day) 

Modeled 

Mn**  oxidation  (1/day) 

Modeled 

Denitrification  (1/day) 

0.00 

Nimrod  Hydraulic 

Depth  (ft) 

2.00 

2.56 

3.57 

Velocity  (ft/sec) 

0.56 

0.45 

0.27 

Start  of  reach  (km) 

100.74 

97.85 

96.08 

Slope 

0.001 

Headwater  Temperature  (*C)  25.88 


Flow  (cfs) 

200.00 

DO  (mg/i) 

7.00 

CBOD  (mg/i) 

1.00 

Sulfide  (mg/i) 

0.016 

Fe**  (mg/i) 

1.00 

Mn'^*  (mg/i) 

1.20 

pH 

6.50 

Org-N  (mg/i) 

0.32 

NH*-N  (mg/i) 

0.28 

NO3-N  (mg/i) 

0.03 

Ambient 

Temperature  (*C) 

25.88 

constants 

pH 

6.50 

Reaction 

CBOD  (1/day) 

0.15 

rates 

SOD  (g/m*  day'*) 

2.00 

K2  option 

3.00 

Org-N  hydrolysis  (1/day) 

0.2 

Org-N  settling  (1/day) 

0.1 

NH4-N  (1/day) 

1.75 

(Continued) 


*  Parameters  were  constant  throughout  reaches  unless  otherwise  Indicated. 

**  CBOD  Is  5  day  CBOD.  (Sheet  1  of  3) 
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Table  7  (Continued) 


Field 

Site 

Nimrod 
(Cent. ) 


Rough 

River 


Canyon 


Reach 


Data  Tvoe 

Parameter 

1 

2 

3 

4  _ 

5 

Fe**  (l/day) 

Hodeled 

Mn+*  (l/day) 

Modeled 

Sulfide  (l/day) 

1.00 

0.50 

0.10 

Denitrification  (l/day) 

2.25 

Hydraulic 

Depth  (ft) 

5.50 

6.00 

9.00 

10.70 

3.00 

Velocity  Cft/sec) 

0.57 

0,30 

0.26 

0.33 

0.62 

Start  of  reach  (km) 

143 . 88 

141.14 

138.89 

135.67 

133.58 

Slope 

0.0023 

Headwater 

Flow  (c£s) 

200.00 

DO  (mg/i) 

8.60 

CBOD  (mg/i) 

1.00 

Fe**  (mg/f) 

0.89 

FeS  (mg/f) 

2.50 

Mn*^*  (mg/f) 

2.03 

pH 

7.30 

Org-N  (mg/f) 

0.64 

NH*-N  (mg/f) 

1.06 

NOj-N  (mg/f) 

0.01 

Ambient 

Temperature  (*C) 

15.20 

15.30 

15.40 

16.50 

17.30 

constants 

pH 

7.30 

7.30 

7.20 

7.30 

7.20 

Reaction 

CBOD  (l/day) 

0.20 

rates 

SOD  (g/m*  day"') 

1.80 

Kj  option 

3.00 

Org-N  hydrolysis  (l/day) 

0.20 

Org-N  settling  (l/day) 

0.10 

NH^-N  (l/day) 

0.20 

0.20 

0.20 

0.20 

0.70 

Fe**  (l/day) 

Modeled 

FeS  (l/day) 

2.00 

FeS  settling  (m/day) 

0.00 

Mn**  (l/day) 

0.05 

0.05 

0.05 

0.03 

0.50 

Denitrification  (l/day) 

3.50 

Hydraulic 

Depth  (ft) 

1.30 

Velocity  (ft/sec) 

1.04 

Start  of  reach  (km) 

487.63 

486.02 

484.73 

482.32 

477.65 

Slope 

0.00144 

Headwater 

Flow  (cfs) 

114.00 

DO  (mg/f) 

8.80 

CBOD  (mg/f) 

1.00 

Sulfide  (mg/f) 

0.04 

Fe**  (mg/f) 

0.012 

Mn^*  (mg/f) 

0.168 

pH 

7.50 

Org-N  (mg/f) 

0.00 

NH*-N  (mg/f) 

O.ll 

NOj-N  (mg/f) 

0.04 

Ambient 

Temperature  (*C) 

17.20 

17.70 

19.10 

20.60 

22.50 

constants 

pH 

7.50 

7.70 

7.80 

7.80 

7.90 

(Continued) 


(Sheet  2  of  3) 


Table  7  (Concluded) 


Field 


Canyon  Reaction  CBOD  (1/day)  0.15 

(Cont.)  rates  SOD  (g/m*  day**)  2.00 

K2  option  3.00 

Org-N  hydrolysis  (1/day)  0.20 

Org-N  settling  (1/day)  0.10 

NH*-N  (1/day)  1.80 

Fe+*  (1/day)  Modeled 

Mn**  (1/day)  Modeled 

Sulfide  (1/day)  25.00 

Denitrification  (1/day)  0.00 


(Sheet  3  of  3) 
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AMMONIA  NITROGEN.  mg/I  DISSOLVED  OXYGEN.  mg/I 


AMMONIA  NITROGEN.  mg/I  DISSOLVED  OXYGEN,  mg/I 


a.  Dissolved  oxygen 


c.  Anunonium  nitrogen 


b.  Dissolved  manganese 


RIVER  KILOMETERS 


d.  Organic  nitrogen 


Figure  8.  Predicted  and  observed  water  quality, 
Nimrod  Dam  tailwater  (Continued) 


NITRATE  NITROGEN.  mg/I 


NITRATE  NITROGEN,  mg/I  AMMONIA  NITROGEN,  mg/I  DISSOLVED  OXYGEN.  mg/I 


AMMONIA  NITROGEN,  mg/I 


490  485  480  475  470  465 

RIVER  KILOMETERS 

a.  Dissolved  oxygen 


RIVER  KILOMETERS 
c .  Anunonlum  nitrogen 


490  485  480  475  470  465 

RIVER  KILOMETERS 

b .  Dissolved  manganese 


d.  Nitrate  nitrogen 


Figure  10.  Predicted  and  observed  water  quality, 
Canyon  Dam  tallwater 


manganese),  from  being  reaerated  through  the  air-water  interface.  For  these 
conditions,  it  seems  reasonable  to  expect  that  there  is  less  opportunity  for 
hydrous  oxides  of  manganese  to  accumulate  on  the  substrate,  thus  lowering  the 
capacity  of  the  pooled  reaches  to  remove  manganese. 

104.  Final  results  for  ammonium  are  shown  in  Figures  7c  (Greeson) ,  8c 
(Nimrod) ,  9c  (Rough  River) ,  and  10c  (Canyon) .  The  ammonia  oxidation  rate 
(l.e.,  nitrification  rate)  was  initially  modeled  using  Equation  27  for  all 
four  sites.  However,  results  at  Nimrod,  Rough  River,  and  Canyon  (Figures  11a, 
12c,  and  13,  respectively)  showed  that  the  nitrification  rate  was  calculated 
to  be  too  high  (Rough  River)  or  too  low  (Nimrod  and  Canyon) .  Modeled  nitrifi¬ 
cation  rates  (i.e..  Equation  27)  were  retained  at  the  Greeson  site,  which  had 
low  concentrations  (Figure  7c).  The  nitrification  rates  were  calibrated 
against  observed  ammonia  nitrogen  data,  yielding  the  rates  shown  in  Table  7 
and  the  concentrations  shown  in  Figures  7c,  8c,  9c,  and  10c.  The  low  rates 
used  for  the  pooled  reaches  of  Rough  River  are  probably  related  to  the  reasons 
discussed  above  for  manganese . 

105.  The  inability  to  predict  the  nitrification  rate  at  three  of  the 
four  sites  may  be  due  to  the  scatter  in  the  data  used  to  develop  Equation  27. 
Nitrification  rates  also  depend  upon  bacteria  and  conditions  for  bacterial 
growth,  thus  making  prediction  of  rates  difficult.  However,  the  predicted 
nitrification  rates  were  generally  close  to  the  calibrated  rates  with  only 
slight  adjustments  required  during  calibration. 

106.  Results  for  organic  nitrogen  are  shown  in  Figures  7d  (Greeson),  8d 
(Nimrod) ,  and  9d  (Rough  River) .  Boundary  conditions  for  organic  nitrogen 
(Table  7)  were  estimated  as  the  difference  between  total  Kjeldahl  nitrogen 
(TKN)  and  ammonia  nitrogen  for  each  site  (Nix  et  al.  1991).  TKN  was  not  mea¬ 
sured  at  Canyon  tailwater,  so  it  was  not  modeled  at  this  site.  Default  values 
(Brown  and  Barnwell  1987)  for  organic  nitrogen  hydrolysis  and  settling  were 
used  for  all  applications  (Table  7).  Final  results  for  organic  nitrogen  com¬ 
pared  favorably  with  observed  concentrations  except  at  Nimrod,  where  the  con¬ 
centrations  increased  downstream  instead  of  decreasing  as  would  be  expected. 
This  indicates  that  there  may  be  a  source  of  organic  nitrogen,  such  as  algae, 
that  is  not  being  accounted  for  in  the  modeling  effort. 

107.  Results  for  nitrate  nitrogen  are  shown  in  Figures  7e  (Greeson),  8e 
(Nimrod),  9e  (Rough  River),  and  lOd  (Canyon).  Initially,  TWQM  had  no  mecha¬ 
nism  for  denitrification,  and  the  results  for  nitrate  nitrogen  did  not  compare 
favorably  with  observed  values,  especially  at  Rough  River  (see  Figure  12d) . 
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a.  Ammonium  nitrogen 


b.  Dissolved  sulfide 


Figure  11.  Initial  predictions  versus 
observed  water  quality,  Nimrod  Dam 
tallwater 
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AMMOMIA  NITROGEN.  mg/I 


125 


US  140  135  130 

RIVER  KILOMETERS 

a.  Dissolved  oxygen 


RIVER  KILOMETERS 
c .  Ammonium  nitrogen 


Figure  12.  Initial  predictions 

Rough  River 


NITRATE  NITROGEN.  mg/I 


145  140  135  130 

RIVER  KILOMETERS 

b.  Dissolved  manganese 


d.  Nitrate  nitrogen 


observed  water  quality, 


RIVER  KILOMETERS 


Figure  13.  Initial  prediction 
versus  observed  ammonium  nitro¬ 
gen,  Canyon  Dam  tailwater 

Consequently,  denitrification  was  added  as  a  loss  mechanism  for  nitrate  nitro¬ 
gen  in  TUQM.  Once  denitrification  was  included  in  the  model,  the  final 
results  (Figures  8e,  Nimrod,  and  9e,  Rough  River)  for  nitrate  nitrogen  were 
improved  from  the  initial  applications.  Greeson  and  Canyon  tailwaters  did  not 
exhibit  denitrification,  which  could  be  related  to  the  stream  sediment  type. 
Both  sites  have  cobble  beds  where  the  anoxic  sediment  zone  necessary  for  deni¬ 
trification  could  be  farther  below  the  sediment -water  interface.  According  to 
Hill  and  Sanmugadas  (1985),  denitrification  rates  are  correlated  to  stream- 
sediment  characteristics. 

108.  Final  results  for  dissolved  sulfide  are  shown  in  Figures  8f 
(Nimrod)  and  lOe  (Canyon).  Nimrod  and  Canyon  were  the  only  sites  where  sul¬ 
fide  was  successfully  measured,  although  the  odor  of  sulfide  wa^  detectable  at 
the  other  sites.  Initial  results  for  sulfide  at  Nimrod  (Figure  11b)  showed 
sulfide  oxidizing  too  quickly  through  the  model  study  reach.  By  reducing  the 
sulfide  oxidation  rate  more  than  an  order  of  magnitude  from  the  default  value 
(25  day'^) ,  a  better  comparison  to  the  observed  data  was  obtained  (Figure  8f ) . 
However,  the  predicted  sulfide  approaches  zero  concentration  at  the  end  of  the 
study  reach,  while  the  observed  data  appear  to  remain  at  0.01  mg/i.  As  men¬ 
tioned  in  Part  III,  the  data  inferred  that  the  oxidation  rate  of  HS~  at  Nimrod 
was  on  the  low  side  of  reported  ranges.  Also,  most  of  the  sulfide  for  the  pH 
conditions  at  Nimrod  is  H2S,  which  volatilizes.  Adjustments  to  the  reaeration 
rate  K2  may  have  Improved  sulfide  predictions.  Sulfide  was  measured  at  the 
first  station  at  Canyon  (40  /ig/i)  but  declined  to  near  zero  at  the  other 
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stations  (Nix  et  al.  1991).  The  default  value  for  sulfide  oxidation  was  used 
at  this  site,  and  the  results  compared  favorably  with  observed  data 
(Figure  lOe). 

109.  Final  results  for  DO  are  shown  In  Figures  7a  (Greeson) ,  8a 
(Nimrod) ,  9a  (Rough  River) ,  and  10a  (Canyon) .  The  only  model  source  of  DO  for 
each  application  was  stream  reaeratlon.  Reaeratlon  Is  especially  high  at  the 
Greeson  site  (Figure  7a) ,  where  DO  releases  were  low  due  to  this  site  operat¬ 
ing  for  hydropower.  For  all  applications  except  the  Greeson  site,  stream 
reaeratlon  was  calculated  using  the  O' Connor -Dobbins  (1958)  equation  option  In 
TVQM.  Measured  stream  reaeratlon  rates  at  Greeson  were  obtained  by  Wilhelms 
(1987)  using  a  hydrocarbon  gas  tracer  method  (Rathbun  et  al.  1977,  1978). 
However,  the  Tslvoglou  and  Wallace  (1972)  formulation  for  reaeratlon  (l.e., 
the  energy  dissipation  method)  was  also  found  to  yield  good  results  at  Greeson 
(Dortch  and  Hamlin  1988) . 

110.  The  sinks  of  DO  for  all  applications  Included  SOD  and  CBOD. 
Although  Mn*^,  Fe^^,  sulfide,  and  ammonia  oxidation  are  also  sinks  of  DO,  they 
had  little  Individual  effect  on  DO  concentrations,  but  their  cumulative  effect 
was  not  negligible.  The  low  Individual  effects  were  due  to  the  relatively  low 
release  concentrations  at  all  the  sites. 

111.  Boundary  conditions  for  CBOD  measurements  were  estimated  from  the 
10-day  BOD  studies  of  In-pool  samples  conducted  by  Nix  et  al.  (1991).  The 
deepest  In-pool  sample  was  used  to  estimate  release  CBOD  for  all  sites.  The 
CBOD  was  estimated  by  subtracting  accountable  oxygen  demands  (l.e.,  oxidation 
of  reduced  substances)  from  the  total  oxygen  demand.  The  decrease  In  reduced 
substances  and/or  the  Increase  In  the  oxidized  forms,  such  as  sulfate,  was 
converted  to  DO  equivalents  through  stoichiometric  relationships.  Estimated 
release  concentrations  for  CBOD  at  all  sites  were  approximately  1  mg/i  except 
Roug^  River,  which  was  2  mg/i.  SOD  rates  were  calibrated  (Table  7)  at  Greeson 
and  Canyon;  measured  SOD  values  (Nix  et  al.  1991)  were  used  for  Rough  River 
and  Nimrod.  Zero  SOD  was  used  at  the  Greeson  site  for  all  reaches  except  the 
last  reach  due  to  the  rocky,  cobble  substrate  conditions.  The  last  reach  Is 
characterized  as  a  slower  moving  stream  with  silty  substrate  and  a  flatter 
slope;  thus,  SOD  was  specified  In  this  reach. 

112.  With  major  sources  and  sinks  of  DO  accounted  for,  all  the  pre¬ 
dicted  DO  measurements  compared  favorably  with  the  observed  data  except  at 
Rough  River  (Figure  12a) ,  where  data  Indicated  there  may  be  another  sink  of  DO 
not  accounted  for  In  the  modeling  effort.  Based  on  examination  of  the  total 
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and  dissolved  Iron  measurements,  the  data  indicated  that  FeS  might  be  present, 
especially  since  sulfate  Increased  In  the  dovmstream  direction.  Sulfate  also 
Increased  with  time  in  the  10-day  BOD  studies  for  Rough  River.  The  decision 
was  made  to  include  FeS  as  a  model  variable  at  this  site.  The  FeS  boundary 
conditions  were  estimated  as  the  difference  between  total  and  dissolved  iron, 
yielding  a  value  for  FeS  of  2.5  mg/i.  Final  results  of  Rough  River  DO  are 
shown  in  Figure  9a  and  compare  favorably  with  observed  data.  The  reaeration 
through  the  mill  dam  is  also  evident  in  Figure  9a;  this  reaeration  was  modeled 
using  the  equation  of  Gameson  (Brown  and  Barnwell  1987). 

113.  The  other  data  sets  were  examined  for  the  possible  presence  of 
FeS,  but  only  the  Nimrod  site,  in  addition  to  Rough  River,  indicated  the  pos¬ 
sible  presence  of  FeS.  A  substantial  difference  in  total  and  dissolved  iron 
occurred  at  Nimrod,  which  could  be  attributed  to  FeS.  However,  there  was  no 
Increase  in  sulfate  in  the  tailwater.  If  FeS  was  present,  it  must  have  oxi¬ 
dized  rapidly  as  it  passed  through  the  release  structure  and  the  first  200  m 
of  the  tailwater. 

114.  Dissolved  iron  (i.e.,  Fe'*'^)  was  modeled  for  Nimrod  and  Rough  River 
(Figures  8g  and  9f,  respectively);  very  little  iron  was  detected  at  the  Gree- 
son  and  Canyon  sites.  The  model  shows  reduced  iron  decreasing  rapidly  to 
values  near  zero,  whereas  the  data  indicate  that  dissolved  iron  decreases  much 
more  slowly  and  does  not  reach  zero.  Nimrod  data  actually  indicate  that 
dissolved  iron  does  not  decrease  after  traveling  through  the  first  few  kilo¬ 
meters  of  tailwater.  As  discussed  earlier,  it  is  suspected  that  fine  particu¬ 
late  iron  passed  through  the  filter,  thus  yielding  an  inaccurate  measurement 
of  dissolved  iron. 
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PART  V:  CONCLUSIONS  AND  RECOMMENDATIONS 


115 .  The  studies  reported  herein  were  an  attempt  to  obtain  an  Improved 
understanding  of  water  quality  processes  in  reservoir  tallwaters.  The  focus 
was  on  DO  and  the  removal  and/or  oxidation  of  reduced  substances  (e.g.,  dis¬ 
solved  iron  and  manganese ,  sulfide ,  and  ammonia)  present  in  deep  releases  from 
anoxic  reservoir  strata.  The  intent  of  the  analyses  was  to  develop  kinetic 
descriptions  and  rates  for  constructing  a  model  of  water  quality  in  reservoir 
tailwaters.  An  attempt  was  made  to  develop  relationships  to  predict  kinetic 
rate  coefficients  based  on  ambient  stream  conditions. 

116 .  Much  information  was  gained  from  the  dissolved  manganese  and 
£umnonla  data  collected  in  this  study.  It  was  possible  to  develop  relation¬ 
ships  to  predict  the  removal  rate  of  dissolved  manganese  through  adsorption 
onto  the  streambed.  These  relationships  are  based  on  the  type  of  streambed 
(i.e.,  sediment  or  cobble)  and  the  stream  hydraulic  conditions,  which  affect 
the  diffusive  sublayer  thickness  and  the  bed/water  column  mass  transfer  rate. 
Similarly,  a  relationship  for  the  nitrification  rate  was  developed  using 
stream  hydraulic  conditions.  The  predictive  relationship  for  nitrification 
was  not  as  strongly  correlated  to  hydraulic  conditions  as  the  relationship  for 
manganese  removal.  Both  relationships  were  developed  from  results  of  four 
field  study  sites;  thus,  only  through  studies  at  other  sites  can  the  relation¬ 
ships  be  further  improved. 

117.  Data  limitations  precluded  the  development  and/or  improvement  of 
predictive  relationships  for  reduced  iron,  sulfide,  and  iron  sulfide  oxidation 
rates.  Methods  for  measuring  reduced  (i.e.,  dissolved)  iron  and  dissolved 
sulfide  were  not  adequate  to  accurately  account  for  the  oxidation  and/or 
removal  of  these  substances  from  the  tailwater.  Additionally,  the  data  sug¬ 
gest  that  iron  sulfide  was  present  in  the  releases  at  several  of  the  sites, 
further  complicating  the  studies.  A  limited  amount  of  data  was  used  to  assess 
existing  equations  for  predicting  iron  and  sulfide  oxidation  rates.  These 
results,  although  far  from  conclusive,  did  suggest  that  existing  rate  equa¬ 
tions  could  be  used  as  a  rough  estimate  for  oxidation  rates. 

118.  Future  studies  should  further  address  the  oxidation/removal 
kinetics  of  iron,  sulfide,  and  iron  sulfide.  Proper  methods  for  sampling, 
handling,  and  analysis  must  be  prescribed  to  better  quantify  concentrations  of 
ferrous  and  ferric  iron  and  iron  sulfide.  For  example,  Nix  et  al.  (1991) 
recommended  that  a  0.1-/i  (rather  than  O.Ub-fi)  filter  be  used  to  separate 
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dissolved  and  particulate  Iron.  The  sulfide  analysis  technique  of  Hongzhang 
(1982),  which  was  used  during  the  latter  stages  of  the  field  studies  by  Nix  et 
al.  (1991),  shows  much  promise  for  measuring  low  concentrations  of  sulfide.  A 
technique  must  be  Identified  to  distinguish  the  amount  of  particulate  Iron 
that  Is  Iron  sulfide  versus  ferric  Iron.  After  prescribing  Improved  sampling 
and  analysis  procedures,  field  studies  similar  to  the  ones  reported  In  this 
study  should  be  undertaken  to  yield  Information  concerning  kinetic  rates . 
Laboratory  kinetic  rate  studies  would  also  greatly  expedite  the  development  of 
general  predictive  relationships  for  processes  that  are  not  dependent  on 
stream  substrate  and  hydraulics. 

119.  The  results  of  the  10-day  BCD  studies  (Nix  et  al.  1991)  Indicate 
that  little  of  the  oxygen  demand  of  hypolimnetic  water  is  caused  by  the  decom¬ 
position  of  organic  carbon.  However,  these  results  were  not  conclusive 
because  of  the  analytical  problems  discussed  above  for  iron  and  sulfide.  I£ 
future  studies  are  conducted  to  investigate  iron  and  sulfide,  10-day  BOD 
studies  should  be  conducted  to  more  accurately  delineate  any  DO  demand  arising 
from  organic  carbon  synthesis.  These  studies  should  be  conducted  with  and 
without  seeded  bacteria.  Some  consideration  should  also  be  given  to  detecting 
methane  and  any  associated  DO  demand. 

120.  The  model  developed  in  this  work  (TWQM)  is  considered  to  have 
performed  well  for  DO,  dissolved  manganese,  nitrogen  species,  and  dissolved 
sulfide.  The  relationships  for  predicting  dissolved  manganese  removal  rates 
are  considered  to  be  applicable  for  tailwaters  that  have  experienced  deep 
releases  containing  reduced  manganese.  For  sites  that  do  not  have  such 
releases  but  are  considering  a  deep  release  (such  as  opening  a  bottom  sluice 
that  Is  rarely  used) ,  it  may  take  some  time  (on  the  order  of  weeks)  for  the 
tailwater  to  acclimate  and  form  a  substrate  that  effectively  sorbs  dissolved 
manganese.  The  predictive  equations  for  manganese  removal  do  not  account  for 
any  acclimation  period. 

121.  Although  the  data  for  hydrogen  sulfide  were  limited,  the  model 
should  provide  reasonable  estimates  since  the  formulations  are  based  in  part 
on  volatilization,  which  is  fairly  well  understood.  Although  the  model  con¬ 
tains  simplistic  algorithms  for  dissolved  iron  and  iron  sulfide,  the  lack  of 
adequate  field  data  prohibited  the  development  of  general  relationships  for 
kinetic  rates  or  proper  evaluation  of  model  results.  As  discussed  above,  more 
work  Is  required  to  better  resolve  the  kinetics  of  dissolved  iron  and  iron 
sulfide . 
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122.  It  is  recommended  that  future  extensions  of  the  model  include  the 
capability  to  predict  pH  changes  in  the  tailwater.  This  could  be  accomplished 
by  building  in  CO2,  alkalinity,  and  carbonate  equilibria.  The  escape  of  CO2 
through  the  air-water  Interface  can  be  modeled  like  other  gas  exchange  pro¬ 
cesses,  such  as  reaeration  and  volatilization  of  H2S.  As  CO2  escapes,  pH 
increases.  Since  pH  affects  other  processes,  such  as  the  rate  of  oxidation  of 
reduced  iron,  incorporating  pH  would  advance  the  utility  of  the  model. 
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APPENDIX  A:  1989  NIMROD  TATLWATER  DATA 


A1 


Table  A1 

Lake  Nimrod  High  Flow  (200  cfs) 


Temp. 

DO 

Specific 

Station 

Time 

”C 

mp/i 

eh_ 

Conductance . 

us /cm* 

Aueust  23. 

1989 

A 

1634 

26.0 

7.0 

6.5 

51 

B1 

1615 

27.0 

26.2 

6.5 

48 

B2 

1548 

25.5 

6.8 

6.6 

48 

B3 

1519 

27.0 

6.7 

6.5 

48 

C 

1443 

28.0 

6.7 

6.5 

46 

D 

1405 

26.5 

5.9 

6.5 

46 

August  24. 

1989 

A 

1045 

29.0 

7.0 

6.5 

49 

B1 

1015 

29.0 

6.7 

6.5 

50 

B2 

955 

27.0 

6.6 

6.5 

48 

B3 

920 

27.0 

6.5 

6.5 

48 

C 

815 

22.5 

6.1 

6.5 

47 

D 

815 

21.5 

5.8 

6.4 

46 

Turbidity 

Alkalinity  CO2 

TOC 

TKN 

NH4-N 

N03N 

Station 

_ lim _ 

- bb/1 

-  mg/i 

m&ll 

mg/i 

mgJA 

mg/i 

August  23. 

1989 

A 

36 

16 

6 

7. 

0.6 

0.28 

0.03 

B1 

24 

19 

5 

7. 

0.5 

0.22 

0.04 

B2 

23 

14 

5 

5. 

0.6 

0.20 

0.07 

B3 

23 

14 

5 

6. 

0.6 

0.18 

0.07 

C 

22 

15 

4 

6. 

0.6 

0.14 

0.10 

D 

16 

14 

6 

5, 

0.5 

0.07 

0.14 

Aueust  24. 

1989 

A 

29 

16 

5 

6. 

0.7 

0.26 

0.02 

B1 

30 

16 

5 

6. 

0.8 

0.25 

0.03 

B1  (R) 

27 

16 

6 

6. 

0.6 

0.26 

0.04 

B2 

28 

16 

5 

6. 

0.7 

0.25 

0.08 

B3 

21 

15 

6 

6. 

0.7 

0.22 

0.08 

C 

19 

15 

6 

6. 

0.5 

0.19 

0.11 

D 

18 

14 

6 

5. 

0.5 

0.10 

0.15 

(Continued) 


*  fis  -  fi  mhos/cm. 
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Table  A1  (Concluded) 


Station 

S‘2 

Hill 

SOi* 

as2i 

ci‘ 

fflS/1 

Total  Fe 
mg/i 

Diss.  Fe 
_ mg/i _ 

Total  Mn 
mg/f 

Diss.  1 
mg/i 

August  23. 

1989 

A 

16 

2.7 

2.0 

2.6 

1.0 

1.30 

1.2 

Bl 

17 

2.9 

1.9 

2.1 

0.8 

1.00 

1.0 

B2 

11 

2.9 

2.0 

1.9 

0.7 

0.90 

0.8 

B3 

10 

3.0 

2.1 

1.7 

0.7 

0.88 

0.7 

C 

9 

3.2 

2.0 

1.6 

0.7 

0.71 

0.6 

D 

10 

3.0 

2.1 

1.3 

0.7 

0.41 

0.3 

August  24. 

1989 

A 

17 

2.7 

1.9 

2.4 

0.9 

1.20 

1.1 

Bl 

12 

3.2 

2.0 

2.2 

0.9 

1.20 

1.1 

Bl  (R) 

11 

2.7 

2.0 

2.2 

0.9 

1.20 

1.1 

B2 

11 

2.8 

2.0 

2.1 

0.8 

1.10 

1.0 

B3 

10 

1.9 

3.0 

1.9 

0.8 

0.95 

0.8 

C 

11 

2.9 

2.0 

1.8 

0.8 

0.82 

0.7 

D 

10 

2.8 

2.0 

1.4 

0.8 

0.42 

0.3 

Note:  B1  (R)  is  a  repeated  sampling  of  station  Bl. 


A4 


Table  A2 

Lake  Nimrod  Profiles 


Temperature  DO  Specific 

Depth,  ro  °C _  ffig/1  pH  Conductance,  us/cm 

Aupust  23.  1989.  0745  hr 


0 

28.3 

5.4 

6.7 

40 

1 

28.3 

5.3 

6.6 

40 

2 

28.1 

4.1 

6.6 

40 

3 

27.9 

3.3 

6.6 

40 

4 

27.4 

2.2 

6.4 

41 

5 

26.6 

0.7 

6.3 

42 

6 

26.2 

0.1 

6.2 

44 

7 

26.0 

0.0 

6.2 

47 

8 

25.9 

0.0 

6.2 

49 

9 

25.7 

0.0 

6.2 

55 

10 

25.0 

0.0 

6.2 

64 

Aueust  23.  1989. 

1845  hr 

0 

28.0 

6.4 

6.8 

41 

1 

28.2 

4.4 

6.7 

41 

2 

27.6 

2.1 

6.5 

41 

3 

26.5 

0.2 

6.4 

44 

4 

26.3 

0.1 

6.3 

46 

5 

26.2 

0.0 

6.3 

46 

6 

26.1 

0.0 

6.2 

47 

7 

26.0 

0.0 

6.2 

52 

8 

25.7 

0.0 

6.2 

61 

9 

25.5 

0.0 

6.2 

61 

10 

24.8 

0.0 

6.2 

71 

Aueust  24.  1989. 

0650  hr 

0 

27.2 

3.0 

6.4 

42 

1 

27.2 

3.1 

6.4 

42 

2 

27.2 

3.0 

6.4 

42 

3 

27.2 

2.8 

6.3 

42 

4 

26.8 

0.6 

6.2 

42 

5 

26.2 

0.0 

6.2 

46 

6 

26.0 

0.0 

6.2 

49 

7 

25.9 

0.0 

6.2 

50 

8 

25.8 

0.0 

6.1 

55 

9 

25.6 

0.0 

6.1 

57 

10 

24.9 

0.0 

6.2 

67 

(Continued) 
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Table  A2  (Concluded) 


Depth,  m 

Temperature 

"C 

DO 

l&Zl 

Aueust  24.  1989. 

7.0 

_EH 

Specific 

Conductance.  «s/cm 

0 

31.2 

1400  hr 

6.8 

41 

1 

30.9 

6.8 

6.8 

41 

2 

30.1 

6.4 

6.7 

41 

3 

28.2 

2.9 

6.5 

42 

4 

27.8 

2.9 

6.1 

42 

5 

27.5 

0.2 

6.1 

41 

6 

26.6 

0.0 

6.1 

47 

7 

26.4 

0.0 

6.1 

50 

8 

26.3 

0.0 

6.1 

51 

9 

26.2 

0.0 

6.1 

53 

10 

25.7 

0.0 

6.1 

61 
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Table  A3 


Sample 


N-1-0 

N-1-1 

N-1-2 

N-1-3 

N-1-4 

N-1-5 

N-1-5R 

N-1-6 

N-1-7 

N-1-8 

N-1-9 

N-1-10 


Sample 


N-1-0 

N-1-1 

N-1-2 

N-1-3 

N-1-4 

N-1-5 

N-1-5R 

N-1-6 

N-1-7 

N-1-8 

N-1-9 

N-1-10 


Lake  Nimrod  In-pool  Water  Quality 


Free 

Total 

Dissolved 

Total 

Dissolved 

Lab 

CO2 

Fe 

Fe 

Mn 

Mn 

S'2 

pH 

mp/i 

mc/i 

me/i 

me/i 

me/i 

ae/i 

Aueust 

23.  1989 

6.64 

4 

0.82 

0.41 

0.09 

0.00 

6 

6.60 

4 

0.79 

0.35 

0.09 

0.01 

7 

6.60 

6 

0.84 

0.45 

0.10 

0.01 

7 

6.54 

7 

0.93 

0.47 

0.12 

0.01 

6 

6.47 

8 

0.95 

0.48 

0.14 

0.02 

8 

6.33 

10 

1.4 

0.57 

0.45 

0.30 

9 

6.33 

9 

1.4 

0.55 

0.38 

0.24 

10 

6.36 

10 

1.5 

0.57 

0.47 

0.30 

10 

6.41 

11 

2.0 

0.48 

0.90 

0.79 

10 

6.39 

10 

2.1 

0.43 

1.1 

1.0 

14 

6.40 

12 

2.6 

0.90 

1.5 

1.4 

14 

6.39 

13 

3.2 

1.4 

1.9 

1.8 

30 

nh;" 

NO3-N 

so;2 

c/' 

TOC 

DOC 

TKN 

Alkalinity 

Turbidity 

mp/i 

me/i 

me/i 

mg/i 

me/i 

me/i 

me/i 

me/i 

NTU 

Aueust 

23.  1989 

0.02 

0.03 

3.90 

2.2 

4.4 

3.1 

0.4 

11 

0.03 

0.02 

2.89 

2.0 

4.2 

3.0 

0.8 

11 

10.0 

0.05 

0.03 

2.99 

2.1 

4.2 

3.2 

0.6 

11 

11.0 

0.08 

0.03 

3.20 

2.2 

4.1 

3.3 

0.6 

11 

14.0 

0.08 

0.03 

3.36 

2.1 

5.3 

3.3 

0.5 

11 

11.0 

0.14 

0.03 

3.98 

2.0 

4.4 

3.6 

0.8 

12 

20.0 

0.13 

0.02 

3.64 

1.9 

4.5 

3.4 

0.8 

12 

21.0 

0.13 

0.02 

3.10 

2.0 

4.1 

3.2 

0.6 

12 

21.0 

0.20 

0.02 

2.94 

2.0 

4.2 

3.4 

0.7 

14 

31.0 

0.23 

0.05 

3.54 

2.4 

3.9 

3.8 

0.7 

14 

30.0 

0.32 

0.04 

3.27 

2.4 

4.0 

3.4 

0.7 

17 

35.0 

0.42 

0.06 

2.75 

2.5 

4.7 

3.3 

0.08 

18 

39.0 

APPENDIX  B 


:  MODEL  OPERATION 


B1 


Introduction 


1.  The  Tailwater  Quality  Model  (TWQM)  is  a  user-friendly,  interactive, 
menu-driven,  PC-based  model  of  water  quality  in  reservoir  tailwaters.  TWQM 
computes  the  one -dimensional  (ID),  longitudinal,  steady-state  distribution  of 
water  quality  along  the  stream  axis,  beginning  at  the  reservoir  release  and 
progressing  downstream  for  any  number  of  miles  as  specified  by  the  user.  The 
model  was  developed  to  evaluate  downstream  water  quality  characteristics 
resulting  from  reservoir  releases.  Examples  of  the  types  of  issues  affecting 
tailwater  quality  that  might  be  addressed  include: 

a.  Changes  in  the  reservoir  outlet  elevation  and/or  release  flow. 

b.  Changes  in  release  reaeration  characteristics,  such  as  adding  hydro- 
power  (i.e.,  add-on  hydropower). 

c.  Changes  in  release  water  quality  due  to  degradation  of  in-pool  water 
quality. 

d.  Changes  in  the  downstream  channel,  such  as  reregulation. 

TWQM  focuses  on  the  problematic  water  quality  variables,  such  as  low  dissolved 
oxygen  and  reduced  substances  (e.g. ,  dissolved  iron,  dissolved  manganese, 
ammonium,  and  sulfide),  but  includes  biochemical  oxygen  demand,  major  nutri¬ 
ents,  and  algae.  The  model  can  also  be  used  to  evaluate  the  downstream  dis¬ 
tance  required  for  recovery  of  water  quality  to  some  standard.  TWQM  is 
designed  for  ease  of  application  with  limited  input  requirements  (due  mainly 
to  the  steady-state  assumption).  A  first-time  user  should  be  able  to  apply 
the  model  within  a  few  days .  An  experienced  user  should  be  able  to  apply  the 
model  within  a  day,  given  that  all  the  data  are  in  hand. 

2.  TWQM  is  composed  of  a  number  of  programs  (i.e.,  pre-  and  post¬ 
processor  programs  and  a  numerical  model)  that  are  menu  driven.  TAILWTR.FOR 
is  the  numerical  one -dimensional  tailwater  model  based  on  the  US  Environmental 
Protection  Agency  riverine  water  quality  model  QUAL2E.  TAILWTR.FOR  predicts 
estimates  of  constituent  concentrations  in  a  river  reach  based  on  the  input 
parameters  given.  The  main  batch  file  that  executes  all  the  programs  that 
make  up  TWQM  is  TWQB.BAT.  Descriptions  of  each  program  and  its  output  are 
given  in  Table  B1 .  The  order  in  which  TWQB.BAT  calls  these  files  is  shown  in 
Figure  Bl.  The  first  several  programs  executed  by  TWQB  are  pre-processor 
programs  that  assist  the  user  in  using  TWQM  (giving  general  information  and 
instructions)  and  also  assist  the  user  in  creating  an  input  data  set  for 
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Table  B1 


File 

TWQB.BAT 

TEXPLAIN.FOR 

TMENU.FOR 

TWQI .  FOR 

TFIX.FOR 

TAILWTR.FOR 

TREFORM.FOR 

TPREDICT . DAT 

TPFILE.FOR 
and  TPFILEPL.FOR 

TWGRAPH.BAT 
and  TWPLOT.BAT 

TPGRAPH.BAT 
and  TPPLOT.BAT 


TWQM  File  Description 


_ Description _ 

Batch  program  that  runs  all  files  contained  in  the  Tailwater 
Quality  Model. 

Prints  screens  containing  general  information  and 
instructions  for  TWQM. 

Allows  the  user  to  select  which  water  quality  constituents 
are  to  be  modeled  and  input  headwater  boundary  conditions 
for  those  constituents  selected. 

Creates  the  input  data  file  required  by  TAILWTR.FOR.  Also 
creates  the  file  containing  observed  data,  OBSERVED.DAT, 
which  is  used  by  K4GRAPH.BAT  to  graph  observed  data. 

Allows  the  user  to  alter  an  existing  input  data  file. 

Determines  the  tailwater  quality  of  dam  discharges .  Creates 
a  user-named  output  file,  and  the  files  PREDICT.DAT  and 
PREDICT. HYD. 

Reads  PREDICT.DAT  and  PREDICT. HYD  and  combines  them  into  one 
file,  TPREDICT.DAT. 

Data  file  for  all  graphing  (graph  viewed  on  monitor)  and 
plotting  (hard  copy  of  graph)  files.  Contains  output  for 
all  elements  for  selected  constituents  from  the  last  run  of 
TAILWTR.FOR. 

Creates  the  files  used  by  TPGRAPH.BAT  (*.TPL)  and 
TPPLOT.BAT  (*.PLT)  for  graphing  and  plotting  TAILWTR.FOR 
output . 

Batch  files  which  call  the  graphing  and  plotting  programs 
written  for  GRAPHER  software  (Golden  Software,  Inc.  1988). 
Uses  files  with  AXS,  GRF,  and  PLT  extensions. 

Batch  files  which  call  graphing  and  plotting  software 
written  for  TPLOT  software  (Benton  1987) .  Uses  files  with 
TPL  and  PLT  extensions. 
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TWQB  STRUCTURE 


FILES  CREATED 

T1TLE.DAT 

TAILWTR  DATA  FILE* 
OBSERVED.DAT 

TAILWTR  OUTPUT  FILE* 
PREDICT.DAT 
PREDICT. HYD 

TPREDICT.DAT 

*.TPL 

*.PLT 


These  files  can  have  any  name  up  to  eight  characters  in  length 
and  an  extension  up  to  three  characters  long. 

Figure  Bl.  Order  in  which  TWQB.BAT  calls  up  files 


TAILWTR.FOR.  Output  from  TAILVTR.FOR  is  plotted  using  commercial  software 
such  as  GEIAPHER  (Golden  Software,  Inc.  1988)*  or  TPLOT  (Benton  1987). 


Tailwater  Model  Input  Screens 

3.  The  TWQM  is  a  user-friendly,  interactive  model.  Screens  have  been 
developed  to  assist  the  user  with  model  operations  and  to  prompt  the  user  for 
appropriate  information.  The  initial  screens  contain  cautionary  comments 
about  the  model  and  input  requirements.  If  the  user  chooses  to  create  a  new 
input  data  set  for  TAILWTR.FOR,  he  is  prompted  for  the  data  required.  In  many 
instances  these  screens  already  contain  default  values  for  the  data.  To 
select  ALL  of  the  default  values  shown  on  a  screen,  answer  no  (N)  when  asked 
if  changes  need  to  be  made.  To  select  some  of  the  default  values  shown  on  the 
screen  and  enter  new  values  for  the  others,  enter  yes  (Y) .  Enter  the  new 
values  where  desired  and  press  ENTER  to  go  to  the  next  value.  To  select  the 
default  value  shown,  press  ENTER.  The  TWQM  is  very  specific  on  how  the  data 
are  input  with  respect  to  integer  and  real  numbers.  Unless  prompted  other¬ 
wise,  all  numerical  input  should  be  real.  A  note  to  this  effect  will  appear 
on  all  screens.  In  Instances  where  integer  input  is  required,  a  note  will 
appear  indicating  so. 

4.  The  user  is  required  to  input  information  in  TMENU,  TWQI ,  TFIX, 
TAILWTR,  and  the  graphics  program  GRAFMENU.BAT  (optional).  In  all  cases  the 
required  input  will  be  an  answer  to  a  direct  question  appearing  on  the  screen. 
Types  of  entries  required  by  the  user  are  real  or  integer  numbers,  yes  or  no 
answers,  data  file  names,  or  a  short  descriptive  text  (i.e.,  river  name). 

5.  In  the  following  section,  screens  are  shown  in  the  order  that  they 
appear  to  the  user  when  explanatory/cautionary  comments  are  required  for  clar¬ 
ification,  or  when  user  input  is  required.  A  short  description  of  each 
screen,  along  with  information  concerning  the  input,  is  also  included. 


*  See  References  at  the  conclusion  of  the  main  text. 
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Screens  in  TEXPLAIN.FOR 


6.  This  program  does  not  require  any  user  input,  but  the  information 
contained  on  its  two  screens  is  VERY  important.  The  first  screen  contains 
general  information  on  program  TWQB.BAT  (Figure  B2) .  It  also  contains 
instructions  for  exiting  TWQM  in  case  of  an  error  by  pressing  the  control 
(CTRL)  key  and  'C'  simultaneously.  It  must  be  noted  that  this  should  be  done 
only  when  the  user  desires  to  start  all  over  again.  TWQM  is  structured  so 
that  a  user  is  given  opportunities  to  change  answers  if  incorrect  values  are 
entered.  'Control  C'  should  be  used  only  as  a  last  resort  since  the  user  will 
lose  any  information  that  a  particular  program  in  TWQM  has  not  written  to  a 
file. 


•  O 

•  **********  TWQB.BAT  **********  • 

•  o 

•  TWQB  is  the  main  batch  progam  that  submits  programs  to:  * 

•  o 

•  1)  create  a  new  input  data  set  for  TWQM  or  alter  an  * 

°  old  data  set,  • 

•  2)  run  TWQM, 

•  3)  and  plot  simulation  results.  * 

•  o 

•  If  a  mistake  is  made  at  anytime  during  the  batch  run,  * 

•  "control  C"  out  of  TWQB. 


PRESS  "ENTER"  TO  CONTINUE 


iAAAAAAAAAAAAdAAAAAAAddAAAAAAAAAdAdAAdAAAdAAAAAAAAAAAdAAAAAAAAAAAAAAAAAAAAAAAAAi 

Figure  B2.  TEXPLAIN  first  screen 
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7.  The  next  screen  in  TEXPLAIN  contains  more  cautionary  comments  (Fig¬ 
ure  B3) .  The  most  important  one  is  that  release  values  for  headwater  boundary 
conditions  (see  Headwater  Data  Screen)  are  required  before  running  TWQH.  If 
these  values  are  not  available,  they  can  be  estimated  from  in-pool  profiles 
using  the  selective  withdrawal  model,  SELECT  (Davis  et  al.  1987).  If  SELECT 
must  be  run,  TWQM  should  be  exited  at  this  time  by  pressing  the  CTRL  key  and 
'C'  simultaneously. 


o  Q 

*  *******  CAUTIONARY  COMMENTS  ******* 

*  o 

°  1)  Before  beginning  TWQB,  values  for  headwater  boundary  conditions  ' 

*  are  needed  whether  they  are  observed  historical  or  estimated  data  * 

*  obtained  from  running  SELECT.  If  you  begin  TWQB  before  you  have  * 

*  estimated  boundary  conditions,  "control  C"  now  and  • 

*  run  SELECT  to  get  estimated  values.  ° 

o  • 

*  2)  Once  TWQB  has  been  submitted,  it  is  important  to  observe  whether  “ 

*  your  response  to  questions  requires  a  real  or  integer  number.  * 

*  You  will  get  a  chance  to  correct  mistakes  so  pay  attention  * 

*  in  case  of  incorrect  responses!!!!  • 

*  • 

•  3)  The  input  data  set  contains  parameters  that  are  not  required  • 

•  by  TWQM  (i.e.  ENDATA9  parameter)  but  are  particular  to  QUAL2E.  • 

•  Since  QUAL2E  was  the  basis  of  TWQM  development,  these  parameters  * 

•  are  still  read  by  the  program.  Future  improvements  to  TWQM  will  ’ 

•  include  removing  unnecessairy  parameters.  * 

•  • 

’  PRESS  "ENTER"  TO  CONTINUE  TO  MENU. FOR!!!! 


Figure  B3 .  TEXPLAIN  second  screen 
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Screens  in  TMENU.FOR 

8.  The  first  screen  in  TMENU  contains  general  information  concerning 
the  program  (Figure  B4) .  Note  that  if  incorrect  responses  are  entered  on  the 
following  table  (Figure  B5) ,  there  will  be  an  opportunity  to  correct  them 
before  proceeding  to  the  next  screen. 


o  o 

*  **********  TMENU.FOR  **********  • 

*  o 

*  TMENU.FOR  is  the  first  program  in  a  series  of  programs  • 

*  that  create  the  input  data  set  for  TWQM.  TMENU.FOR  • 

*  allows  the  user  to  select  the  water  quality  constituents  * 

’  they  wish  to  model.  It  also  sets  up  the  headwater  * 

*  boundary  input  of  the  constituents  selected.  * 

®  e 

*  The  following  table  of  state  variables  contains  water  • 

*  quality  constituents  that  can  be  modeled  by  TWQM.  The  * 

*  constituents  with  Y  preceding  them  will  be  modeled  by  TWQM,  * 

*  If  you  wish  to  change  the  constituents  being  modeled,  answer  * 

*  Y  or  N.  If  you  make  a  mistake,  you  CAN  go  back  and  correct  * 

*  the  answer.  Press  "Enter"  to  continue  to  next  constituent  * 

*  or  if  the  answer  does  not  need  to  be  changed.  * 


•  PRESS  "ENTER"  TO  CONTINUE  ! ! ! ! 

•  • 

•  • 

dAddddddddddddddddddAdddddddddddddAdddddddAdAddddAdAdddAddAdAAdddAdddAddAdAddddi 

Figure  B4.  TMENU  first  screen 
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9.  The  next  screen  in  TMENU  lists  all  constituents  that  TAILWTR  can 
model  (Figure  B5) .  The  constituents  that  will  be  modeled  are  selected  from 
this  list  by  entering  'Y'  at  the  flashing  cursor.  Choices  made  here  are  very 
important  since  they  determine  not  only  which  processes  will  be  simulated  but 
also  what  input  data  are  required. 


•  ANSWER  CONSTITUENT 


N  CONSERVATIVE  1  (mg/1) 

N  CONSERVATIVE  2  (mg/1) 

N  TEMPERATURE  (C) 

Y  BOD  (mg/1) 

N  CHL  ALGAE  (rag/1) 

N  PHOSPHOROUS  CYCLE  AS  P  IN  mg/1 

(ORGANIC- P;  DISSOLVED- P) 

Y  NITROGEN  CYCLE  AS  N  IN  mg/1 

(ORGANIC-N;  AMMONIA-N;  NITRATE-N;  NITRITE-N  ) 

Y  DISSOLVED  OXYGEN  (mg/1) 

N  ANC  (mg/1) 

Y  IRON  (mg/1) 

Y  MANGANESE  (mg/1) 

Y  SULFIDE  (mg/1) 

N  SULFATE  (mg/1) 

ARE  THERE  ANYMORE  CHANGES  TO  BE  MADE  (Y/N)? 


Figure  B5.  TMENU  second  screen 
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10.  Using  the  title  information  screen,  the  user  enters  two  titles  for 
this  simulation  (Figure  B6) .  For  example,  if  a  user  decides  to  decrease  the 
release  flow  from  a  project  by  half,  one  of  the  titles  could  be  "Nimrod 
Reservoir  -  Flow  Halved."  The  other  title  could  be  used  to  identify  project 
name  and  date  of  simulation.  These  titles  can  be  up  to  80  characters  long  and 
are  the  first  two  lines  in  the  input  data  file  for  TAILWTR. 

O  • 

•  • 

•  icic'k’kie'kicicicieieiciciciHcicifie-k-k'k-kit'k'k'k'kic-k-kieic'k-k'k-k-k-k-k-k'k'k'k'kic-kie-k  • 

•  *  *  • 

•  *  TITLE  INFORMATION  * 

•  *  *  • 

•  'k'k'k'k-k-k'k'k'k'kicic-k'k-k'kie-kic-k'k-kick'k-k-kic-k'k'k-k-k'k'k'k'k-k'k'k-k'k-k’k'k'k-k-k'k  • 

•  e 

'  ENTER  TITLE  1  (i.e.  STREAM  QUALITY  MODEL  WITH  FLOWS  DOUBLED) 

“  NIMROD  89 

•  e 

•  ENTER  TITLE  2  (i.e.  DATA  SET  ONE  OR  SCENARIO  ONE) 

•  STEADY  STATE  (200  CFS) 


AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA^ 

Figure  B6.  Title  information  screen 
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11.  The  next  screen  in  TMENU  (Figure  B7)  consists  of  a  list  of  the 
constituents  that  were  selected  earlier  (Figure  B5)  to  be  modeled.  Answer  yes 
(Y)  to  those  constituents  for  which  the  headwater  boundary  conditions  are 
known  and  no  (N)  to  those  for  which  boundary  conditions  are  not  knovm.  On 
subsequent  screens  the  observed  boundary  condition  values  and  estimated  (using 
SELECT)  boundary  conditions  are  entered  for  all  constituents  being  modeled 
(Figure  B8) . 


•  ANSWER  CONSTITUENTS  TO  BE  MODELED 


Y  BOD  (mg/1) 

Y  NITROGEN  CYCLE  AS  K  IN  mg/1 

Y  DISSOLVED  OXYGEN  (mg/1) 

Y  IRON  (mg/1) 

Y  MANGANESE  (mg/1) 

Y  SULFIDE  (mg/1) 

pH  (NOT  MODELED  AT  THIS  TIME,  BUT  INPUT  REQD) 


*  These  are  the  constituents  you  chose  to  model.  * 

*  Answer  Y  for  those  constituents  that  you  have  • 

*  boundary  conditions  for,  and  N  to  those  you  • 

*  do  not.  For  those  you  answered  N  to,  boundary  * 

°  conditions  can  be  estimated  from  SELECT  or  * 

*  can  be  estimated  from  historical  data.  * 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAi. 

Figure  B7 .  Constituents  to  be  modeled 
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•  CONSTITUENTS  WITH  -KNOWN-  BOUNDARY  CONDITIONS 

*  Enter  observed  boundary  condition  value  for  each  constituent.  * 


CONSTITUENT 

DEFAULT 

BOD  (mg/1) 

1.00 

ORGANIC  N  (mg/1) 

0.32 

AMMONIA  N  (mg/1) 

0.28 

NITRITE  N  (mg/1) 

0.00 

NITRATE  N  (mg/1) 

0.03 

DISSOLVED  OXYGEN  (mg/1) 

7.00 

IRON  (mg/1) 

1.00 

MANGANESE  (mg/1) 

1.20 

SULFIDE  (mg/1) 

.00016 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAl 

Figure  B8 .  Known  boundary  conditions 
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12.  TWQI  is  the  main  interactive  program  in  TWQM  and  serves  as  a 
preprocessor  that  generates  the  input  data  file  for  TAILV/TR.  Since  TAILWTR 
was  developed  using  QUAL2E,  the  data  file  generated  by  TWQI  is  very  similar  to 
a  QUAL2E  input  data  file. 

13.  TWQI  prompts  the  user  for  most  input  required  to  generate  a  data 
file  for  TAILWTR.  Some  data  such  as  evaporation  coefficients  and  certain 
constants  are  included  in  TWQI  but  are  not  user  definable  and  will  not  appear 
on  any  menus .  There  is  usually  a  default  value  for  the  data  that  the  user 
must  enter.  The  notable  exceptions  are  input  dealing  with  flow  and  location 
(river  mile).  A  value  MOST  be  entered  for  these  before  any  simulation  can  be 
performed  even  though  a  default  value  of  0.0  appears  on  the  screen. 

14.  TWQI  is  very  specific  concerning  input  format.  The  importance  of 
reading  and  following  the  directions  at  the  bottom  of  each  screen  cannot  be 
overemphasized.  Most  values  that  are  input  into  TWQI  must  be  real  numbers. 
However,  there  are  certain  instances  when  integer  values  are  required  by  TWQI. 
Unless  indicated  otherwise,  all  numeric  input  is  to  be  in  real  numbers. 

15.  The  first  screen  of  TWQI  is  the  title  screen  which  contains  a 
description  of  TWQI  (Figure  B9) .  Upon  pressing  RETURN  (or  ENTER) ,  TWQI 
prompts  the  user  for  the  name  of  the  input  data  file  TWQI  will  generate.  This 
file  will  be  a  new  file;  the  name  of  an  existing  file  must  not  be  used.  If 
the  file  name  of  an  existing  file  is  used,  TWQI  will  stop.  This  precaution 
prevents  a  user  from  erroneously  writing  over  an  existing  data  file.  The  file 
name  must  not  exceed  eight  characters  and  may  have  an  extension  of  up  to  three 
characters . 


•  • 

•  ******************************************************  • 


* 

* 

* 

TTTTTTT 

WWW  WW  WWW 

QQQQQQQ 

IIIIIIII 

* 

* 

TTT 

WWW  WW  WWW 

QQ  QQ 

nil 

* 

* 

TTT 

WWW  WW  WWW 

QQ  QQ 

nil 

* 

* 

TTT 

WWW  WW  WWW 

QQ  QQQ 

nil 

* 

* 

* 

TTT 

wwwwwwww 

QQQQQQQQQ 

nnnn 

* 

* 

* 

* 

////////////////////////////////////////////// 

* 

* 

•  ****************************************************  • 

•  • 

•  TWQI  IS  AN  INTERACTIVE  PROGRAM  DESIGNED  TO  CREATE 

’  THE  INPUT  DATA  SET  OF  THE  TAILWATER  QUALITY  MODEL, 

•  TWQM,  DEVELOPED  BY  THE  US  ARMY  ENGINEER  WATERWAYS 

•  EXPERIMENT  STATION,  VICKSBURG,  MS 

•  o 

•  PRESS  RETURN  TO  CONTINUE 

•  o 

•  ENTER  NAME  OF  DATA  SET  YOU  WISH  TO  CREATE 

•  NIMROD89.DAT 

•  « 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAi 

Figure  B9 .  Title  screen 
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16.  The  next  screen  echoes  the  data  input  earlier  in  TMENU  (Figure  B5) 
dealing  with  what  constituents  will  be  modeled  (Figure  BIO) .  This  information 
should  not  be  changed  unless  it  was  entered  incorrectly  in  TMENU. 


TITLEOl 

NIMROD  89 

TITLE02 

STEADY  STATE  (200  CFS) 

TITLE03 

NO 

CONSERVATIVE  1  (mg/1) 

TITLE04 

NO 

CONSERVATIVE  2  (mg/1) 

TITLE05 

NO 

TEMPERATURE  (C) 

TITLE06 

YES 

BOD  (mg/1) 

TITLE07 

NO 

CHL  ALGAE  (mg/1) 

TITLE08 

NO 

PHOSPHOROUS  CYCLE  AS  P  IN  mg/1 

TITLE09 

(ORGANIC-P;  DISSOLVED-P) 

TITLEIO 

YES 

NITROGEN  CYCLE  AS  N  IN  mg/1 

TITLEll 

(ORGANIC-N;  AMMONIA-N;  NITRATE-N; 

NITRITE -N  • 

TITLE12 

YES 

DISSOLVED  OXYGEN  (mg/1) 

TTTLE13 

pH  (NOT  MODELED  AT  THIS  TIME,  BUT 

INPUT  REQD* 

TirLE14 

NO 

ANC  (mg/1) 

TITLE15 

YES 

IRON  (mg/1) 

TITLE16 

YES 

MANGANESE  (mg/1) 

TITLE17 

YES 

SULFIDE  (mg/1) 

TITLE18 

NO 

SULFATE  (mg/1) 

THE  CONSTITUENTS 

WITH  "YES"  PRECEDING  THEM  ARE  TO  BE  MODELED 

BY  TWQM, 

DO  YOU 

WISH  TO  MAKE  CHANGES  (Y/N)? 

• 

Figure  BIO.  Currently  modeled  constituents 
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17.  Screen  2  contains  the  program  control  data  (Figure  fill).  The 
screen  is  self-explanatory.  If  no  (N)  is  entered  for  the  input  metric  ques¬ 
tion,  the  input  must  be  in  the  non-metric  units  specified.  If  no  (N)  is 
entered  for  the  output  metric  option,  the  output  will  be  in  non-metric  units. 
One  item  of  interest  is  the  value  for  the  average  slope  of  the  riverbed. 

These  data  are  required  by  the  algorithm  that  determines  the  computational 
element  length.  A  default  value  corresponding  to  a  slope  of  1  ft  to  10,000  ft 
is  shown.  This  value  can  be  changed,  but  a  slope  of  0  ft/ft  is  not  allowed. 


^eeeeeeeeeeeeegee§eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee£ 

n  n 

n  *****************************************************  n 


n 

n 

R 

a 

a 

a 


* 

* 

* 

* 


SCREEN  TWO 

DATA  TYPE  1:  PROGRAM  CONTROL  DATA 


* 

* 

* 

* 


***************************************************** 


R 

R 

H 

R 

R 

R 


R 

INPUT  METRIC  (Y/N) 

NO 

R 

OUTPUT  METRIC  (Y/N) 

NO 

R 

NUMBER  OF  REACHES 

3.00 

R 

NUMBER  OF  POINT  LOADS  AND  WITHDRAWAL 

0.00 

R 

LATITUDE  OF  BASIN  (DEG) 

34.57 

R 

LONGITUDE  OF  BASIN  (DEG) 

93.10 

R 

STANDARD  MERIDIAN  (DEG) 

90.00 

R 

DAY  OF  YEAR  START  TIME 

203.00 

R 

ELEV.  OF  BASIN  (ELEV) 

480.00 

R 

AVG.  SLOPE  OF  RIVER  BED 

0.001000 

R 

COBBLE  OR  SEDIMENT  BOTTOM  (C  or  S) 

S 

R 


R 

R 

R 

R 

R 

R 

n 

R 

R 

R 

R 

R 


R  H 
R  ARE  THERE  ANY  CHANGES  TO  BE  MADE  (Y/N)?  n 
R  R 


4eeeeegee€eeee§6Si§§§ee6Sgeeee§geeeeeeeeeeeeeeeeeeeeeeeeeee§eeeeeg§eeeeeeeeeeee¥ 


Figure  Bll.  Program  control  data  screen 
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18.  The  next  screen  (Figure  B12)  asks  the  user  if  changes  to  the  theta 
values  are  required.  Theta  values  are  used  to  correct  for  the  effects  of 
temperature  upon  reaction  rates.  If  no  (N)  is  entered,  TWQI  proceeds  to 
screen  four,  and  TAILWTR  uses  default  theta  values  (Brown  and  Barnwell  1987). 
If  the  answer  is  yes  (Y) ,  a  list  of  the  default  theta  values  appears,  and 
changes  can  be  made  to  any  value  (Figure  B13) .  Once  all  changes  are  made,  the 
theta  values  are  printed  to  the  screen  one  last  time. 

o  o 

•  *^MHIr*^hk**********************************************  * 


*  * 

*  SCREEN  THREE  * 

*  DATA  TYPE  IB:  THETA  VALUES  * 

*  * 


★★^Hr************************************************* 

DO  YOU  WISH  TO  CHANGE  THE  TEMPERATURE  CORRECTION 
CONSTANTS  (THETA  VALUES)  FOR  THE  RATE  COEFFICIENTS  (Y/N)? 


AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAi 

Figure  B12.  First  theta  screen 
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RATE  CODE 

THETA  VALUE 

RATE  CODE 

THETA  VALUE 

CBOD  DEC 

1.0470 

DISP  SRC 

1.0740 

BOD  SETT 

1.0240 

ALG  GROW 

1.0470 

OXY  TRAN 

1.0240 

ALG  RESP 

1.0470 

SOD  RATE 

1.0600 

ALG  SETT 

1.0240 

ORGN  DEC 

1.0470 

SULF  OXR 

1.0500 

ORGN  SET 

1.0240 

ANC  DEC 

1.0000 

NH3  DECA 

1.0830 

ANC  SETT 

1.0240 

NH3  SRCE 

1.0740 

ANC  SRCE 

1.0000 

N02  DECA 

1.0470 

IRON  OXR 

1.1400 

FORG  DEC 

1.0470 

SULF  RED 

1.0000 

PORG  SET 

1.0240 

MANl  OXR  T<25 

1.0540 

MAN2  OXR  T>25 

1.1050 

ENTER  NEW  VALUES  (REAL  NUMBERS)  AT  "FLASHING  CURSOR" 
PRESS  "ENTER"  TO  CONTINUE 


AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAi. 

Figure  B13.  Theta  changes  screen 
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19.  Screen  4  (Figure  B14)  contains  the  headwater  boundary  conditions 
entered  in  TMENU.  This  screen  appears  to  the  user  in  two  stages.  First,  the 
screen  title  block  and  a  prompt  for  entering  project  release  flow  appear  on 
the  screen.  Once  a  value  for  flow  has  been  entered,  the  rest  of  the  headwater 
boundary  conditions  appear  on  the  screen.  At  that  point,  any  or  all  of  the 
data  appearing  may  be  changed. 


•  • 

•  ****************************lfr*******************************  • 


* 

SCREEN 

FOUR 

* 

* 

HEADWATER  DATA 

* 

************************************************************ 

CONSTITUENT 

VALUE 

CONSTITUENT 

VALUE 

FLOW  (cfs) 

200.000 

ORGANIC  N  (mg/1) 

0.32000 

TEMPERATURE  (F) 

78.60000 

AMMONIA  N  (mg/1) 

0.28000 

DISSOLVED  OXYGEN  (mg/1) 

7.00000 

NITRITE  N  (mg/1) 

0.00000 

BOD  (mg/1) 

1.00000 

NITRATE  N  (mg/1) 

0.03000 

IRON  (mg/1) 

1.00000 

ORGANIC  P  (mg/1) 

0.00000 

MANGANESE  (mg/1) 

1.20000 

DISSOLVED  P  (mg/1) 

0.00000 

SULFIDE  (mg/1) 

0.00016 

ANC  (mg/1) 

0.00000 

SULFATE  (mg/1) 

0.00000 

CONSERVATIVE  1  (mg/1) 

0.00000 

pH 

6.50000 

CONSERVATIVE  2  (mg/1) 

0.00000 

CHL  ALGAE  (mg/1) 

0.00000 

*  DO  YOU  NEED  TO  MAKE  ANY  CHANGES  (Y/N)? 

•  • 

Figure  B14.  Headwater  boundary  conditions  screen 
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20.  Following  the  headwater  data  screen  is  the  data  type  2,  reach 
identification  screen  (Figure  B15).  Each  reach  may  be  given  its  own  name 
(maximum  of  15  characters),  or  the  same  name  may  be  used  for  all.  The  value 
of  the  river  kilometer  or  mile  at  the  head  of  the  reach  must  be  greater  than 
that  at  the  end  of  the  reach.  This  screen  is  repeated  once  for  each  reach. 
After  the  first  reach,  the  river  kilometer/mile  at  the  head  of  the  next  reach 
is  automatically  set  to  the  river  kilometer/mile  at  the  end  of  the  previous 
reach  and  cannot  be  changed  by  the  user. 


•  • 
•  *****************************************************  ' 


*  * 

*  SCREEN  FIVE  * 

*  DATA  TYPE  2:  REACH  IDENTIFICATION  DATA  * 

*  REACH  1  * 

***************************************************** 

RIVER  NAME  FOR  SIMULATION:  NIMROD  1 

RIVER  MILE  AT  HEAD  OF  REACH  62.600 

RIVER  MILE  AT  END  OF  REACH  60.800 


ARE  THERE  ANY  CHANGES  (Y/N)?  Y 

ENTER  ALL  VALUES  AS  REAL  NUMBERS 
PRESS  "ENTER"  TO  CONTINUE 


Figure  B15.  Reach  identification  data 
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21.  After  screen  5,  TWQI  prompts  the  user  for  the  number  of  point  loads 
or  withdrawals  within  that  reach  (Figure  B16) .  This  number  cannot  exceed  the 
value  entered  earlier  on  the  control  data  screen  (Figure  Bll)  for  the  total 
number  of  point  loads  and  withdrawals.  If  it  does,  a  warning  flag  informs  the 
user  that  the  number  of  point  loads  and  withdrawals  within  this  reach  or  the 
total  number  within  all  reaches  has  been  exceeded  and  must  be  corrected. 

22.  Next,  a  prompt  appears  asking  what  is  first  from  the  upstream  end 
of  the  reach,  a  point  load  (P)  or  withdrawal  (W) .  After  entering  P  or  W, 
prompts  for  the  river  kilometer/mile  location  of  the  point  load  or  withdrawal 
and  its  identification  name  appear.  These  three  steps  are  repeated  for  all 
point  loads  and  withdrawals  within  a  reach. 

OddddddddddddddAdAdAdAAAdddddddddddddddAAdAddAdAddddddaddddddddiddddddddddddiddC 

o  • 

•  • 

•  REACH  1 

o  • 

“  ARE  THE  FOLLOWING  OR  OTHER  TYPES  OF  POINTLOAD  SOURCES  PRESENT 

•  IN  THIS  REACH  (Y/N)?  N 

P  • 

•  A:  POINT  LOADS  (eg.  WASTEWATER  TREATMENT  PLANT  OUTFLOWS) 

'  B:  WITHDRAWALS 


Figure  B16.  Point  load/withdrawal  locations 
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23.  Screen  6  (Figure  B17)  contains  the  hydraulic  data  (such  as  depth, 
velocity,  and  Manning's  n)  for  this  reach.  Values  must  be  entered  for  all 
three  hydraulic  variables.  If  multiple  reaches  are  simulated,  the  values  for 
the  previous  reach  are  displayed  as  the  defaults  of  the  current  reach.  These 
values  can  be  selected  by  simply  pressing  ENTER,  or  new  values  may  be  entered 


•  o 

•  **■***■**********★*********■***★***************■****'*****  • 


• 

• 

• 

• 

o 

• 

*  * 

*  SCREEN  SIX  * 

*  DATA  TYPE  5:  HYDRAULIC  DATA  * 

*  REACH  1  * 

***************************************************** 

o 

• 

• 

o 

o 

o 

• 

• 

HYDRAULIC  DATA  TYPE 

DEFAULT 

o 

• 

o 

• 

• 

HYDRAULIC  DEPTH  IN  FEET: 

2.000 

• 

• 

• 

• 

• 

AVERAGE  VELOCITY  IN  ft/s: 

.5600 

o 

• 

• 

• 

MANNINGS  n: 

0.100 

• 

• 

• 

• 

• 

ARE  THERE  ANY  CHANGES  (Y/N)?  Y 

o 

• 

• 

• 

• 

ENTER  ALL  VALUES  AS  REAL  NUMBERS 

• 

• 

«  • 

Figure  B17.  Hydraulic  data 
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24.  Coefficients  for  DO/BOD  used  by  TAILWTR  are  listed  on  screen 
7  (Figure  B18).  Also  listed  are  choices  for  the  method  of  calculating  the 
reaeration  coefficient  K2  .  TWQM  provides  the  same  eight  options  as  QUAL2E 
(Brown  and  Barnwell  1987),  allowing  the  user  to  select  estimation  or  reading 
in  of  K2  values.  These  options  are  as  follows: 


•  • 

•  *****************************************************  • 
•  *  *  • 

•  *  SCREEN  SEVEN  * 

•  *  DATA  TYPE  6:  DO/BOD  COEFFICIENTS  * 

^  *  REACH  1  * 


REACTION  COEFFICIENT  TYPE 

CBOD  DECAY  RATE  (1/day) 
CBOD  SETTLING  RATE  (1/day) 
SOD  RATE  (g  02/sqft-day) 

K2  OPTION 


DEFAULT 

0.150 

0.000 

0.200 

3.000 


ANY  CHANGES  TO  REACTION  COEFFICIENTS  (Y/N)? 


K2  OPTION 
1 
2 

3 

4 


K2  COMPUTATIONAL  METHOD  K2  OPTION 
READ  IN  VALUE  5 

CHURCHILL  6 

OCONNER  AND  DOBBINS  7 

OWENS  AND  GIBBS  8 


K2  COMPUTATIONAL  METHOD 
THACKSTON  AND  KRENKEL 
LANGBIEN  AND  DURUM 
K2  -  aQ**b 
TSIVOGLOU-WALLACE 


Figure  B18.  DO/BOD  coefficients 


Kr>  Option  1 


25.  Option  1  allows  the  user  to  read  in  estimated  or  measured  K2 
values . 

Ko  Option  2 

26.  Churchill,  Elmore,  and  Buckingham  (1962)  developed  the  following 
equation  for  K2  at  20  °C  based  on  measured  stream  reaeration  data. 


K,  =  5.026  u°  ®®®  d'^  ®”  2.31 


where 


reaeration  coefficient,  1/day 
average  stream  velocity,  ft/sec 
average  stream  depth,  ft 
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K-,  Option  3 

27.  O'Connor  and  Dobbins  (1958)  proposed  equations  based  on  the  turbu 
lence  characteristics  of  a  stream.  For  streams  that  have  low  velocities  and 
isotropic  conditions,  the  equation  is  given  as 


K, 


(Di„u)0-5 

dl-5~ 


(B2) 


28.  For  streams  with  higher  velocities  and  nonisotropic  conditions,  the 
equation  is  given  as 


K, 


480  “ 

"  - TTTS — 


2.31 


(B3) 


where 

D„  -  molecular  diffusion  coefficient,  ftV^ay  and  is  found  as 
D„  -  1.91  X  10^  (1.037)f-20 

So  -  slope  of  the  streambed,  ft/ft 

d  “  mean  stream  depth,  ft 

Equation  B2  has  been  found  to  be  generally  applicable  to  most  streams  and  is 
the  default  K2  option  for  TAILWTR. 

Ko  Option  4 

29.  Owens,  Edwards,  and  Gibbs  (1964)  estimated  reaeration  coefficients 
for  shallow,  fast-moving  streams  based  on  six  streams  in  England.  Combining 
their  work  with  that  of  Churchill,  Elmore,  and  Buckingham  (1962),  they 
developed  an  equation  for  streams  having  depths  of  0.4  to  11.0  ft  and  veloci¬ 
ties  of  0.1  to  5.0  ft/sec.  Their  equation  is  written  as 


K, 


9.4  u° 

dl” 


2.31 


(B4) 


where 

u  -  mean  velocity,  ft/sec 
d  -  mean  depth,  ft 
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K,  Option  5 

30.  Thackston  and  Krenkel  (1966)  proposed  the  following  equation  based 
on  their  investigation  of  several  rivers  in  the  Tennessee  Valley  Authority 
system. 


K,  =  10.8  (1  +  FO-5)  ^2.31 

^10  (J 


(B5) 


where  F  is  the  Froude  number .  which  is  given  by 


(B6) 


and  u*  is  the  shear  velocity  (ft/sec)  given  by 


u.  =  v'dS.g 


un  \/g 
1.49d^i®^ 


(B7) 


where 

d  -  mean  depth,  ft 
S,  -  slope  of  the  energy  gradient 
g  -  acceleration  of  gravity,  ft/sec^ 
u  -  mean  velocity,  ft/sec 
n  -  Manning's  coefficient 
Kn  Option  6 

31.  Langbien  and  Durum  (1967)  developed  a  formula  for  K2  at  20  'C: 

K,  =3.3  2.31  (B8) 

^1.33 
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where 


u  -  mean  velocity,  ft/sec 
d  -  mean  depth,  ft 
Ki  Option  7 

32.  Option  7  computes  the  reaeration  coefficient  from  a  power  function 
of  flow  and  is  given  by 


Kj  =  aQ** 


(B9) 


where 

a  -  coefficient  of  flow  for  K2 
Q  -  flow,  ft^/sec 
b  -  exponent  on  flow  for  K2 
Ko  Option  8 

33.  Tsivoglou  and  Wallace  (1972)  based  their  method  on  the  assumption 
that  the  reaeration  coefficient  for  a  reach  is  proportional  to  the  change  in 
elevation  of  the  water  surface  in  the  reach  and  inversely  prop'rtional  to  the 
flow  time  through  the  reach.  Their  equation  is  written  as 


K, 


=  c 


(BIO) 


where 

c  -  escape  coefficient,  1/ft 
Ah  -  change  in  water  surface  elevation  in  reach,  ft 
tf  -  flow  time  through  reach,  days 

34.  Assuming  uniform  flow,  the  change  in  water  surface  elevation  is 


Ah  =  Sg  Ax 


(BID 


where 

Sg  -  slope  of  energy  gradient,  ft/ft 
Ax  -  reach  length,  ft 
and  the  flow  time  through  the  reach  is 
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Ax 

u 


(B12) 


t 


where  u 
35. 


is  the  mean  velocity  in  the  reach  (ft/sec). 
Substituting  Equations  Bll  and  B12  into  BIO  gives 


=  (3,600  X  24)  cSgU 


(B13) 


36.  Equation  B13  is  the  formula  used  for  Option  8  in  TAILWTR.FOR.  The 
constants  3,600  and  24  convert  velocity  to  units  of  feet  per  day.  The  slope 
is  input  directly  for  computing  K2  with  this  option.  The  escape  coefficient 
can  be  estimated  based  on  recommendations  of  TenEch  (1978)  using 

c  -  0.054  ff^  (at  20  “C)  for  15  <  Q  <  3,000  ft^/sec 
c  -  0.110  ft*^  (at  20  “O  for  1  <  Q  <  15  ftVsec 

37.  Any  of  the  methods  listed  can  be  selected  by  entering  their  corre¬ 
sponding  number  (real  number).  If  K2  option  1,  7,  or  8  is  selected,  TWQI 
will  prompt  the  user  to  input  the  appropriate  coefficients  or  values. 
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38.  Screen  8  (Figure  B19)  contains  coefficients  for  nutrient  removal 
kinetics.  The  removal  rate  for  ammonia  oxidation  can  be  specified  or  deter¬ 
mined  by  TAILWTR  by  entering  a  '-1.'  in  place  of  an  actual  value.  The  word 
'modeled'  will  appear  on  the  screen  (Figure  B19)  for  ammonia  oxidation  to 
indicate  that  the  user  chose  for  TAILWTR  to  calculate  the  rate.  All  other 
rates  on  this  screen  are  set  to  default  values  if  no  changes  are  made  or 
specified  by  the  user.  Default  values  for  all  the  rates  are  average  values  of 
the  ranges  specified  by  the  QUAL2E  manual  (Brown  and  Barnwell  1987). 


•  o 

•  *-kiirkirk****************1rk***irk*irk********irk******it****  • 


*  * 

*  SCREEN  EIGHT  * 

*  DATA  TYPE  6A:  NUTRIENT  KINETICS/COEFFICIENTS  * 

*  REACH  1  * 


REACTION  COEFFICIENT  TYPE  DEFAULT 


ORG-N  HYDROLYSIS  RATE  (1/day)  0.400 

ORG-N  SETTLING  RATE  (1/day)  0.010 

NH-3  OXIDATION  RATE  (1/day)  MODELED 

BENTHAL  SOURCE  NH-3  (mg  02/sqft-day)  0.000 

NO- 2  OXIDATION  RATE  (1/day)  1.000 

ORG-P  DECAY  RATE  (1/day)  0.350 

ORG-P  SETTLING  RATE  (1/day)  0.050 

BENTHAL  SOURCE  DISP  (mg  P/sqft-day)  0.000 


*  ANY  CHANGES  TO  ANY  NITROGEN  COEFFICIENTS  (Y/N)? 

•  • 

Figure  B19.  Nutrient  kinetics 
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39.  Screen  9  (Figure  B20)  contains  values  for  coefficients  for  algae 
and  the  arbitrary  nonconservative  constituents.  Default  values  were  obtained 
from  ranges  recommended  in  the  QUAL2E  manual  (Brown  and  Barnwell  1987).  If 
any  changes  are  made  to  these  values ,  the  new  value  should  have  units  that  are 
consistent  with  those  shown. 


•  • 

•  **★***************'***********************************  • 

•  *  *  • 

•  *  SCREEN  NINE  * 

•  *  DATA  TYPE  6B:  ALGAE/OTHER  COEFFICIENTS  * 

•  *  REACH  1  * 

•  *****************************************************  • 

•  • 

*  REACTION  COEFFICIENT  TYPE  DEFAULT 


CHLA  to  ALGAE  RATIO  50.000 
ALGY  SETT  RATE  (ft/day)  3.000 
NON  ALGAL  LIGHT  EXT  COEF  (1/ft)  1.000 
ANC  DECAY  (1/day)  0.050 
ANC  SETT-COEFF  (1/day)  0.005 
BENTHAL  SOURCE  ANC  (mg  ANC/sqft-day)  0.000 


ANY  CHANGES  TO  COEFFICIENTS  (Y/N)? 


AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA^. 

Figure  B20.  Algae  coefficients 
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40.  Screen  9A  (Figure  B21)  contains  removal  rate  values  for  the  three 
reduced  substances  modeled:  sulfide,  iron,  and  manganese.  New  removal  rates 
may  be  entered  for  one  or  all  reduced  substances,  or  TAILWTR  can  determine 
their  removal  rates  (except  for  sulfide) .  New  removal  rates  are  entered  in 
the  same  manner  as  other  data  in  previous  screens.  To  have  TAILWTR  determine 
the  removal  rate  of  a  substance,  answer  yes  (Y)  when  asked  if  there  are  any 
changes.  As  with  ammonia  oxidation,  enter  a  -1.  for  the  removal  rate  of  the 
substance  whose  removal  rate  is  to  be  modeled.  The  word  'modeled'  will  appear 
on  the  screen  in  place  of  an  actual  value  for  iron/manganese  decay  whenever  a 
-1.  was  entered  as  the  removal  rate. 


o  # 

•  ****★*******************★*★★**★★★★*★*************★★'**  • 

•  *  *  • 

•  *  SCREEN  NINE  -  A  * 

•  *  DATA  TYPE  6B:  REDUCED  SUBSTANCES  REMOVAL  RATES  * 

•  *  REACH  1  * 

•  *****************************************************  • 

«  • 

•  REMOVAL  RATE  DEFAULT 


SULFIDE  DECAY  (1/DAY) 
IRON  DECAY  (1/DAY) 
MANGANESE  DECAY  (1/DAY) 


25.000 

MODELED 

MODELED 


ANY  CHANGES  TO  COEFFICIENTS  (Y/N)? 


Figure  B21.  Reduced  substances  removal  rates 
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41.  Screen  10  actually  appears  in  two  screens  and  only  for  simulations 
where  temperature  is  not  being  modeled.  This  screen  contains  information  on 
the  initial  conditions  in  a  reach.  First,  the  initial  pH  and  temperature  of 
the  reach  are  entered  (Figure  B22).  Initial  values  for  the  other  constituents 
are  set  equal  to  their  headwater  boundary  conditions  and  then  displayed  with 
the  pH  and  temperature  values  on  the  second  screen  10  (Figure  B23).  If 
desired,  changes  can  be  made  to  any  of  the  initial  condition  values  at  this 
time.  However,  changes  to  any  variables  besides  temperature  and  pH  will  not 
Impact  the  model  results  since  a  steady  state  solution  is  obtained. 


ic'kieic'kiricicic'k'k’kirk'k'it'k'k'icic'k'kicicic'kifif'kicicic'kiHcicirk'kicieic'irkic'k'k^ic'k^'k’k'k'k'k'k'k'k’k 

*  SCREEN  TEN  * 

*  INITIAL  CONDITIONS  * 

*  REACH  1  * 


ENTER  INITIAL  pH  FOR  THIS  REACH  7. 

ENTER  INITIAL  TEMPERATURE  FOR  THIS  REACH  78.6 


Figure  B22.  Initial  conditions,  Part  1 
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* 

* 

* 


SCREEN  TEN 
INITIAL  CONDITIONS 
REACH  1 


* 

* 

★ 


CONSTITUENT 

VALUE 

TEMPERATURE  (F) 

78,600 

DISSOLVED  OXYGEN  (mg/1) 

7.000 

CBOD5  (mg/1) 

1.000 

IRON  (mg/1) 

1.000 

MANGANESE  (mg/1) 

1.200 

SULFIDE  (mg/1) 

0.000 

SULFATE  (mg/1) 

0.000 

pH 

7.000 

CHL  ALGAE  (mg/1) 

0.000 

CONSTITUENT 

VALUE 

ORGANIC  N  (mg/1) 

0.320 

AMMONIA  N  (mg/1) 

0.280 

NITRITE  N  (mg/1) 

0.000 

NITRATE  N  (mg/1) 

0.030 

ORGANIC  P  (mg/1) 

0.000 

DISSOLVED  P  (mg/1) 

0.000 

ANC  (mg/1) 

0.000 

CONSERVATIVE  1  (mg/1) 

0.000 

CONSERVATIVE  2  (mg/1) 

0.000 

•  DO  YOU  NEED  TO  MAKE  ANY  CHANGES  (Y/N)?  N 

•  « 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAi. 


Figure  B23.  Initial  conditions.  Part  2 
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42.  Screen  11  (incremental  inflows)  is  just  like  screen  10  in  that  it 
appears  in  two  parts.  The  first  of  the  screens  contains  information 
describing  incremental  inflows  (Figure  B24)  and  asks  if  there  are  any  in  this 
reach.  If  there  are  no  incremental  inflows  in  a  reach,  the  user  answers  no 
(N) .  If  yes  (Y)  is  selected,  the  second  screen  will  appear,  and  constituent 
values  for  incremental  inflow  can  be  entered  (Figure  B25). 


*irkirkirkirkirk*irk*1c****irkicicic**it'k*1cit*-k*-kit*irkilrk*’k**icilrk**** 


* 

* 

* 

* 


SCREEN  ELEVEN 
INCREMENTAL  INFLOW 
REACH  1 


* 

* 

* 

* 


INCREMENTAL  INFLOWS  ARE  THOSE  DUE  TO  NONPOINT  SOURCES 
SUCH  AS  RUNOFF  FROM  A  HEAVY  RAINFALL.  SUCH  FLOWS  ARE 
ASSUMED  TO  BE  EVENLY  DISTRIBUTED  OVER  EACH  ELEMENT  IN 
THE  REACH. 

ARE  THERE  INCREMENTAL  INFLOWS  IN  THIS  REACH  (Y/N)? 


AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAi 

Figure  B24.  Incremental  inflows.  Part  1 
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* 

•k 

* 

SCREEN  ELEVEN 

* 

* 

INCREMENTAL  INFLOW  1 

* 

* 

REACH  1 

* 

DATA  TYPE 

DEFAULT 

DATA  TYPE 

DEFAULT 

INFLOW  (cms) 

0.000 

AMMONIA  N  (mg/1) 

0.000 

TEMPERATURE  (C) 

20.000 

NITRITE  N  (mg/1) 

0.000 

DISSOLVED  OXYGEN  (mg/1) 

0.000 

NITRATE  N  (mg/1) 

0.000 

BOD  (mg/1) 

0.000 

ORGANIC  P  (mg/1) 

0.000 

pH 

7.000 

PHOSPHATE  P  (mg/1) 

0.000 

SULFATE  (mg/1) 

0.000 

CONSERVATIVE  ELEMENT  1 

(mg/l) 

0.000 

CHLOROPHYL  (mg/1) 

0.000 

CONSERVATIVE  ELEMENT  2 

(m^l) 

0.000 

ORGANIC  N  (mg/1) 

0.000 

ANC  (mg/1) 

0.000 

•  ANY  CHANGES  TO  THESE  VALUES  (Y/N)7 

•  • 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA^ 

Figure  B25.  Incremental  inflows.  Part  2 
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43.  After  finishing  with  screen  11,  TWQI  returns  to  screen  5  (Figure 
B15)  if  more  than  one  reach  is  being  modeled.  Screens  5  through  11  will  be 
repeated  once  for  each  additional  reach.  After  finishing  the  last  reach, 
screens  12 -P  (Figure  B26)  and  12 -W  (Figure  B27)  will  appear  only  if  there  are 
point  loads  or  withdrawals .  These  screens  will  appear  in  the  order  in  which 
the  point  loads  and  withdrawals  occur  downstream  from  the  headwater.  The 
screen  for  point  loads  (Figure  B26) ,  12-P,  is  very  similar  to  that  for  the 
headwater  data  (Figure  B14) .  It  appears  once  for  each  point  load,  and  concen¬ 
tration  values  for  all  water  quality  constituents  selected  to  be  modeled  can 
be  entered. 

•  *******-k**************************-k**iHe************-kit  • 

*  *  *  • 

°  *  SCREEN  TWELVE  -  P  * 

•  *  DATA  TYPE  11:  * 

•  *  POINT  LOAD  BOUNDARY  CONDITIONS  *  • 

•  *  REACH  1  * 

•  ****************irkirk*****-k**************-k************  • 


POINT  LOAD 

point  2 

CONSTITUENT 

VALUE 

CONSTITUENT 

VALUE 

TEMPERATURE  (C) 

20.000 

ORGANIC  N  (mg/1) 

0.000 

DISSOLVED  OXYGEN 

(mg/1) 

7.000 

AMMONIA  N  (mg/1) 

0.000 

CB005  (mg/1) 

1.000 

NITRITE  N  (mg/1) 

0.000 

IRON  (mg/1) 

0.100 

NITRATE  N  (mg/1) 

0.000 

MANGANESE  (rag/1) 

0.100 

ORGANIC  P  (mg/1) 

0.000 

SULFIDE  (mg/1) 

0.000 

DISSOLVED  P  (mg/1) 

0.000 

SULFATE  (mg/1) 

0.000 

ANC  (mg/1) 

0.000 

pH 

7.000 

CONSERVATIVE  1  (mg/1) 

0.000 

CHL  ALGAE  (mg/1) 

FLOW 

0.000 

(cms) 

CONSERVATIVE  2  (mg/1) 
1.500 

0.000 

•  ANY  CHANGES  (Y/N)? 

44444444444444444444444444444444444444444444444444444444444444444444444444444441 
Figure  B26.  Point  load  boundary  conditions  (screen  12-P) 
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44.  The  only  information  required  on  screen  12 -W  is  the  withdrawal 
flow.  This  value  must  be  entered  as  a  negative  real  number. 

***************************************************** 

*  * 

*  SCREEN  TWELVE  -  W  * 

*  DATA  TYPE  11:  * 

*  WITHDRAWL  BOUNDARY  CONDITIONS  * 

*  REACH  1  * 

***************************************************** 


NAME  OF  WITHDRAWL 
WITHDRAWL(-)  (cms) 

ARE  THERE  ANY  CHANGES  (Y/N)? 


withone 

-2.000 


44444444444444444444444444444444444444444444444444444444444444444444444444444441 

Figure  B27.  Withdrawal  flow  (screen  12-W) 
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A5.  A  screen  containing  local  climatological  data  will  appear  only  when 
temperature  is  being  modeled  (Figure  B28).  The  cloudiness  scale  has  a  range 
from  0  to  1  with  a  value  of  0  indicating  clear  skies.  The  other  meteorologi¬ 
cal  data  required  are  dry  bulb  temperature,  wet  bulb  temperature,  barometric 
pressure,  and  wind  speed. 

OdddddddddddddddddddddddddddddAddddddddddddddddddddddddddddddddddddddddddddddddC 

•  o 

•  'kitit'k'k^^iricieirkieie^icicicicicicicic'k'k'kic^^ieirk-kicikieie'k-k^icieit'k-k'k'k'kie^'kiirk^'k'k-k-k-kie  * 

•  *  *  • 

*  *  LOCAL  CLIMATOLOGICAL  DATA  * 

•  *  *  • 

•  ************************************************************  • 


DATA  TYPE: LOCAL  CLIMATOLOGY 


CLOUDINESS  (O.-l.)  0.00 
DRY  BULB  TEMPERATURE  (C)  20.00 
WET  BULB  TEMPERATURE  (C)  20.00 
BARAMETRIC  PRESSURE  (mb)  980.00 
WIND  SPEED  (m/s)  10.00 


DO  YOU  NEED  TO  MAKE  ANY  CHANGES  (Y/N)? 


Figure  B28.  Local  climatological  data 
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46.  The  next  screen  deals  with  observed  data  (Figure  B29)  at  stations 
in  the  study  reach.  This  is  optional,  and  the  data  entered  are  used  for 
graphing  only.  If  there  are  no  observed  data,  the  user  enters  no  (N) .  If 
there  are  observed  data,  the  user  enters  yes  (T)  and  then  enters  the  number  of 
observed  data  sets.  For  each  observed  data  set,  the  observed  data  screen  will 
appear  (Figure  B30) .  This  screen  is  similar  to  the  screen  for  point  loads 
(Figure  B26)  except  that  in  place  of  flow,  the  river  kilometer/mile  of  the 
sampling  location  is  entered.  All  data  entered  with  this  screen  are  stored  in 
the  file  named  OBSERVED.DAT. 


DO  YOU  HAVE  ANY  OBSERVED  DATA?  (Y/N) 


•  • 

•  •k-kirkltlrklrklcirkltltltlililtltlrlclfklileltlfkieiflrk-kifklc-k-k-kitie'irk-k-kifklrk-klrkiclc-kifkirtrk  * 


* 

* 

OBSERVED 

DATA 

* 

* 

a. 

1rk1t1r1i1rkirkitirkirkit1rkirkit1eit-lt-k-k1c-k1t1t1t1c1rk'irkitirkirk-k-irkic1rk-k1Hrirk'kieit'icirirkitie 

CONSTITUENT 

VALUE 

CONSTITUENT 

VALUE 

TEMPERATURE  (F) 

78.620 

ORGANIC  N  (mg/1) 

0.320 

DISSOLVED  OXYGEN  (mg/1) 

7.000 

AMMONIA  N  (mg/1) 

0.280 

CBODS  (mg/1) 

0.000 

NITRITE  N  (mg/1) 

0.000 

IRON  (mg/1) 

1.000 

NITRATE  N  (mg/1) 

0.030 

MANGANESE  (mg/1) 

1.200 

ORGANIC  P  (mg/1) 

0.000 

SULFIDE  (mg/1) 

0.000 

DISSOLVED  P  (mg/1) 

0.000 

SULFATE  (mg/1) 

0.000 

ANC  (mg/1) 

0.000 

pH 

6.500 

CONSERVATIVE  1  (mg/1) 

0.000 

CHL  ALGAE  (mg/1) 

0.000 

CONSERVATIVE  2  (mg/1) 

0.000 

RIVER  MILE  OF  OBSERVED  DATA  62.60 


•  DO  YOU  NEED  TO  MAKE  ANY  CHANGES  (Y/N)? 

•  • 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAi. 


Figure  B30.  Observed  data,  input  screen 
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47.  The  final  screen  Indicates  that  TWQI  was  completed  successfully  and 
displays  the  name  of  the  data  file  that  has  just  been  created  (Figure  B31) . 

044444d444d4d444d4444d4d4ddddddddd4dd4d444d4d44444Fr 1  12 -  21 • 199044410 : 48  4C 

•  C:\TAILW 


PROGRAM  TWQI  IS  COMPLETE 
NEW  DATA  SET  IS:barry.91 


•  >  • 

44444444444444444444444444444444444444444444444444444444444444444444444444444441 

Figure  B31.  TWQI,  final  screen 
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Screens  in  TFIX.FOR 

48 .  TFIX  allows  the  user  to  make  changes  to  an  existing  TAILWTR  input 
data  file.  Only  certain  information  entered  using  TWQI  may  be  altered  with 
TFIX.  Data  that  cannot  be  changed  using  TFIX  are  variables  such  as  the  number 
of  reaches  and  number  of  point  loads ,  because  new  lines  of  data  would  have  to 
be  entered  into  the  input  data  set.  TFIX  does  not  allow  lines  to  be  added, 
only  changed.  Items  such  as  headwater,  initital,  or  point  load  conditions, 
simulation  date,  and  theta  values  can  be  changed  in  TFIX.  If  the  user  must 
add  a  new  reach  or  point  load,  then  TWQI  should  be  run  or  the  new  information 
manually  entered  into  the  input  data  set. 

49.  To  alter  an  existing  data  file  using  TFIX  instead  of  creating  a  new 
one  with  TWQI,  the  user  types  TWQB  OLD  when  beginning  the  batch  program.  This 
causes  TWQB  to  skip  TEXPLAIN,  TMENU,  TWQI,  and  proceed  directly  to  TFIX.  Once 
TFIX  is  exited,  TAILWTR  is  run  automatically  along  with  the  graphics  options. 
If  OLD  is  not  specified  when  TWQB  is  begun,  then  TFIX  is  not  run. 

50.  The  first  input  for  TFIX  are  the  names  of  the  existing  data  file 
and  the  name  of  the  new  file  to  be  created.  Data  from  the  existing  file  are 
read,  and  can  be  altered  if  necessary.  These  data  are  then  stored  in  the  new 
data  file. 
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51.  TFIX  next  displays  a  screen  that  contains  nine  options  (Figure 
B32) .  The  first  eight  options  correspond  to  different  portions  of  the  TAILWTR 
input  data  file  that  are  to  be  changed.  The  ninth  option  exits  TFIX.  To  make 
a  change,  the  user  simply  selects  the  appropriate  option.  When  more  than  one 
option  is  to  be  selected  during  a  session,  the  user  selects  the  option  with 
the  lowest  number  FIRST.  For  example,  suppose  the  user  must  make  changes  to 
some  of  the  theta  values  (option  3),  headwater  data  (option  6),  and  hydraulic 
data  (option  4).  The  user  should  choose  option  3  first,  then  option  4,  and 
finally  option  6 .  This  is  required  because  of  the  method  TFIX  uses  to  create 
a  new  data  file.  It  reads  the  old  data  file  sequentially  and  copies  the 
unchanged  portions  line  by  line  to  the  new  data  file.  The  first  options  in 
the  list  correspond  to  data  stored  at  the  beginning  of  the  data  file.  These 
data  must  be  written  to  the  new  data  file  before  data  corresponding  to  subse¬ 
quent  options,  so  that  the  structure  of  the  data  file  is  maintained. 


CHOICE  OF  PLACES  TO  EDIT 


1  -  Title  Cards 

2  •  Control  Data 

3  -  Theta  Values 

4  -  Hydraulic  data 

5  -  Biological  and  Chemical  Coefficients 

6  •  Headwater  data 

7  -  Point  Load  data 

8  -  Local  Cllmatalogical  data 

9  -  No  changes 


•  Enter  Choice  (INTEGER) 

•  • 

•  DO  NOT  REPEAT  A  NUMBER  OR  SELECT  ONE  SMALLER 

•  THAN  ANY  NUMBER  PREVIOUSLY  SELECTED 

«  • 

AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAi. 

Figure  B32.  TFIX  options  screen 
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52.  Each  of  the  first  eight  options  has  its  own  screen.  The  cursor 
will  go  to  each  item  on  the  screen.  The  user  makes  changes  by  entering  the 
new  value  when  the  cursor  is  in  the  appropriate  place.  To  accept  the  value 
shown,  the  user  presses  the  ENTER  key.  None  of  these  screens  allows  the  user 
to  go  back  and  make  changes  if  mistakes  were  made  entering  new  values .  TFIX 
will  have  to  be  rerun  to  correct  a  mistake. 


B44 


Screens  In  TAILWTR.FOR 

53.  TAILWTR  requires  only  the  name  of  the  input  data  file  and  the  name 
of  the  output  file.  The  input  data  file  will  usually  be  the  data  file  just 
created  by  TUQI.  However,  the  name  of  any  old  input  data  file  can  be  entered. 
There  are  three  options  concerning  output.  It  can  be  written  to  the  screen, 
written  directly  to  a  printer,  or  written  to  a  file.  The  final  option  is 
suggested  due  to  the  volume  of  output  generated.  To  have  the  output  written 
to  a  file,  enter  the  name  of  the  output  file  at  the  prompt  concerning  the 
output  device.  The  name  of  the  file  can  be  up  to  eight  characters  long  with 
an  extension  of  three  characters. 

54.  When  running,  TAILWTR  will  print  messages  to  the  screen  indicating 
what  tasks  are  being  performed.  If  desired,  these  messages  can  be  turned  off 
by  answering  no  (N)  when  asked  about  the  run  time  messages. 


Plotting  and  Graphing  Options 


55.  Graphical  routines  have  been  developed  for  TWQM  output  using 
GRAPHER  and  TPLOT  software.  A  copy  of  one  of  these  is  required  in  order  to 
use  these  graphics  routines.  After  TAILWTR  has  finished,  a  screen  will  appear 
listing  the  graphing  and  plotting  options  along  with  options  to  exit  TWQM  or 
return  to  its  beginning  (Figure  B33).  Viewing  (Figure  B33)  a  graph  means  that 
the  graph  will  appear  on  the  screen,  and  plotting  (Figure  B33)  a  file  means 
that  the  graph  can  be  plotted  using  an  HP  Series  Il-compatible  laser  jet 
printer.  If  one  of  the  options  corresponding  to  viewing  or  plotting  is 
selected,  a  screen  containing  a  list  of  the  different  constituents  that  can  be 
viewed  or  plotted  will  appear  (Figure  B34) .  To  plot  or  view  graphs  of  TAILWTR 
output  data,  the  user  selects  the  appropriate  option. 
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Figure  B33.  First  graphics  screen 
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Which  would  you  like  to  Graph  (A,  B  .  .  .  H  or  EXIT)  ? 

Figure  B34.  Constituent  options 


56.  When  the  constituent  options  menu  is  exited,  the  first  graphics 
screen  appears  and  another  selection  can  be  made  (Figure  B33) .  The  user 
selects  E  to  exit  from  TWQM. 

57.  The  plotting  options  for  TPLOT  may  be  selected  even  if  a  printer  is 
not  attached  to  the  computer.  This  program  creates  the  plotting  files,  which 
are  dumped  to  a  printer.  One  file  is  created  for  each  graph.  These  files 
are:  TPTPLOTl . OUT ,  TPTPLOBS . OUT ,  FETPLOTl . OUT .  MNTPLOTl . OUT ,  DOTPLOTl . OUT , 
FETPLOBS.OUT,  MNTPLOBS . OUT ,  and  DOTPLOBS . OUT .  When  the  plotting  option  for 
TPLOT  is  selected,  these  files  are  created  but  not  sent  to  a  printer.  Once 
TWQB  is  complete,  these  files  can  be  sent  to  a  printer  by  typing  HP2LJET  fn.ft 
if  an  HP-compatible  laser  jet  II  printer  is  being  used  (Benton  1987).  Other 
files  are  included  with  TPLOT  that  allow  for  other  printers  to  be  used.  For 
more  information,  refer  to  the  TPLOT  user's  manual.^ 
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APPENDIX  C:  KINETIC  RATE  REGRESSIONS 
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Figure  C2.  Regressions  for  ammonia  nitrogen  removal  rate 
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DOi/2 
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D„ 
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Dt 

D° 

D* 

[Fe(II)] 

(FeS) 

FRACT 
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kpss 

K 


Regression  coefficient,  dimensionless 
Coefficient  of  flow  for  K2 
Bottom  planar  area 

Surface  area  of  the  bed  control  volume ,  m^ 

Exponent  of  flow  for  K2 

Escape  coefficient,  1/ft 

Constituent  concentration 

Mn'*’^  observed  at  downstream  station,  mg/i 

Concentration  of  reduced  manganese  or  iron  sulfide  in  the  bed,  mg/i 

Initial  (release)  Mn'*'^  concentration,  mg/i 

Concentration  of  H2S  in  the  stream 

Average  stream  depth,  ft 

Longitudinal  dispersion 

Dissolved  oxygen  half -saturation  constant,  mg/i 

Molecular  diffusion  coefficient  of  Mn'*'^  at  25  ®C, 

6.88  X  10'®  cmVsec 

Molecular  diffusion  coefficient,  ftV^ay 

Numerical  dispersion  for  steady-state  conditions,  L^/T 

Stream  dispersion,  L^/T 

Molecular  diffusion  coefficient  of  Mn"^^  at  stream  temperature 
T  ,  cmVsec 

Molecular  diffusivity  of  DO,  cmVsec 
Molecular  diffusivity  of  H2S,  cmVsec 
Concentration  of  ferrous  iron,  moles/i 
Iron  sulfide  concentration,  mg/i 
Fraction  of  total  dissolved  sulfide  as  H2S 
Gravitational  acceleration,  L/T^ 

Hydraulic  depth,  ft/ft 

Hydrogen  ion  concentration,  moles/i 

Volatilization,  mg 

Heterogeneous  reation  rate,  day"^  (moles/i)'* 

3.0  X  10'^^  moles/i  min'^  at  20  ®C 
First-order  sulfide  oxidation  rate,  day'^ 

Sulfide  oxidation  rate,  day'^ 

Oxidation  rate  at  20  ®C 
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Sc 
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Homogeneous  reaction  rate,  day"^  (moles/i)"^ 

Reaction  coefficient 
Nitrification  rate,  day'^ 

Oxidation  rate  of  reduced  manganese  or  iron  sulfide  in  the 
bed,  1/day 

Overall  reation  rate  for  ferrous  iron 
First-order  removal  rate  for  Mn'*’^ 

Mn'*’^  loss  rate  for  a  specific  field  condition,  1/day 
Volatilization  rate  of  H2S,  day*^ 

Reaction  rate  at  temperature  T 
DO  stream  reaeration  rate,  day'^ 

Reaction  rate  at  20  “C 

Dissolved  (reduced)  manganese  concentration,  mg/i 

Reduced  manganese  concentration  in  sediment  pore  water,  mg/i 

Dissolved  (reduced)  manganese  concentration,  moles/i 

Particulate  (oxidized)  manganese  concentration,  moles/i 

Molecular  weight  of  DO,  g/mole 

Molecular  weight  of  HjS,  g/mole 

Manning's  coefficient 

Hydroxide  ion  concentration,  moles/i 

Point  load 

Partial  pressure  of  oxygen,  atm 
Flow  rate 

Water  surface  or  streambed  slope,  nondimens ional 

Rate  of  change  in  concentration  resulting  from  transformation 
or  chemical  reactions 

Sulfide  concentration  in  the  form  of  HS',  mg/i 

Schmidt  number  (j//D)  ,  nondimeni  ional 

Slope  of  energy  gradient 

Slope  of  streambed,  ft/ft 

Travel  time  to  station,  days 

Ambient  stream  temperature,  “C 

Flow  time  through  reach,  days 

Average  stream  velocity,  m/sec 

Shear  velocity  (/gHs) ,  L/T 

Stream  mean  velocity 

Control  volume  of  fixed  volume 
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Vj,  Volume  of  bed  control  volume,  m® 

V,  Hass  transfer  velocity  across  diffusive  sublayer  separating 

water  and  substrate,  L/T 

Measured  mass  transfer  rate,  L/T 

V,c  Mass  transfer  coefficient  computed  with  Equation  25,  L/T 

w  Loading  rate  (i.e.,  mass  transfer  rate  or  settling  rate)  to  bed, 

m/day 

U  Top  width  of  stream,  m 

Stream  withdrawal  in  TWQM  input 
Wg  Settling  rate  of  FeS,  m/day 

X  Distance  downstream 

X  pH  of  stream 

a  Statistical  error  bound  in  analytical  solution  of  1-D  transport 

Stoichiometric  conversion  factors 
Ah  Change  in  water  surface  elevation  in  reach,  ft 

ASOD  Additional  SOD  resulting  from  manganese  or  iron  sulfide  oxidation  in 

bed,  g  Oa/CmVday) 

Ax  Computational  element  length,  m 

0  Temperature  correction  constant 

1/25  Kinematic  visccosity  of  water  at  25  “C,  6.88  x  10*®  cmVsec 

j/j  Kinematic  viscosity  at  stream  temperature  T  ,  cm^/sec 
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